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ABSTRACT 


The  purpose  of  this  research  to  stxidy  the  phenomena  -of 

multiple  quantum  processes  and  to  consider  possible  applications.  These 
processes  Involve  the  Interaction  of  more  than  a  single  quantum  of  radia¬ 
tion  vlth  an  atomic  system  and  are  nonlinear  In  character.  An  equation 
of  motion  approach  Is  used  In  the  solution  of  the  quantum, mechanical 
problem,  allowing  the  calculation  of  the  observable  quantities  important 
for  the  radiation  processes  considered.  _ _ - — ■ — — 


1,  third 
The 


ThVee  specific  cases  have^en  cbnsldered.  :pa  the  first 
harmonic  u^ene:^ation  is  predated  f^or  a  tko-level  quantum 
magniyfle  of  tlieSnonllnM^ty  and  the  magnti^ude  of  the  thl^d  harmo] 

generated  arb  ej^tiated  in  terms  of  knc^  parametoi^s  of  the  qubptum 
and  Its  surpo^dlngs.  A  shl|;t  of  the  nanural,^anBltlon  freque: 

'of  the  quantum  aj^stem  is  also  predicted.  In  the  ^fecond  case,  the  same 
two-level  s^xem  is  founc^to  ,glfve  rise  to  a  of  p^ajjetric  effect, 
sultable^fbr  aii5)llflcatloj)^^The  threshold^ondltions  pot  parametric 
oscill/rcion  are  found  iff,  ten^  of  the  kl!^own  properties  oK^the  system. 

leXhlrd  case  considered  was  that  Bff  a  three-level  systemA  Its  properties 
aS\a  frequency  c^verter  were  potfif^ed  out  and  a  detailed  studKof  second 
harmhijlc  genq^tion  was  perfpraed,  w\th  the  nonlinearity  and  sedtod 
harmonfc  Pi^er  evaluated  terms  of  the  Imcwn  properties  of  the  s^^i^m. 

Of  pajrti^lar  interestydas  the  predlctlOTV.  th^'  the  optimum  hanajjHic 
genejjdCtion  O^curredy^Jhen  two  of  the  thre^^i^sitions  involAfsdf*"dif f ere 


degenere 


one  linewldth. 


A 


Two  expe^^f^^s  were  perfomstfl  In  the  fi^^,  third  harmonic 
jenerationjjllenomen^^  predicte^f^’’  the  theory  was  ^oAerved.  T 

ory.  Ii 


the  esi^rlment  werfe  In^^^^^r  agreement  with 
’l^f^nt  the  second  h^inonlc  generation  p) 


thi 


omenon  prei 


thea^vwas  observed/dlnd  wa^  found  to  be  IfT  general  agreei 
within  ^e  experimental  uncertainty,  j,' 


From  the  resxilts  of  this  research  it  Is  thus  concluded  that  a 
quantum  mechanical  system^  possesses  nonlinear  as  veil  as  linear  properties 
and  that  these  nonlinear  properties  may  find  application,  especially  In 
the  sub-mllllmeter  and  optical  regions  vhere  suitable  nonlinear  elements 


do  not  presently  exist 


N 
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IBJPBQDOCTIOH 


A.  BACaSOBOOND 

In  1955  Gordon,  Zelger  and  Townes^  announced  the  successful  operation 

of  the  amnonla  beam  maser.  With  this  pioneering  breakthrough  the  field 

of  quantum  electronics  had  its  beginning.  The  achievements  in  this  field 

have  been  many,  starting  with  this  first  ammonia  maser  and  inclxuilng  the 

2 

solid  state  microwave  maser,  proposed  by  Bloembergen  and  successfullj'^ 

3 

operated  by  Scovil,  Peher,  and  Seidel,  the  optical  maser  proposed  by 

k  5  6 

Schawlow  and  Townes  and  operated  by  Malman  in  a  solid,  and  by  Javan 

in  a  gas,  and  most  recently,  the  successful  operation  of  a  semiconductor 
7-9 

Junction  laser. 

The  field  of  quantum  electronics  differs  from  its  predecessor,  classical 
electronics,  in  that  it  makes  use  of  the  internal,  quantized  energy  states 
of  an  atom  or  molecule  rather  than  the  translations^,  states  of  electrons. 

In  genersd,  efforts  to  date  in  the  application  of  these  properties  have 
been  concerned  with  producing  amplification  and  oscillation  using  the 
phenomenon  of  stimulated  emission  from  an  atcmilc  system  which  is  in  an 
Inverted  popiilation  state. 

A  st;idy  of  the  interaction  of  radiation  with  matter  forms  the  basis 
for  the  analysis  of  these  qviant\mi  electronic  systems.  Such  a  study  must 
be  quantum  mechanical  in  nature  in  order  to  account  for  the  discrete  nature 
of  the  atomic  system.  The  phenomenon  of  amplification  by  such  a  quantum 
system  is  a  linear  process  and  can  be  described  as  a  preponderance  of  stimu¬ 
lated  emission  over  absorption.  Such  a  linear  process  is  microscopically 
descrlbable  as  the  interaction  of  a  single  photon  with  the  atonic  system 
while  macroscoplcally,  a  ccmiplex  susceptibility,  independent  of  the  field 
strength,  may  be  used.  The  latter  description  gives  a  complex  polarization 
or  magnetization  proportional  to  the  applied  field. 

^^^The  words  atom  or  atomic  system  are  used  to  mean  atom  or  molecule. 
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A.  BACXCSiOUIlD 

In  1955  Gordon,  Zelger  and  Townee^  announced  the  successful  operation 

of  the  amnonla  team  maser.  With  this  pioneering  breakthrovigh  the  field 

of  quantum  electronics  had  Its  beginning.  The  achievements  In  this  field 

have  been  many,  starting  with  this  first  ammonia  maser  and  Including  the 

2 

solid  state  mlcrovave  maser,  proposed  by  Bloemhergen  and  successfully 

3 

operated  by  Scovil,  Feher,  and  Seidel,  the  optical  maser  proposed  by 

ij.  5  6 

Schavlov  and  Townes  and  operated  by  Maiman  In  a  solid,  and  by  Javan 

in  a  gas,  and  most  recently,  the  successful  operation  of  a  semiconductor 
7-9 

Junction  laser. 

The  field  of  quantvnn  electronics  differs  from  its  predecessor,  classical 
electronics,  in  that  it  makes  use  of  the  internal,  quantized  energy  states 
of  an  atom  or  molecule  rather  than  the  translational  states  of  electrons. 

In  general,  efforts  to  date  in  the  application  of  these  properties  have 
been  concerned  with  producing  amplification  and  oscillation  using  the 
phenomenon  of  stimulated  emission  from  an  atomic  system  which  Is  In  an  . 
inverted  population  state. 

A  study  of  the  interaction  of  radiation  with  matter  forms  the  basis 
for  the  analysis  of  these  quantum  electronic  systems.  Such  a  study  must 
be  quantum  mechanical  In  nature  in  order  to  account  for  the  discrete  nature 
of  the  atomic  system.  The  phenomenon  of  amplification  by  such  a  quantum 
system  is  a  linear  process  and  can  be  described  as  a  preponderance  of  stimu¬ 
lated  emission  over  absorption.  Such  a  linear  process  Is  microscopically 
descrlbahle  as  the  Interaction  of  a  single  photon  with  the  atomic  system 
while  macroscoplcsJly,  a  complex  susceptibility,  independent  of  the  field 
strength,  may  be  used.  The  latter  description  gives  a  complex  polarization 
or  magnetization  proportions^,  to  the  applied  field. 
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Itonllnear  Effects 

In  addition  to  the  linear,  single  quantiam  processes  important  for 
maser  theory,  there  exist  higher  order  processes  Involving,  on  a  mlcrosc  jIc 
scale,  the  interaction  of  more  than  a  single  quantum  of  radiation  with  the 
atomic  system.  Macroscoplcally  these  interactions  exe  found  to  be  described 
by  a  polarization  or  magnetization  which  is  dependent  on  q\Mdratlc  or  higher 
powers  of  the  applied  radiation  field.  This  nonlineeur  dependence  on  the 
field  makes  possible  such  phenomena  as  harmonic  generation,  parametric 
amplification,  limiting,  modulation,  demodulation,  and  other  frequency 
mifclng  effects. 

The  research  described  in  this  report  is  concerned  with  a  study  of 
these  nonlinear  phenomena.  Because  of  their  nonlinear  character  and  quai^tum 
mechanical  aspects  they,  will  be  called  nonlinear  quantum  effects.  They  may 
be  further  classified  as  bulk  nonlinear  effects  since  the  phenomena  are 
not  connected  with  either  a  surface  or  a  Junction  as,  for  example,  in  a 
semiconductor  diode. 

Although  nonlinear  phenomena  have  found  general  application  in 
various  fields  of  science  and  engineering,  little,  if  any,  use  has  been 
made  of  the  bulk  nonlinear  phenomena  considered  here.  In  the  region  of 
the  spectrum  from  dc  to  microwaves  the  nonlinear  properties  of  diodes, 
varactors,  and  electron  beams  have  had  the  most  applications.  This 
results  both  from  the  fact  that  such  systems  have  strong,  relatively 
frequency- Insensitive  nonlinearities  and  that  these  nonlinearities  are 
esisily  descrlbable  in  terms  of  familieur  concepts. 

The  bulk  nonlinearities  of  the  type  considered  here,  on  the  other 
hand,  are  generally  smaller  in  magnitude,  frequency-sensitive,  and  are 
described  in  a  less  familiar  manner  than  their  classical  counterparts. 

With  the  host  of  presently  available,  strong,  nonlinear  elements,  the 
question  might  be  asked  as  to  the  purpose  of  studying  these  nonlinear 
phenomena;  Two  main  reasons  can  be  offered.  First  of  all,  such  a  study 
would  contribute  to  a  further  understsmding  of  phenomena  associated  with 
the  interaction  of  radiation  with  matter.  Of  particular  interest  is  the 
study  of  radiation  processes  where  atomic  coherence  effects,  imposed  by 
strong  fields,  beccane  important.  A  second  reeison  comes  from  the  point  of 
view  of  applications.  Although  the  nonlinear  properties  of  the  presently 
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used  diode,  etc.,  are  relatively  frequency-independent  below  the  microwave 
range,  their  efficiencies  do  fall  off  with  increasing  frequency  above  this 
range.  Thus  applications  in  the  submillimeter  to  optical  frequency  ranges 
are  msirglnal,  if  possible  at  all. 

The  bulk  nonlinear  properties  described  in  this  report  are  found  to 
depend  on  the  natural  transition  frequencies  of  the  atom.  Since  atoms 
exhibit  resonances  throughout  the  electromagnetic  spectrum  these  nonlinear 
properties  will  also  occur  throughout  the  spectnom.  Further,  it  is  shown 
in  this  report  that  the  magnitude  of  the  nonlinearity  is,  among  other 
things,  dependent  on  the  Boltzmann  factor  and  Increases  in  magnitude  with 
increasing  frequencies.  Thus  at  high  frequencies  these  nonlinear  effects 
may  be  expected  to  play  an  Important,  if  not  dominant,  role. 

Scope  of  the  Problem 

The  purpose  of  this  reaearch  is  to  study  these  nonlinear  effects  and 
to  consider  possible  applications. 

The  first  task  in  a  study  of  these  nonlinear  effects  is  to  find  a 
method  of  formulation  of  the  problem.  Such  a  formulation  should  be  valid 
for  the  large  field  strengths  encountered  in  exciting  the  nonlinearities 
and  must  be  capable  of  predicting  coherent  emission  as  well  as  absorption 
for  the  higher  order  processes.  For  such  an  analysis  the  usual  approach 
of  time-proportional  transition  probabilities  and  instantaneous  quantum 
Jumps  between  unperturbed  energy  states,  used  in  solving  maser  problems, 
is  not  applicable  inasmuch  as  the  assumptions  implicit  in  the  derivation 
of  the  method  are  not  valid.  In  perticuleur,  the  assumption  of  weak  field 
strengths  used  in  deriving  the  time-proportional  transition  probabilities 
is  not  valid. 

From  such  an  analysis  we  sho\ild  like  to  be  able  to  evalmte  the  magnitude 
of  the  nonlinear  effect  in  terms  of  known  or  measurable  properties  of  the 
system  such  as  natviral  frequencies,  linewidths,  and  matrix  elements.  The 
effects  of  temperature,  concentration  and  pop\alation  should  also  be 
considered. 

Finally,  possible  applications  of  these  nonlinear  phenomena  should 
be  considered.  Included  among  these  considerations  should  be  the  effects 
of  the  material's  surroundings,  whether  a  cavity  resonator,  propagating 
circuit,  or  free  space. 
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Reeulta 

The  results  pf  the  research  described  in  this  report  may  be  summarized 
as  follows: 

1.  A  study  lias  been  made  of  the  Interaction  of  a  radiation  field  with 
an  atomic  system  using  an  equation  of  motion  formulation  including  the  effects 
of  relaxation.  Tills  approach  directly  relates  macroscopic,  observable 
quantities  to  their  microscopic  sourqes,  predicting  nonlinear  effects  as 
well  as  linear  ones.  The  effects  of  atomic  coherence  are  included,  and 

the  formulation  Is  valid  for  lArge  as  well  as  small  field  strengths.  From 
this  ana'iysis  varj.ous  nonlinear  effects  may  be  predicted  and  their  magnitude 
may  be  calculated  in  terms  of  known  parameters  of  the  system. 

2.  Detailed  theoretical  analyses  have  been  carried  out  for  three 
specific  cases. 

a.  Harmonic  generation  in  a  two-level  system.  From  a  consideration 
of  the  interaction  of  an  applied  radiation  field  with  a  two-level  electric 
dipole  system,  third  harmonic  generation  is  predicted. Two  values  of 
the  frequency  of  the  applied  radiation  field  are  found  which  result  in  a 
strong  effect,  namely  when  applied  frequency  or  its  third  harmonic  is 
near  the  natural  transition  frequency  of  the  atcan. 

A  shift  of  the  natural,  trsuisition  frequency  due  to  the  strong  rf  fields 
is  predicted  from  the  theory.  This  shift  is  the  electric  dipole  eq\iivalent 
of  the  Bloch-Siegert,  shift  predicted  by  Bloch^^  for  the  case  of  magnetic 
resonsince. 

The  magnitude  of  the  induced  third  order  polarization  is  calculated 
in  terms  of  the  known  parameters  of  the  system  euid  the  power  generated  by 
such  a  means  is  calculated,  assimiing  a  resonant  cavity  structure.  The 
magnitude  of  the  nonlinearity  is  found  to  be  independent  of  whether  or  not 
the  system  is  in  an  inverted  population  state,  although  population  inversion 
is  found  to  alter  the  efficiency  throu^  a  modification  of  the  source 
impedance.  The  dependence  of  the  nonlinear  effect  on  other  properties  of 
the  atmlc  system  Is  also  considered. 

12 

b.  Parsmetric  effects  in  a  two-level  system.  Parsmetrlc  oscillation 
an^  amplification  are  considered  for  a  two-level  electric  dipole  system. 

The  frequency  conditions  for  such  an  effect  are  that  the  pump  frequency 

be  equal  to  one-half  the  sxmi  of  the  signal  plus  the  idler,  where  the  idler 
may  be  either  a  cavity  mode  or  the  Interpal  transition  of  the  molecule. 
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An  analysis  Is  carried  out  for  the  latter  case  and  the  start-osclllatlon 
condition  is  evaluated  in  terms  of  the  parameters  of  the  atomic  system 
and  the  surrounding  cavity. 

c.  Harmonic  generation  in  a  three-level  system.  A  three-level 

system  in  which  two  natural  frequencies  are  nearly  degenerate  is  found 

to  present  a  nonlinearity  suitable  for  second  harmonic  generation.  Second 

harmonic  generation  Is  predicted  and  Its  magnitude  is  evaluated  fran  the 

general  theory.  Of  particular  Interest  is  the  prediction  that  the  effect 

will  occur  not  when  the  transitions  are  degenerate  but  when  their  values 

differ  by  a  frequency  corresponding  to  one  llnewldth.  This  second  harmonic 

generation  effect  predicted  by  the  theory  was  subsequently  observed  by 
13 

Kelllngton  and  Independently  by  the  author. 

3.  Two  experiments  at  microwave  frequencies  have  been  performed  to 
verify  various  aspects  of  the  phenomena  predicted  by  the  theory. 

a.  Third  harmonic  generation  in  a  two-level  electric  dipole  system 

l4 

was  observed  using  ammonia  gas.  The  fundamental  and  harmonic  frequencies 
were  8.5  hMc  and  25.5  WUc,  respectively,  and  the  inversion  transition  at 
24  IsMc  was  the  two-level  system.  This  experiment  provided  a  direct  con¬ 
firmation  of  the  effect  predicted  by  the  theory  and  gave  results  in  general 
agreement  with  the  theory. 

b.  The  phenomenon  of  second  harmonic  generation  in  a  three-level 
system  was  observed  using  the  paramagnetic  resonance  levels  in  the  ground 
state  of  ruby.  This  effect  was  independently  observed  and  reported  by 

13 

Kellington,  who  also  used  ruby  under  approximately  the  same  operating 
conditions . 

B.  BRIEF  HiarORY  OF  THE  STOUY  OF  MULTIPLE  QUAHTOM  EFFECTS 

Until  recently,  interest  in  the  phenomenon  of  multiple  quantxmi 
transitions  has  centered  on  the  stxidy  of  these  effects  arising  in  various 
types  of  absorption  spectra.  The  first  published  reports  on  the  observa¬ 
tion  of  this  effect  came  from  the  work  of  Kastler's^^'^^  group,  resulting 
from  their  stvidy  of  the  hyperfine  spectrum  of  sodium.  In  their  experiments 
they  observed  the  appearance  of  additional  spectral  lines  not  predicted  by 
the  theory.  These  lines  were  present  in  the  spectnjm  at  high  rf  power 
levels  but  were  not  present  at  low  powers,  while  the  frequencies  of  these 
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new  lines  corresponded  to  the  aJLgehralc  mean  of  adjacent  predicted  lines . 

An  example  of  the  type  of  spectrum  observed  Is  shown  In  Fig.  1.1.  In  this 

case  the  selection  rule  for  the  hyperflne  transition  Is  An  *  ±  1  ,  where 

m  Is  the  magnetic  quantum  number.  With  further  study  they  found  that 

the  phenomenon  responsible  for  these  additional  lines  was  a  double  quajitum 

absorption  for  which  ^  =  2  .  This  absorption  process  occurred  when  the 

applied  frequency  was  equal  to  one-half  of  the  frequency  difference  between 

levels  differing  In  m  by  two  units.  Physically  this  process  corresponded 

to  an  absorption  of  two  photons  with  the  system  going  directly  from  a  level 

17 

m  to  a  level  m  +  2  .  Winter  showed  that  the  absorbed  power  for  this 

two-quantum  process  was  proportional  to  the  square  of  the  applied  power, 

explaining  the  absence  of  the  effect  at  low  power  levels. 

In  subsequent  reiiorts,  Margerie,  Brossel  euid  Winter  '  reported 

three  and  four  quantim  absorptions  as  well  as  multiple  quantum  absorptions 

involving  combinations  of  two  applied  frequencies. 

Other  experimental  results  of  absorption  studies  where  multiple  quantm 

19  20 

phenomena  were  observed  have  been  described  by  Kush,  Hughes  and  Geiger, 

21  22 

Wolga  and  others.  Sorokin,  et  al.,  found  both  double  and  triple  quantum 
absorptions  in  some  of  their  paramagnetic  resonance  work  on  solids.  In  all 
these  studies  the  additional  lines,  due  to  multiple  quantum  processes,  were 
present  only  at  high  rf  field  strengths  and  in  absorption  rather  than  emission 
spectra. 

Theoretical  studies  of  these  multiple  quantum  absorption  processes  have 

23  2h  25  26 

been  done  by  Salwen,  Hughes  and  Grabner,  and  Winter.  '  These  studies 

involve  the  use  of  time-dependent  perturbation  theory  ceurried  to  higher 

order  and  predict  multiple  quantum  absorptions.  None  of  these  analyses, 

however,  considers  nor  predicts  any  form  of  emission  related  to  or  resulting 

from  these  higher  order  absorptions.  It  is  precisely  this  latter  effect 

which  is  responsible  for  such  phenomena  as  harmonic  generation. 

Following  the  proposal  €Uid  successful  operation  of  the  solid  state 

27 

maser,  detailed  theoretical  studies  of  maser  action  were  made  by  Javan, 

28  29 

Clogston  and  Yatsiv.  ^  These  disciuslons  pointed  out  that  the  rate 
equation  approach,  used  in  earlier  analyses  of  the  maser,  failed  to 
describe  some  of  the  quantum  mechanlcstl  aspects  of  the  radiation  process . 

In  psirtlcular  they  foiond  that  under  some  circumstances  double-quantum 
processes  could  play  an  important  role  in  maser  operation.  Yajima  and 
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Magnetic  Field 


FIG.  1.1 — Besonance  curves  obtained  by  Oagnac  et  eO.. 

In  the  upper  curve,  for  low  power  levels,  the 
four  expected  ^  >  ±  1  transitions  are 
observed.  In  the  lower  curve,  for  high  power 
levels,  the  three  dai  =  ±  2  ,  two-quantum 

transitions  appear  as  the  very  sharp 
resonances. 


30-32 

Etoimoda  observed  maser  operation  where  these  doxible  quantum  processes 

33  3^  35 

play  the  dominant  role.  Anderson  and  Suhl  '  at  the  same  time  analysed 
the  operation  of  two  new  forms  of  quantum  mechanical  amplifier,  the  reaction 
field  amplifier  and  the  ferromagnetic  amplifier,  both  of  which  operated 


without  the  need  for  population  inversion. 


(1) 

.36 


iU-ong  this  same  line  a 


proposal  for  a  Roman  maser  was  made  by  Javan  and  was  experimentally 

37 

verified  by  Winjter.  The  operation  and  analysis  of  these  last  three 
devices  is  closely  related  to  the  subjects  discussed  in  this  report. 

Recently  there  has  been  renewed  interest  in  the  phenomenon  of  multiple 
quantum  processes,  this  time  in  the  optical  frequency  range.  Franken, 
et  al.,  generated  optical  hsurmonica  by  focusing  the  output  of  an  optical 


39 


40 


maser  on  a  piece  of  qiiartz.  Oiordmaine,  Maker,  et  al.,  and  Terhune, 
li-l 

et  al.,  followed  with  other  experimental  results  of  harmonic  generation. 
k2  ^3 

Kaiser  and  Garrett,  and  recently  Abella,  have  observed  optical  double¬ 
photon  absorption  processes.  Another  form  of  nonlinear  process,  coherent 

Raman  scattering  at  optical  frequencies,  has  been  observed  by  research 

45 

groups  at  Hijghes.  *  The  microwave  modulation  of  li^t  has  been 

46 

accomplished  by  Kamlnow,  and  others,  using  the  electro-optic  effect 
which  is  another  form  of  bulk  nonlinear  effect. 

Theoretical  treatments  of  some  of  these  nonlineeur  optical  phenomena 

47  48  49  50  51 

have  been  given  by  Braunsteln,  Bloembergen  '  ^  and  KLelnman.  *  Of 

48 

these,  the  paper  by  Bloembergen  is  perhaps  the  most  thorough,  discussing 
these  higher  order  effects  from  a  point  of  view  that  shows  the  coherent 
nature  of  the  phenomena.  The  approach  used  by  Bloembergen  is,  in  many 
respects,  similar  to  that  used  by  the  author.  The  details  will  be 
considered  in  following  chapters. 

The  original  suggestion  for  the  use  of  these  multiple  quantum  effects 

52 

to  produce  harmonic  generation  is  due  to  E.  T.  Jaynes.  Jaynes  pointed 
out  the  generality  and  importance  of  multiple  quantum  effects  and  B\iggested 
subharmonic  pumping  of  a  two-level  system  as  a  means  of  generating  high 
power  levels  at  high  frequencies.  This  suggestion  of  Jaynes  provided  the 
initial  incentive  for  the  re8eeu:ch  described  in  this  report. 


Although  these  devices  have  sometimes  been  included  under  the 
broad  heading  of  masers,  their  operation  is  essentially  parametric  in 
nature. 
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CBhFUBR  II 


THE  mTESACnOir  OF  RADIATION  UITH  MAmR 

‘  In  this  chapter  ve  shall  study  the  Interaction  of  a  radiation  field 
with  an  atcm  or  molecule.  A  study  of  this  type  will  form  a  basis  for 
analyzing  multiple  quantum  or  nonlinear  phenomena  as  well  as  lineeir; 
single  qviant\mi  phencaaena.  Two  approaches  are  presented.  The  first,  using 
the  concept  of  transition  probabilities,  is  developed  to  include  both 
single  and  do\ible  quantum  absorption.  The  second  is  a  zoacroscoplc  approach 
relating  the  macroscopic  properties  of  a  material  to  its  microscopic 
properties . 

A.  RADIATION  FIELD 

In  order  to  be  rigorously  correct  the  radiation  field  should  be 
quantized.  However,  if  we  are  interested  in  large  field  strengths,  corre¬ 
sponding  to  large  photon  occupation  numbers  for  the  radiation  fields,  we 
may  consider  the  fields  classically.  In  such  an  approach  the  field  is 
described  by  E  and  H  ,  its  electric  and  magnetic  fields,  or  by  its 
vector  and  scalar  potentials  jA  and  0  .  Use  of  this  approach  neglects 

the  effect  of  spontaneous  emission,  which  will  not  be  important  for  the 
microwave  frequency  range  but  may  play  a  role  in  optical  phenomena. 

B.  ATCMIC  SXSTQt 

The  atomic  system  must  be  considered  quantum-mechanically.  The  eqiaatlon 
governing  the  behavior  of  the  quantum  system  is  Schrodinger ’ s  equation 

ifi|^lt>»Wlt>  ,  (2.1) 

where  |t)  is  the  wave  function  of  the  system  and  is  the  Hamiltonian. 

If  there  are  no  applied  radiation  fields  present,  then  the  Hamiltonian 
consists  of  the  sum  of  the  kinetic  and  potential  energies  of  the  atom. 
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We  sball^  for  the  time  being,  neglect  the  Interactions  between  atoms  themselves 
(spin-spin),  and  between  the  atoms  and  the  lattice- (spin-lattice).  This 
Hsmlltonian  is  denoted  and  has  stationary  elgensolutlons  |n)  and 
corresponding  energies  satisfying  the  equation 

JJq  In)  -  E^  |n>  .  (2.2) 

We  shall  assume  that  these  energies  and  their  corresponding  functions  or 
kets  are  known.  The  functions  {|n)^  are  assumed  to  be  a  ccanplete  set. 

If  we  turn  on  the  radiation  field,  then  it  will  interact  with  the 
atom  and  cause  a  change  in  the  Bsmlltonlan  of  the  system.  The  new 
Hamiltonian  may  be  written  in  the  form 

^  (2.3) 

where  is  the  term  due  to  the  interaction  of  the  atom  with  the  radiation 

field.  We  shall  be  Interested  in  the  case  where  the  perturbation  results 
from  a  time  harmonic  field  or  sum  of  such  fields  and  hence  we  may  write  it 

^  0  icopt  p*  -iopt 

Jt'  =  )  (^e  +3r  e  '•^  )  .  (2.4) 

The  quantity  is  real  and  the  superscript  -t  refers  to  the  particular 

✓ 

frequency  present.  For  a  given  applied  field,  W  is  proportional  to  the 
field. 

The  matrix  element  of  an  operator  Q  in  the  energy  representation  is 
defined  as 


=  <n|Q|m> 


(2.5) 
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C.  TRARSinOH  HtQBABILm  MBTBOD  OP  80UJTI0N 
1.  Flret  Order  ProceaseB 

In  solving  the  radletlon  problem  we  asavnae  that  the  total  wave  function 
Is  a  superposition  of  the  unperturbed  eigenstates 

-1  E„t/h 

It)  =  )  ajj(t)  e  |n>  ,  (2.6) 


where  coefficient  of  state  n  .  Substituting  (2.6)  and 

(2.3)  into  (2.1)  gives 


1  fi 


z 


a  W  .  e 
n  nk 


i(\ 


E^)t/fi 


(2.7) 


The  general  approach  to  the  solution  of  this  set  of  equations  is  to  assume 
that  the  frequency  of  the  applied  signal  satisfies  the  condition  htu  wi 
and  that  the  system  is  initially  in  the  state  m  ,  a^(0)  =1  .  If  it 
is  further  assumed  that  we  consider  the  solution  for  times  small  eno\agh 
that  the  approximation  a^  =  1  ,  ^  '  n  /  m  ,  is  velid,  then  we 

have  fran  (2.4)  and  (2.7) 


i  ft 


-  “•)» 


(2.8) 


where  =  (E^^  -  Ejn)/h  and  the  nonresonant  term,  involving  the  sum 

frequency  (flj^  +  w)  ,  has  been  dropped.  Upon  Integrating  (2.8)  subject 
to  the  initial  condition  ^^^(O)  =0  ,  we  have 


a^(t) 


i(0j^  - 

e 


1 


(2.9) 
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The  probability  of  occupancy  of  state  k  is  then  given  by 


ain^  -  <D)t/2 
(“k.  - 


(S.IO) 


If  we  assume  that  the  final  atate  is  not  well  defined  but  rather  is  given 
by  a  density  of  final  states,  then  in  order  to  find  the  probability  of  a 
transition  to  one  of  these  states  we  must  integrate  over  the  density  of 
final  states.  We  shall  assume  that  the  final  state  is  not  sharp  on  account 
of  line  broadening  end  tl’iat  the  density  of  final  states  is  given  by  a  line 
shape  function  g(v)  ,  where 


J  g(v)  dv  =  1  .  (2.11) 

0 

For  a  Lorentzlan  line  we  have 


g(v) 


2  T^ 


1  +  (v  - 


(2.12) 


where  is 

the  system  in 


55 

the  spin-spin  relaxation  time, 
one  of  these  states  i;s  given  by 


The  probability  of  finding 


Prob.  =  J  l\(v)|^  g(v)  dv  .  (2.13) 

0 


If  we  further  assume  that  the  strength  of  the  perturbation  is  small,  then 
we  can  find  a  time  t  large  enough  that 


sln^  (n  -  <o)p/2  rtt 

- - - »  —  6(0  -  oj)  ,  (2.14) 

(tl  -  03)  2 
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consistent  with  the  assuiiptlon  that  at  the  time  t  ,  ^  ^  * 

By  ecnblning  (2*10 .)>  (2.12),  (2.13)/  and  (2.1^)  and  Integrating,  we  have. 


Prob. 


— 


2  T. 


«■  1  -  “f 


(2.15) 


The  transition  probability  per  unit  time  is  given  by 


2  T, 


mk 


1  .  -  0.)^ 


(2.16) 


Since  is  a  Hermitlan  matrix,  we  have  ,  and 


W  ,  =  W. 
mk  km 


(2.17) 


thus  stating  that  the  probability  of  the  system  going  from  k  to  m  is 
the  same  as  its  going  from  m  to  k  .  If  we  assume  that  and 

the  number  of  atomic  systems  in  levels  k  and  m  are  n^^  and  n^  , 
respectively,  then  the  power  absorbed  by  such  a  system  is  given  by  the 
product  (transition  probability) (energy  per  photon) (net  difference  in 
ixjpxilation) 


-  (“,  -  • 


2  T, 


1  ♦ 


(2.18) 


If  the  macrosystem  is  in  e  nomal  population  state,  \  ^  ® 

the  power  absorbed  is  positive.  Naser  action  occurs  when  an  Inverted 
population  is  achieved,  (Oj^  -  n^)  >  0  ,  and  the  power  absorbed  is 

negative,  implying  emission.  This  single  quantvmi  process  is  linear  in 
the  sense  that  the  power  sbsoxbed  or  emitted  is  proportional  to  the  square 
of  the  field  and  hence  to  the  Incident  power. 
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The  asBiauption  leading  to  the  derivation  of  this  relation  are  that 
t  «  Tg  in  order  that  (2.l4)  be  consistent  with  (2.12)  and  that 
t/ft  «  1  in  order  that  la^^j  «  1  .  Ccnibinlng  these  gives  the 
relation 


(2.19) 

If  we  look  upon  as  the  coherence  time  of  the  atom,  (coherence  being 

disturbed  by  "collisions"  occurring  with  a  meeui  time  Tg)  then  (2.19) 
states  that  (WVii)'"  ,  which  is  a  measure  of  the  time  it  takes  to  cause 
a  transition,  must  be  greater  than  the  coherence  time.  In  other  word.s  the 
assumptions  leading  to  (2.l6)  and  (2.l8)  are  satisfied  only  if  the  pertui^a- 
tlon  is  weak  enough  so  as  not  to  introduce  coherence  effects  in  the  atomic 
system. 

2.  Higher  Order  Processes 

The  equations  derived  in  the  preceding  section  assumed  energy  conserva¬ 
tion,  01  3  ,  and  the  absorption  or  emission  of  a  single  quantum  of 

radiation.  These  resvilts  may  be  extended  to  processes  involving  more  than 
a  single  quantum  by  carrying  the  calculation  to  higher  order.  These  higher 
order  processes  may  Involve  quanta  from  different  radiation  fields  or  more 
than  one  quantum  from  a  single  radiation  field.  In  order  to  be  definite 
we  shall  consider  the  latter  and  calculate  a  double -quantum  absorption 
process.  The  Important  aspects  of  multiple  quantum  processes  will  be 
evident  from  this  specific  calculation. 

Let  us  assume,  as  before,  that  the  system  is  initially  in  state  m  , 
aj^(O)  =  1  .  Then  the  equations  for  the  ,  to  first  order,  are  given 

by  (2.9); 


«  1 


h 


SubBtltutlng  these  first  order  solutions  into  (2.7)  gives,  for  to 
second  order, 


axx  .nh, 

»(n_  -  “) 

tm 


(2.20) 


where  we  have  assvned  that  now  2cd  st  and  that  the  frequency  (o  is 

not  near  any  of  the  natural  resonances,  Q.  emd  0  .  Upon  integrating 

Kn  znn  2 

(2. 20)  under  the  above  assumption  and  finding  the  quantity  >  we 

have 


mn  nk 


h(0  -  m) 

'  nm  \ 


sln^(n^  -  2cD)t/2 


(2.21) 


After  comparing  with  the  first  order  case,  (2.10),  we  see  that  the  expressions 
are  similar,  with  (flj^  -  m)  being  replaced  by  (ftj^  -  ape)  and  the  matrix 
element  Ji  .  being  replaced  by  the  sum 


z 


}i 

mn  nh 

-  “) 
nm 


(2.22) 


By  performing  the  integration  over  final  states,  we  have  for  the  transition 
probability 


11 


Z 


ft(Q  -  a>) 
L  n  isn 


2 


1 . 


(2.23) 


and  for  the  power  absorbed 


P  .  ■  (n  -  n,  )  2tto  W  , 

abs  '  m  k'  mk 


(2.2k) 
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The  important  aspects  of  this  double  photon  absorption,  quantitatively 
described  by  (2.23)  and  {2.2^-),  may  be  seen  by  considering  a  single  term 
from  the  sum  appearing  in  these  expressions.  This  term  describes  the  coupling 
of  levels  m  and  k  vlth  the  level  n  acting  as  an  intermediate  or 
"virtual"  state.  The  term  virtual  is  used  because  the  intermediate 
state  (quantum  system  in  level  n  ,  radiation  field  with  N  -  1  photons) 
does  not  conserve  energy  vlth  the  initial  state  (quanttnn  system  in  level 
m  ,  radiation  field  with  N  photons).  The  amount  by  which  energy  is  not 
conserved  is  given  by  the  energy  dencaajnator,  "  ^)  •  total 

probability  of  a  transition  from  state  m  to  k  involves  a  summation  over 
all  possible  virttial  levels,  taking  into  account  the  lack  of  conservation 
of  energy  through  the  energy  denominator. 

The  followltig  conclusions  can  thus  be  reached  regarding  these  two- 
photon  processes: 

(1)  There  must  be  some  intermediate  quantum  state  n  which  is 

connected  both  to  the  Initial  and  to  the  final  state.  This  state  may  be 

a  third  level  or  either  the  initial  or  final  state  itself.  In  the  latter 

case  a  term  of  the  form  or  ,  would  be  involved. 

ami  kk 

(2)  The  strength  of  the  effect  will  depend  on  how  close  the  virtual 
states  come  to  conserving  energy;  the  closer  they  sure  to  conserving  energy, 
the  stronger  the  effect.  The  second  order  effect  will  be  of  the  order  of 
(ii  /h(n  -  Oi))  times  the  first  order  effect. 

(3)  The  ixjwer  absorbed  depends  on  the  fourth  power  of  and, 

since  is  proportional  to  the  applied  radiation  field,  the  power 

absorbed  is  proportional  to  the  fourth  power  of  the  field  or  to  the 
square  of  the  applied  power. 

Second  order  processes  involving  two  different  frequencies  and  higher 
order  processes  involving  more  than  two  qxianta  are  considered  in  a  similar 
manner.  The  third  order  process,  for  example,  will  involve  two  virtual 
states. 

3.  Example 

As  an  example  consider  a  three-level  system  with  unequally  spaced 
levels,  (Fig.  2.1).  Let  the  energy  eigenvalues  be  ,  and 

define  the  natiural  freqxiencies  "  ®i)/*  >  further 

suppose  that  there  are  matrix  elfflnents  of  the  pert\irbation  connecting 
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FIG.  2.1 — Energy  level  diagram  of  a  three-level  system.  First  order 
absorptions  occur  at  and  .  Second  order  absorp¬ 

tion  for  20  *  Is  raown  wlth'^a  "virtual"  level  depicted 
by  the  dotted  lire. 


Relative 

Absorption 


FIG.  2. 3- -Theoretical  absorption  spectrum  at  high  rf  power  levels.  The 
double  quantum  absorption  appears  midway  between  the  two  first 
order  absorptions,  here  shown  saturated. 
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levels  1  and  2,  and  2  and  3#  but  levels  1  and  3  are  not  directly  connected 
in  first  order.  An  observation  of  the  spectrum  of  this  system  would,  to 
first  order,  give  two  lines  centered  at  and  as  shown  in  Fig. 

2.2.  Levels  3  and  1,  although  not  connected  in  first  order,  are  connected 
in  second  order  via  level  2.  The  power  absorbed  in  this  case,  from  (2.24), 
is 


1  +  -  203)^ 


(2.25) 


which  shows  a  peak  at  203  =  =  (0^2  ^21^  *  algebraic  mean  of 

the  two  first-order  lines.  This  is  shown  in  Fig.  2.3*  From  (2.25)  we  see 


that  although  the  double  quantum  abeoiT>tion  line  always  occurs  at  the 


frequency  03  = 


its  strength  will  depend  upon  how  close  to  degeneracy 


work. 


and 

(1) 


Jlg^  are.  Upon  comparing  this  example  with  Kastler's  original 
we  see  that  this  is  precisely  the  effect  he  sawr 


We  may  ai3mniarize  this  approach  by  the  following  obseinrations ;  First 


of  all,  for  the  field  strengths  usually  encountered  in  experimental 
conditions  the  transition  probability  method  correctly  predicts  single 


quantum  absorption  and  emission  processes  as  verified  by  experimentation. 
When  carried  to  higher  order,  the  theoiy  predicts  miiltiple  quantum 
absorption  processes  which  have  been  experimentally  observed.  The 
assumption  on  which  the  derivation  is  based,  Eq.  (2.10),  does,  however. 


limit  the  general  validity  of  the  theoiy  to  field  strengths  satisfying 
this  condition.  Further,  it  is  not  clear  from  this  point  of  view  how 
a  process  such  as  the  absorption  of  two  photons  at  frequency  co  and  the 
coherent  emission  of  a  single  photon  at  the  frequency  2m  would  be  handled. 
Rather  than  pursue  this  approach  further,  we  shall  turn  to  the  macroscopic 
approach,  which  avoids  many  of  the  problems  inherent  in  the  vise  of  tran¬ 
sition  probabilities. and  which  is  used  throu^  the  rest  of  this  paper. 


D.  MACROSCOPIC  APPROACH 

We  may  begin  a  study  of  the  macroscopic  approach  by  writing  down 


Mucvell’s  equations  in  tte  presence  of  matter: 

7  •  2  =  4«p  ,  7  •  B  =  0  ,  (2.27) 

along  with  the  constitutive  irelatlona 

L  =  E  +  WP  ,  B  =  H  +  4jtM  .  (2.28) 

The  characteristics  of  the  matter  are  embodied  In  the  quantities  F  and 
M  ,  the  polarization  and  magnetization^  respectively.  In  particular, 
the  relations 


P  =  £(I)  >  M  =  M(h)  ,  (2.29) 

define  the  response  of  the  material  to  the  applied  fields. If  the 
dependence  Is  known,  then  a  solution  to  Maxwell's  eqimitionB  may  in  principle 
he  found.  It  Is  generedly  assumed  that  2  proportional  to  E  ,  and 
M  to  H  ,  giving  the  tensor  suBceptlbillties,  and  X^  : 

P=:XJ  ,  M  =  XJ  .  (2.30) 

From  these  susceptlhllltles  the  familiar  dielectric  constant  and 
permeability  are  defined  as 

6  =  1  +  4rtX  ,  p  =  1  +  4jtX  ,  (2.31) 

e  m 


where  now 


£»  €B  ,  B  =  pH 


(2.32) 


more  general,  relation  would  be  J  =  £(E  ,  g)  and  M  »  M(E  ,  g)  . 

We  shall  for  simplicity  assume  the  simpler  relations,  (2.29);  where  the 
polarization  depends  only  on  the  electric  field  and  the  magnetization  depends 
only  on  the  magnetic  field. 
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The  above  relations  form  the  basis  for  the  solutions  to  Maxwell's  equations 
for  the  linear  case. 

It  Is  more  generally  true,  at  least  we  may  assume  so,  that  the  polar¬ 
ization  depends  not  only  linearly  on  K  but  on  higher  powers  as  well.  By 
using  a  similar  argument  for  the  magnetization  we  may  write  these  more 
general  dependences  as 


P  =  E  +  X^^)  E  E  +  X^^^  E  E  E  +  .  .  . 

(2.33) 

M  =  X^^^  H  +  X^^^  H  H  +  X^^^  H  H  H  +  .  .  .  , 


where  the  new  tensors  are  of  higher  order.  This  more  general  constitutive 
relation  is  seen  to  give  rise  to  nonlinear  effects.  For  example,  if  E 
varies  as  cos  oot  ,  the  term  E  E  will  have  a  component  at  the  frequency 
2aj  and  through  (2.28)  and  (2.26)  radiation  fields  will  be  set  up  at  the 
frequency  2ta  . 

Eqmtlons  (2.26)  through  (2.33)  constitute  a  description  of  the 
macroscopic  electromagnetic  field  for  both  the  linear  and  nonlinear  cases. 

It  remains,  however,  to  determine  the  response  of  the  medium  to  the  applied 
fields,  Eqs.  (2.29)  or  (2.33).  In  particular  we  should  like  to  determine 
this  response  analytically  from  the  microscopic  properties  of  the  atomic 
systems  making  up  the  medium.  The  basis  for  such  a  calculation  follows 
from  the  relations 


P  =  Njg,M  =  Hm  ,  (2.3li) 

where  P  and  M  as  before  sure  the  macroscopic  polarization  and  magnetization, 

N  is  the  number  of  atoms  per  unit  voliaie  and  p  and  m  are  the  average 

~  (2) 

or  expectation  ved.uee  of  the  polarization  euid  magnetization  per  atom. 

The  average  Implied  by  the  double  beur  constitutes  both  a  statistical 

.V  ^^Wr  example,  the  second  term  in  the  expans  ion.  for  P  would  be 
P'^^^  =  X^jj^  Ej  Ej^  and  the  third  term  would  be  Pj-^^  =  • 

^^^8ee  J.  H.  Van  Vleck,  reference  5^,  for  a  proof  of  this  relation. 
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average  and  a  quantun  macbanical  average,  the  latter  coming  as  a  result 
of  the  inherent  statistical  nature  of  quantum  phenomena.  The  problem 
thus  is  reduced  to  a  quantum  statistical  evaluation  of  £  and  m  in 
response  to  the  radiation  fields  present. 

E.  EVALUATION  OF  THE  MACROSCOPIC  PROPERTIES:  THE  OEISITY  MATRIX  APPROACH^ 

The  density  matrix  formulation  provides  a  direct  means  for  formulating 
quantum  statistical  problems.  The  density  matrix  may  be  characterized  by 
its  equation  of  motion, 

ifip  =  ,  p]  +  —  (p  -  P®)  ,  (2.35) 


and  the  prescription  for  finding  the  expectation  value  of  an  observable, 


(2) 


(Q)  =  Trace  (pQ) 


(2.36) 


In  these  expressions  p  is  the  density  operator,  31  is  the  total  Hamiltonian, 
p®  is  the  value  of  p  at  equilibrium,  and  t  is  the  relaxation  time 
associated  with  the  return  to  eqxiilibrium.  The  value  of  t  will  depend 
on  the  particular  element  of  the  density  matrix  considered.  In  general, 
when  using  the  energy  represisntation,  t  may  be  divided  into  longitudinal 
or  spin-lattice  relaxation  times  sussociated  with  the  diagonal  components 
of  p  and  transverse  or  spln’-spln  relaxation  times  associated  with 
the  off-diagonsd  conqjonents  of  p  .  The  longitudinal  relaxation  times, 
usually  denoted  by  Tj^  ,  describe  the  chswacteristlc  time  in  which 
the  spin  system  and  the  lattice  reach  equllibrliun  with  each  other. 

This  process  involves  the  exchamge  of  energy  between  the  spin  system 
and  the  lattice.  The  transverse  relaxation  times,  usually  denoted  by 
T2  ,  describe  the  characteristic  time  in  which  the  spin  system  reaches 
equilibrium  within  itself.  This  relaxation  process  is  important  in  determining 

^^^For  an  alternates  approach  to  this  calculation  see  Slater, Chapter 
6,pp.  154-157,  Kramers,^  Chapter  8,  pp.  48o,  and  Bloembergen.  Briefly, 
this  method  involves  a  calculation  of  the  perturbed  wave  function  due  to  the 
applied  fields  and  the  evaluation  of  the  polarization  and  magnetization  from 
these  wave  functions.  This  method  assvmies,  however,  that  the  frequency  of 
the  radiation  field  is  not  near  any  natural  transition  frequency  and  hence  the 
field  causes  no  transitions.  The  method  used  here  is  not  limited  by  this 
assumption. 

'  'These  relations  are  derived  in  Appendix  A. 
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In  Eq.  (2.36)  the  trace  ie 

(2.37) 

,  (2'.  38) 

where  |n)  is  any  complete  set  of  functions.  It  is  usually  found  convenient 
to  define  this  set  of  functions  as  eigenfunctions  of  the  unperturbed 
Hamiltonian,  ,  where 

ii^ln)  =  E^|n>  .  .  (2.2) 

In  this  representation  the  time  average  value  of  the  diagonal  terms,  , 

when  normalized  such  that  £  =  1  ,  may  be  interpreted  as  the  probability 

of  occupancy  of  level  n  .  The  off-diagonal  elements,  p^  ,  are 
Interpreted  as  giving  a  measure  of  the  amount  of  coupling  between  levels 
m  and  n  ,  where  the  coupling  is  proportional  to  the  magnitude  of  the 
element.  The  procedure  for  finding  the  values  of  the  observables  is  first 
to  solve  the  density  matrix  Eqs.  (2.35)  under  the  particular  conditions 
of  the  problem.  From  this  we  have  a  knowledge  of  all  the  elements  p 

nm 

The  magnitxide  and  phase  of  an  observable  quantity  is  dependent  on  the 
elements  of  p  through  (2.36)  euid  is  generally  dependent  on  both  the 
diagonal  and  off-diagonal  components  of  p  : 

"I  %  Sm 

m,n 

(2.39) 

'T'' 

Pnn  }  Prim  * 

nn  nn  nm  mn 

n  n 

mj^n 

The  matrices  of  the  observables  are  assumed  to  be  constants.  Independent 
of  time,  and  hence  the  time  dependence  of  the  observable  is  found  from  the 


the  effects  of  coherence  in  the  atomic  system, 
taken  of  the  product  of  the  matrices  defined  as 


p^  =  <n|p|m> 


nm 


Q  =  <nlQ|m> 
nm 
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time  depmdence  of  p  .  examining  the  equations  of  motion  for  p  , 
we  can  see  vhat  the  time  dependence  is,  and  in  particular  we  can  see  how 
linear  and  nonlinear  processes  arise. 

Let  us  write  Eq«  (2.33)  in  matrix  form  in  the  following  way: 


n  a  m 


ifip  +  ^  (p  -  P*  ) 
^nn  T  '^nn 


'■•fn  -  I®-'") 


n  m 


ifip 


ms  +  (E  -  E  )p  +— SE=)(  (p  -p  ) 

nm  D  n'  nm  _  nm  '  mm  nn 


(2.41) 


Prom  these  equations  we  can  see  first  of  ed.1  that  the  "natural  frequency" 
of  the  diagonal  components  is  zero  while  for  the  off-diagonal  term  p^ 
it  is  (Ej^  -  Ejj)/*  •  Clearly  a  given  ccarponent  will  be  most  strongly 
driven  when  the  frequency  of  the  driving  term  is  near  its  natural  frequency. 

1.  Linear  Effects 

II  —  I  1  ■■  IIW  ■  ^ 

Direct  coupling  of  levels  n  and  m  occurs  via  the  term  (p  -  P  ) 
®  /  nm'  mm  '^nn' 

and  is  strongest  when  varies  with  time  as  exp  i(E  -  E  )t/ft  .  As 

nm  m  n  / 

a  result  of  the  direct  coupling  term,  p^  is  proportional  to  and 

hence  to  the  applied  field.  Evaluation  of  an  observable  through  Eq. 

(2.39)  also  gives  a  linear  dependence  on  the  applied  field.  Hence  the 
direct  coupling  terms  are  seen  to  give  rise  to  linear  effects,  among 
which  are  those  described  in  Eqs.  (2.32). 

2.  Honlinear  Effects 

In  addition  to  these  direct  coupling  terms  there  are  additional  terms 
of  the  form  p_^  which  describe  the  indirect  coupling  of  levels  n 
and  m  via  level  <6  .  These  terms  lncl\]de  the  right  hand  side  (rhs) 
of  (2.4o)  and  the  last  two  terms  on  the  rhs  of  (2.4l).  A  term  of  this 
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form  describes  the  fact  that  if  levels  -t  and  m  are  coupled,  /  0  > 
and  if  -t  is  edao  coupled  to  n  ,  0  »  then  n  and  m  will  also 

be  coiffled.  If  is  proportional  to  ,  then  is  proportional 

,  and  hence  an  observable  dependent  .on  will  be  quadratic  in  the 

applied  fields.  A  tern  of  this  form  thus  gives  rise  to  a  nonlineeu*  de¬ 
pendence  on  the  field  and  accounts  for,  tha  second  term  in  Bii.  (2.33).  If 
p_£^  is  nonzero  because  of  indirect  coupling  to  s(»e  level  q  , 
then  p^  is  proportional  to  ^4q  is  thus  a  third  order 

quantity  in  the  radiation  field.  An  observable  dependent  on  p^  would 
thus  be  proportional  to  the  cube  of  the  field  and  would  give  rise  to  third 
order  nonlinearities. 

We  thus  see  that,  in  this  particular  formalism,  nonlinear  effects 
are  handled  by  means  of  a  consideration  of  indirect  coupling  mechanisms. 
Processes  of  any  order  may  be  handled  by  this  means. 

The  usual  method  of  solution  of  such  a  problan  may  be  summarized  as 
follows: 

(1)  For  the  particular  problem  at  hand  allow  for  all  the  radiation 
fields  which  are  present  by  Including  appropriate  terms  in  the  Hamiltonian. 

(2)  Assume  solutions  for  the  components  of  the  density  matrix  of  the 
form  of  harmonic  series  of  the  frequencies  present. 

(3)  Find  a  steady  state  solution  by  equating  terms  with  the  same 
time  dependence.  This  results  in  a  set  of  algebraic  equations  which  are 
then  soluble  by  standard  methods.  In  practice,  simplifications  can 
usually  be  made  which  will  make  these  solutions  easy  to  obtain. 

Once  the  density  matrix  is  found,  then  through  (2.3^)  and  (2.36)  the 
macroscopic  properties  are  known  and  Maxwell's  equations  (2.26)  through 
(2.29)  may  be  solved  under  the  particular  conditions  of  the  problem. 

These  equations  then  form  a  self-consistent  means  of  solving  the  general 
radiation  problem. 

3.  Examples 

Let  us  briefly  consider  two  examples  using  the  density  matrix  approach, 
a.  Spin  1/2  System 

Assume  that  a  spin  I/2  system  is  in  a  dc  magnetic  field  and  that  the 
natural  frequency  of  the  transition  is  given  by  (E-  -  E.,  )/li  =  H  .  Choose 
the  representation  to  be  the  one  in  which  S  and  S  are  diagonal. 
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«bere  the  aetrlcee  of  the  x  and  z  coeiponente  of  the  magnetic  dipole 
operator  are 


--C  3  ’  “‘‘“(o-O 


(2.42) 


If  an  rf  magnetic  field  cos  oot  is  applied  in  the  x-direction,  then 
the  interaction  Brailtonlan  is 


^  ®1 


,  (2.43) 


and  hence  the  nonzero  matrix  elements  of  are 


^  * 

^12  “^21  =  - 


(2.44) 


Substituting  (2.44)  into  (2.4o)  and  (2.4l),  using  the  longitudinal  and 
transverse  relaxation  times,  and  Tg  ,  defined  in  Appendix  A, 

and  employing  harmonic  balance  gives  the  following  solutions  for  p  : 


•i  Tg  (Pu  '  ^22^  ^ 

2»  (1  -  iTgS) 


(2.45) 


^21  ^  ^12 


(2.46) 


^11  “  ^11 


(pjl  -  P^)  ^  ^  Tg  [sin  2ci3t  +  TgS  cos  2a3t] 

4fi^a)(l  +  T?  6^) 


(2.47) 


P22  -  1  -  P^j_ 


(2.48) 
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and 


P 


o 

11 


1  + 


T^Tg  \?  ^ 

ft^d  +  6^) 


(2.49) 


where 


6  =  ft  -  o)  .  (2.50) 

The  qviantitles  and  p|g  are  the  thermal  equilibrium  values  of  p 

while  p^j^  and  Pgg  are  the  dynamic  equilibrium  values  or  the  average 
values  with  the  rf  flelda  present.  Equation  (2.4-9)  thus  describe®!  the 
effects  of  saturation.  Prcm  (2.39)>  (2.42),  (2.45)  and  (2.46)  we  find 
the  expectation  value  of  the  x  component  of  the  magnetic  dipole  moment 


<'‘x>  •  “  <<>12  ♦  “ai’ 


2  .f  Iff!  /  ^  \ 

^  ®1  ^2^^11  ”  *^22^ 


„  ^  (sin  ojt  +  T„  6  cos  cut)  .  (2.51) 

fid  +  B^)  2 


The  power  absorbed  is  given  by 


abs 


/-  ■  K  \ 


•  ('•u  -  4)  *“ 


time  average 

\  2fi/  1  +  T?(ft  - 


T|(ft  -  a))2 


(2.52) 


By  multiplying  by  the  number  of  particles  N  emd  noting  that  N(p°j^  -  p^g) 
=  n^  -  Og  ,  we  see  that  (2*52)  is  the  same  as  (2.l8),  calculated  by 
transition  probabilities,  where 
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If  ve  now  evaluate  {u  )  ,  the  expectation  value  of  the  z-coaiponent^ 

of  the  fiipole  moment,  we  find  from  (2.39),  (S.te),  (2.47)  and  (2.48) 


^  "  ‘*<^11  ■  ^22^ 


(2.53) 


2(Pii  -  ®1  '^2 


=  u(p..  -  p„«) - *5 - ^  p  (sin  2t0b  +  T„6  cos  2ajt) 

fi^a)(l  +  B^)  2 


which  shows  a  dc  component  as  well  as  a  ccaaponent  at  the  frequency  2cu  . 
This  component  at  2cd  conies  as  a  result  of  the  indirect  coupling  of  level 
1  to  itself  In  second  order  via  level  2,  and  similarly  for  level  2  via  1. 
When  this  latter  ccmiponent  of  u  is  Introduced  Into  Maxwell's  equations 
it  will  give  rise  to  radiation  fields  at  the  frequency  2ai  ,  an  effect 
not  predicted  by  the  method  used  to  derive  Eq.  (2.l8).  Hence  we  see  that 
the  density  matrix  formvilatlon  gives  a  more  general  approach,  giving  more 
Information  than  do  transition  probabilities, 
b.  Three-level  3ystem 

For  the  second  example  let  us  consider  the  three-level  system 
previously  discussed  in  section  C.2,  this  time  using  the  density  matrix 
approach.  To  be  definite,  let  us  assume  magnetic  dipole  transitions  and 
an  applied  field  cos  ojt  .  The  nonzero  perturbation  terms  are 

/ 

^12  =  -  Pi2  ^ 


^21  “  "  ^21  ^ 

(2.54) 

/ 

^23  “  ■  ^23  ® 

3*32  “  "  ^32  ® 
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By  aseumiLng  that  -  X^^  -  X^^  , 

then  to  first  order  we  have  the  aolutlons^^^ 


.  «  ^*^11  ~  ^22^  ^12  ^ 
^  26  -  oi)  +  i/Tg] 

X  -  '  ^33^  >^23  h 

26  [(fl^g  -  m)  +  i/Tg] 


(a. 55) 


(2.56) 


and 


X^2  *  0  >  first  order 


(2.57) 


A  calculation  of  the  magnetic  sment  using  these  values  would  predict 
first  order  absorption  at  oi  =  flg^^  and  an  =  n^g  .  By  keeping  second 
order  terms  we  find  that  levels  1  and  3  are  coupled  since  level  2  is 
coupled  to  both  1  and  3  in  first  order.  Simplifying  the  exact  expression 
somewhat  by  assianlng  that  (fl^  -  ta)  »  l/Tg  and  2to  =  gives 


i  ^2  ^12  ^23  ^ 
l6^  6 


(2.58) 


where  8  =  flg^  -  (JJ  .  Finally,  using  this  value  for  X^^  ,  we  find  that 
X^g  and  Xg^  take  on  nonzero  third  order  terms.  In  the  case  of  X^ 
this  occurs  as  a  result  of  the  fact  that  levels  1  and  3  are  connected 
(in  second  order)  and  3  is  in  turn  connected  to  2,  A  similar  argument 
applies  for  Xg^  .  These  third  order  terms  are 


^12  ^32  ^23 


(2.59) 


— rry 

detailed 


See  Chapter  V  for  a  derivation  of  these  equation^  and  for  a  more 
explanation  of  this  problem. 
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and 


(3)  ^  ^2  ^12  ^21  ^23  ^  /  0 

‘23  “  ^^^33 


(2.60) 


where  the  time  dependence  Is  still  e^^  .  Calculation  of  the  magnetic 

moment  from  the  relation 


<H>  =  (^2.2  ^21  *  ^21  ^12^  ^^23  ^32  *  ^32  ^23^  '  (2*6l) 

and  subsequently  evaluating  the  power  absorbed,  gives  Eq.  (2.25)  where  we 
have  set  2ai  =  *  Th'is  not  too  unexpectedly  we  arrive  at  the  same 

answer  as  was  given  by  the  transition  probability  method.  One  important 
point  can  be  made,  however.  If  the  selection  rules  for  the  transitions 
were  such  that  was  not  zero,  then  from  (2.39)  and  (2.58)  we  immedi¬ 

ately  find  a  nonzero  component  of  magnetization  at  2a>  and,  as  mentioned 
before,  this  will  generate  radiation  fields  at  the  new  frequency.  A  more 
detailed  analysis  of  this  particular  case  is  given  in  Chapter  V. 


F.  GEKERAL  ASPECTS  OF  NONLINEAR  QUANTUM  EFFECTS 

From  the  formulation  of  the  problem  presented  above  we  can  draw 
some  general  conclusions  about  the  various  nonlinear  effects  which  are 
possible. 

a.  The  strength  of  the  nonlinear  effect  measured  by  the  magnitude  of  the 
nonlinear  polarization  or  magnetization  will  be  proportional  to  the 
magnitude  of  the  higher  order  components  of  the  density  matrix.  These, 
in  turn,  will  Involve  combinations  of  first  order  terms  in  the  density 
matrix.  Hence  the  higher  order  effects  will,  in  general,  be  large  when 
the  corresponding  linear  effects  are  large.  The  latter  are  found  to  be 
largest  when  the  matrix  elements  connecting  the  various  levels  are  large 
and  when  energy  is  nearly  conserved,  i.e.,  when  the  applied  frequency  is 
near  a  natural  transition  frequency  (if  the  applied  signal  is  resonant 
with  a  natural  frequency  then  the  strength  of  the  transition  is  Inversely 
proportional  to  the  llnewldth  of  the  transition).  We  then  see  that 
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nonlinear  effects  will  be  strongest  when  the  frequencies  of  the  racllation 
fields  are  near  the  natural  frequencies  of  the  transitions  involved  and 
when  these  transitions  are  strong.  The  dependence  on  the  proximity  to 
natural  resonant  frequencies  points  out  the  fact  that  for  a  given  material 
the  nonllneeu*  effects  will  be  frequency-sensitive.  This  is  in  contrast 
to  a  diode,  for  examplej  where  the  nonlinearity  is  relatively  frequency  in¬ 
sensitive. 

b.  In  general,  the  higher  the  order  of  the  effect  (the  larger  the  number 
of  photons  involved),  the  smaller  is  the  magnitude  of  the  effect.  In 
practice  second  and  third  order  nonlinear  effects  will  probably  find  the 
most  application. 

c.  The  total  energy  of  the  field  and  atomic  system  must  be  conserved  in 
the  overall  Interaction.  This  may  be  achieved  either  by  the  radiation 
field  and  the  molecule  separately  conserving  energy,  or  by  an  increase 
in  energy  of  one  system  being  compensated  by  a  corresponding  decrease  in 
the  energy  of  the  other.  In  the  latter  case  the  change  in  the  energy  of 
the  quantiun  system  must  correspond  to  the  difference  between  two  of  the 
eigenenergies  E  ,  defined  in  (2.2). 

d.  The  type  of  nonlinearity,  i.e.,  whether  even  or  odd,  and  the  number 
of  photons  Involved,  depends  upon  the  number  of  levels  involved  and  on 
the  selection  rules  for  the  pertinent  transitions.  The  selection  rules 
are,  in  turn,  determined  by  the  symmetry  properties  of  the  atomic  system 
and  its  surroundings.  For  example,  if  a  crystal  has  a  center  of  inversion, 
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then  only  odd  order  effects  due  to  electric  dipole  transitions  are  allowed. 

e.  In  any  nonlinear  effect  there  will  in  general  be  contributions  to  the 
effect  from  many  combinations  of  levels.  It  may  be  the  case,  however, 
that  the  contribution  of  two  or  three  of  the  levels  dominates,  in  which 
case  the  essential  properties  may  be  found  from  a  consideration  of  the 
simpler  two-  or  thfee-level  system. 

G.  EXAMPLES  OF  NONLINEAfi  EFTTICTS 

We  may  consider  briefly  some  of  the  more  important  effects  possible 
in  two-  and  three-level  systems.  Sane  of  these  effects  are  considered  in 
detail  later  in  this  report,  some  are  considered  by  other  authors,  and 
some  have  not  yet  been  investigated. 
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1.  Two-level  Sbrteai 
a.>  Hamanlc  Oeoeretitm 


Prom  an  examination  of  the  density  matrix  Bqts.  (2.4o)  and  (2.4l)  for 
the  specific  case  of  a  two-level  system,  it  can  easily  be  seen  that  the 
off-diagonal  congKJnente  contain  only  odd  powers  of  the  applied  field,  while 
the  diagonal  components  contain  even  powers.  From  (2.39)  we  see  then, 
that  if  a  dipole  moment,  which  may  act  as  a  source  of  radiation,  has 
dlagonsJ.  ccxnponents,  then  generation  of  even  harmonics  is  possible, 
whereas  if  it  contains  only  off-diagonal  components,  thent  only  odd 
harmonic  generation  is  possible.  The  question  of  diagonal  and  off-diagonal 
components  of  the  dipole  is  determined  by  the  type  of  interaction  (magnetic 
or  electric  dipole)  and  upon  the  symmetry  of  the  Hamiltonian.  Second 
harmonic  generation  is  briefly  discussed  in  part  E.3  of  this  chapter  and 
third  harmonic  generation  is  considered  in  detail  in  Chapter  III. 
b.  Parametric  Processes 

There  are  several  types  of  parametric  processes  possible  in  a  two- 
level  system.  They  may  be  separated  into  two  classes  depending  on  whether 
the  dipole  operator  possesses  a  diagonal  component  (^^^  or  ji^g)  ,  or 
not. 

If  the  system  possesses  either  a  or  a  then  the  two  processes 

shown  in  Figs.  2.4  and  2.^  are  possible. 


FIG.  2. 4- -Parametric  process  withi  three  frequencies  present; 
or  pgg  nonzero. 
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FIG.  2.5“”Fara3netrlc  process  'with  two  frequencies  present; 

Ull  or  ^22 


In  the  first,  shown  in  Fig.  2.4,  a  strong  field  is  applied  at  the 

frequency  oo^  ,  and  two  radiation  fields  at  frequencies  cUg  and 

are  generated  by  the  nonlinear  action.  This  is  a  threshold  effect  and 

requires  =  a)_  +  cu-  in  order  to  conserve  energy.  Such  a  system  is 
i.  d  5  05 

essentially  considered  by  Suhl.  The  second  effect,  shown  in  Fig.  2.5, 
involves  only  two  radiation  fields,  at  and  cOg  ,  and  is  usually 

referred  to  as  a  Raman  effect.  From  section  F  above,  we  have  =  ojg 
+ (Eg  -  •  Neglecting  the  effects  of  spontaneous  emission,  this  is 

also  a  threshold  effect.  Javan  has  considered  this  as  a  form  of  Raman 
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maser. 

For  systems  with  a  dipole  mcsnent  containing  off-diagonal  components 
two  parametric  processes  are  also  possible  and  are  shown  in  Figs.  2.6 
and  2.7* 
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FIG.  2.6 — Parametric  process  with  four  frequencies  present. 


PIG.  2.T--Parametric  process  with  three  frequencies  present. 


In  the  first,  ^Ig.  2.6,  four  radiation  fields  are  present,  their 
frequencies  satisfying  (Oj^  +  .  Here,  for  example,  fields 

applied  at  and  will  generate  fields  at  and  cdj^  when  the 
magnitude  of  the  applied  fields  exceeds  a  threshold  value  determined  by 
the  parameters  of  the  problem.  In  the  second  case  the  quantum  system 
resonance  replaces  one  of  the  radiation  fields  and  the  frequencies  must 
satisfy  03^  +  ojg  =  CB^  +  (Eg  -  E^)/fi  .  This  is  also  a  threshold  effect 
with  a  field  at  generated  when  the  fields  applied  at  and  cOg 
exceed  a  given  threshold.  In  the  last  two  cases  considered  it  is  possible 
for  oXj^  to  equal  Og  ,  reducing  the  actual  number  of  fields  present. 

The  second  of  these  effects  is  considered  in  Chapter  IV  while  the  first 
has  not  been  considered  in  the  literature. 

In  the  case  of  the  two-level  system  we  find  that  the  presence  of  a 
diagonal  component  of  the  dipole  operator  allows  processes  involving  one 
less  radiation  field  than  for  systems  with  only  off-diagonal  matrix  elements. 
This  may  be  seen  by  comparing  the  processes  described  by  Pigs.  2.4  and  2.6 
and  those  shown  in  Figs.  2.5  and  2.7* 

There  are  several  possible  applications  of  these  parametric  effects. 

As  in  the  case  of  classical  p«arametric  systems,  these  quantim  systems  may 
act  as  amplifiers  when  operated  below  threshold  and  as  oscillators  when 
operated  above  threshold.  The  two  cases  shown  in  Figs.  2.5  and  2.7  have 
a  slightly  different  character  than  classical  peurametric  effects.  Here 
the  quantum  system  acts  as  the  "idler,"  removing  the  necessity  for  providing 
an  electroi-fi^netlc  field  at  this  frequency. 

2.  Three-level  System 

The  nvnnber  of  nonlinear  effects  possible  in  a  three-level  system 
exceeds  those  possible  in  a  two-level  system  as  a  result  of  the  additional 
number  of  combinations  of  levels  possible.  We  sliall  thus  pci'nt  out  only 
some  of  the  important  possible  applications.  The  first  is  a  three- 
frequency  mixing  process,  shown  in  Pig.  2.8. 


FIG.  2.8— Three-frequency  mixing  process; 
**12  ^  **23  ^  **13  nonzero. 

-  - 


H»e,  the  application  of  any  two  frequencies,  ,  o:^  or  » 
where  +  oig  =  oj^  >  will  generate  fields  at  the  third  frequency  if 
all  the  transitions  are  allowed.  A  special  example  of  this  is  the  case 
in  second  harmonic  generation  where  ’  considered  in  Chapter  V. 

Such  a  mixing  process  could  also  act  as  an  up-  or  down-conveirter. 

The  reverse  of  this  process,  namely  where  fields  at  are  applied 
and  parametric  oscillations  at  oa-  and  ox,  occur  above  a  given  threshold, 

d  X.  -o 

has  been  considered  as  a  form  of  amplifier  by  Anderson. 

Javan,  and  Shlmoda  and  Yajima^  have  considered  a  process  similar 
to  that  of  Fig.  2.8,  in  which  only  two  radiation'  fields  are  present.  This 
is  shown  in  Pig.  2.9  where  now  oa^^  =  oa^  +  (Eg  -  E^)/fi  and  only  the  1-3 
and  2-3  transitions  need  be  allowed.  This  effect  is  very  similar  to  that 
considered  in  Pig.  2.5  where  now  the  coupling  to  the  third  level  replaces 
the  self-coupling  or  p^g)  .  This  is  also  a  type  of  Raman  effect. 


PIG.  2.9“-Parametric  oscillator; 


*^13  ^23 


nonzero . 


Finally,  let  us  consider  the  three-frequency  process  shown  in  Fig. 

2.10.  E 

“3 


- 

FIG.  2.10 — Ihree-frequency  parametric  process;  p^g  and  Pg^  nonzero. 
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Here  It  le  aeeumed  that  there  exlet  nonzero  dipole  nments  connecting 
levels  1  end  2  and  2  and  3*  !Rie  application  of  radiation  fields  at 
and  (1^  will  cause  a  field  to  be  generated  at  the  frequency  ,  when 
the  level  of  the  applied  fields  exceeds  a  given  threshold.  If  oi^  »  , 

then  this  process  Is  seen  to  provide  a  fom  of  parametric  amplification 
where  the  pump  frequency  Is  near  that  of  the  signal.  This  has  not  been 
previously  considered. 

H.  SUMHAHY 

In  this  chapter  nonlinear  as  well  as  linear  effects  have  been  discussed 
and  their  analysis  has  been  presented  from  the  transition  probability  and 
density  matrix  points  of  view.  The  latter  approach  Is  the  more  general 
and  Is  best  suited  for  the  types  of  problems  where  multiple  quantum 
processes  are  Involved.  The  nonllnearltles  encountered  make  possible 
harmonic  generation  and  various  forms  of  parametric  processes  which  may 
be  used  for  amplification  and  mixing. 

Since  the  natural  frequencies  of  atoms  extend  from  radio  frequencies 
to  optical  and  above,  these  forms  of  nonlinear,  multiple  quantum  effects 
will  similarly  occur  throughout  the  spectrum.  The  transitions  involved 
may  include  paramagnetic,  rotational,  vibrational,  and  electronic  energy 
levels . 
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CHAPTER  III 


THEORY  OF  THIRD  HARMONIC  GENERATION 
IN  A  TWO-LEVEL  QUANTUM  MECHANICAL  SYSTEM 

In  this  chapter  we  shall  consider  the  theory  of  third  harmonic 
generation  In  a  two-level  quantum  mechanlcaLL  system.  Such  a  system  may 
consist  of  only  two  levels  as  In  the  case  of  a  spin  l/2  magnetic  dipole^ 
or  a  particular  pair  of  levels  In  the  CEuse  of  a  more  ccmipllcated  atomic 
structure.  The  classical  interpretation  of  the  problem  Is  presented 
briefly,  followed  by  a  detailed  quantum  mechanical  analysis  In  which 
various  aspects  of  the  Interaction  are  discussed.  The  results  of  the 
analysis  will  be  applied  to  the  specific  case  of  third  harmonic  generation 
In  a  gas  where  the  Inversion  transition  In  Is  used  as  an  example. 

A  conq)arison  of  these  predictions  with  the  experiment  perfonned  on 
NH^  are  presented  in  Chapter  VI. 

A.  GENERAL  CONSIDERATIONS 

The  problem  to  be  examined  is  the  interaction  of  an  electromagnetic 
field  of  frequency  cn  ■  2jrv  with  a  quantum  mechanical  system  consisting 
of  a  pair  of  levels  with  energies  E^  and  Eg  and  a  natural  frequency 
G  =  (E^  -  E^)/h  .  Such  a  process  may  take  place  via  either  an  electric 
dipole  or  magnetic  dipole  Interaction  (electric  quadrupole  and  higher 
order  Interactions  will  not  be  considered)  and  the  order  of  the  interaction, 
i.e.,  the  number  of  photons  involved,  will  be  determined  by  the  type  of 
Interaction  Involved.  These  may  be  summarized  as  follows: 

1.  Electric  Dipole  Transitions 
a.  Induced  Electric  Dipole 

When  the  Interaction  la  pvure  electric  dipole  In  character,  provided 
neither  the  upper  nor  lower  state  possesses  a  permanent  dipole  mennent, 
then  only  odd  order  Interactions  may  occur,  I.e.,  only  an  odd  number  of 
photons  may  be  absorbed  or  emitted.  This  type  of  Interaction  will  occur 
when  eeush  state  has  either  even  or  odd  parity. 
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13.  Permanent  Electric  Dipole 

If  either  of  the  etatee  posseasee  a  permanent  dipole  moment  or  a  dc 
electric  field  Is  applied  (creating  a  permanent  electric  moment  via  a 
mixing  of  states  then  even  order  Interactions  are  allowed  ^  well 
as  the  odd  order  ones.  Permanent  moments  may  exist  when  a  state  has  neither 
even  nor  odd  paurity. 

2.  Magnetic  Dipole 

If  the  transitions  between  states  are  due  to  pure  magnetic  dipole 
Interactions j  then  three  possible  situations  may  be  considered: 

a.  H  _  Perpendicular  to  H, 

rr  ac 

If  the  rf  magnetic  field  Is  orthogonal  to  the  dc  field  (assuming  a 
free  spin  and  no  crystalline  fields)  then  the  only  transitions  which  az% 
allowed  are  of  the  form  Am  >  i  1  ,  where  m  Is  the  magnetic  quantum 
nimiber.  In  order  to  go  from  one  level  to  the  other,  a  net  change  In  m 
of  1  Is  required.  This  Is  seen  to  occur  only  when  an  odd  nmber  of 
photons  is  involved. 

b.  Parallel  to 

When  the  rf  magnetic  field  is  parallel  to  the  dc  field,  the  photons 
carry  no  angular  momentum.  In  this  case,  no  transitions  sure  allowed  between 
states  which  have  different  m  values. 

c.  Arbitrary  Relative  to 

With  the  rf  field  in  an  arbitrary  direction,  it  may  be  divided  into 
components  parallel  to  and  perpendicular  to  ;  these  photons  carry 

0  and  ±  1  units  of  angular  momentum,  respectively.  Transitions  between 
states  can  then  occur  with  any  number  of  photons. 

The  phenomenon  of  second  harmonic  generation  involves  the  absorption 
of  two  photons  and  the  emission  of  a  single  photon  at  twice  the  frequency. 
From  the  discussion  above,  such  a  process  Is  seen  to  occur  In  either  a 
magnetic  dipole  system  of  type  (2c)  or  in  an  electric  dipole  system 
possessing  a  x>ennanent  dipole  mcxnent  (lb).  Third  hamonlc  generation  may 
occur  In  magnetic  cases  (2a)  end  (2c)  and  for  either  electric  dipole  case. 

For  the  two-level  s^rstem,  the  electric  dipole  transitions  have  the 
distinct  advantage  over  their  magnetic  counterparts  in  that  the  value  of 

^^'In  the  absence  of  state  mixing  and  for  states  of  either  even  or 
odd  parity,  we  have  the  quantity  (ijxjl)  =  0  ,  which  means  zero  permanent 
electric  dipole  moment.  If  a  dc  field  Is  applied  or  the  state  has  neither 
even  or  odd  parity,  then  we  have  (l|x|l)  /O  ,  and  a  permanent  dipole  exists. 
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the  electric  dipole  moaeat,  (in  C*0>8*  unite  where  E  »  E)  ,  is  typically 
100  tines  that  of  the  magnetic  dipole  moment.  Since  the  output  power 
generated  varies  as  some  high  power  of  the  dipole  moment,  (eighth  for  third 
harmonic,  sixth  for  second  harmonic)  magnetic  dipole  interactions  are  of 
little  interest  compared  to  their  electric  counterpart.  In  this  chapter 
we  shall  consider  third  harmonic  generation  in  an  electric  dipole  system. 
Magnetic  dipole  transitions  will  be  considered  later  in  Chapter  V,  pertain¬ 
ing  to  three-level  systems  in  which  other  factors  are  made  to  compensate 
for  the  smallness  of  the  dipole  moment. 

B.  PHXBICAL  PIGIUHE  OF  THE  EROBliM 

At  this  point,  before  launching  into  a  detailed  quantum  mechanical 
study  of  the  problem,  it  is  perhaps  appropriate  to  consider  it  from  a 
physical  point  of  view.  Let  us  assume  that  we  have  an  atom  in  which  there 
is  a  single  electron  bound  to  the  nucleus  by  some  form  of  potential  well. 
For  simplicity,  assvoae  the  electron  to  be  a  loceilized  particle  with 
coordinate  x  (rather  than  a  smeared  out  wave  function)  and  assume  that 
the  potential  is  expanded  about  its  equilibrium  position.  From  texts  on 
classical  mechanics^^^  the  potential  (assumed  one-dimensional)  may  be 
written 


(3.1) 


where  the  first  term  is  an  additive  constant  to  the  zero  of  energy  which 
may  be  set  equal  to  zero  without  loss  of  generality  and  (dV/dx)xQ  =  0 

^^^For  example,  see  Ref.  57,  P*  319* 


-  39  - 


since  the  point  Is  asstmed  to  be  a  point  of  eq[ulllbrlum.  Taking 
Xq  =  0  ,  we  have 

V(x)  -  I  x^  +  i  Vg  x3  +  ^  V3  x^  +  .  .  .  ,  (3.2) 

where  the  constants  *  ^2  '  ^3  derived  from  the  coefficients 

in  the  Taylor  expansion  of  V(x)  .  We  shall  be  Interested  in  the  qualita¬ 
tive  details  of  V  and  not  in  its  quantitative  nature  since  in  practice 
such  detailed  knowledge  is  not  usually  available.  In  the  usual. case,  for 
small  excursions,  the  first  term  predominates  in  the  expansion  and  one 
obtains  simple  periodic  motion  with  a  natural  frequency  given  by  the 
relation 


01^ 


(3.3) 


If  the  particle  has  a  charge  e  ,  then  the  application  of  an  oscillating 
electric  field  E  cos  cA  will  cause  the  charge  to  oscillate  at  the  frequency 
0)  in  response  to  this  drive.  The  equation  of  motion  for  the  particle, 
assuming  only  the  first  term  in  the  expansion  for  V  to  be  kept,  is 

d^x  2  eE 

— 5  >»•  0^  X  =  —  cos  oA  .  (3*^) 

dt  m 


The  solution  is  easily  seen  to  be 


X 


e£ 

—  cos  oA 
m 

2  2 
%  ■  “ 


(3.5) 


where  the  particle  oscillates  at  the  frequency  of  the  applied  electric 
field  and  the  strength  of  the  oscillation  is  proportional  to  the  applied 
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field  ai>d  depends  on  the  relatlm  of  the  exiled  to  netrual  frequencies 
through  the  resonance  denoninator.^^^ 

The  polarization  or  dipole  moment  induced  hy  ^he  notion  of  this  charge 
is 


2 

e  E  cos  (ut 
p  =  ex  =  ---g - ^ 

m(oQ  -  0)  ) 


giving  the  susceptibility-  per  atcm 


e 


2 


X  = 


(3.6) 


(3.7) 


If  now  the  magnitude  of  the  fundamental  excitation,  measured  by 
^  t  Eq.  (3.5)^  le  increased,  the  higher  order  terms  in  the  expansion  of 
V(x)  become  increasingly  Important  and  must  be  Included.  In  order  to 
be  specific,  let  us  assume  that  the  potentiail  function  V(x)  is  symmetric, 
Vg  =  0  ,  and  the  next  nonvanishing  term  in  (3.2)  is  x^/4  ,  By 

including  this  in  the  equation  of  motion  for  the  charge,  Eq.  (3.^)>  have 

2  3 

d  X  V-  x^  eE 

- +  CD^  X  +  — * -  =  COS  cut  .  (3.8) 

dt*^  °  mm 

To  first  order  in  the  perturbation  E  the  solution  is  simply  that 
given  by  Eq.  (3.5).  Upon  substituting  the  first  order  solution  into 

3 

(3.8),  expending  the  term  containing  x  ,  and  using  the  trigonometric 
identity  cos^cut  *  i/h  (3  cos  oat  +  cos  3tut)  ,  we  find  that  in  order  to 
satisfy  the  equation  to  the  order  of  we  must  include  in  x  a  term 

^^^Por  completeness  there  should  be  a  term  in  the  differential 
equatlto  (3.4)  which  is  proportional  to  x  and  which  describes  the 
damping,  both  dissipative  and  radiative.  This  term  will  remove  the 
infinite  response,  for  x  in  Eq.  (3.5)  when  cu  ■  cUq  . 
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at  •  Btononlc  balance  yields 


X 


3 


1  ^  COB  3tct 

m  (9  <0^  -  a>^)(a>^  -  a>^)^ 


(3.9) 


This  component  of  x  gives  rise  to  a  term  In  the  polarization  at  the 
frequency  3o)  given  by 


P(3aj) 


e  X, 


^/e\^  cos  3cut 

h\mj  (90)^  - 


(3.10) 


which  when  incorporated  into  Maxwell's  equations  gives  rise  to  third 
harmonic  fields.  We  thus  have  a  simple  classical  picture  of  the  non¬ 
linearity. 

From  a  quantum  mechanical  point  of  \riew,  the  discrete  charge  located 
at  the  point  x  is  replaced  by  the  charge  density  e  ♦*(x)  i|f(x)  ,  where 
t(x)  is  the  wave  function  in  the  gchrodinger  picture.  In  this  picture 
the  motion  of  the  center  of  gravity  cf  the  charge  density  is  equivalent 
to  the  motion  of  the  particle.  Wv*  may  imagine  the  motion  of  this  charge 
cloud  oscillating  in  response  to  the  applied  fields  as  setting  up  a  dipole 
moment  of  the  form 


P 


I 


t*e  X  dx 


(3.11) 


The  characteristics  of  this  motion  will  determine  the  linear  and  nonlinear 
character  of  p  .  The  exact  motion  of  the  charge  cloud,  and  hence  p  , 
is  determined  by  the  time  evolution  of  the  wave  function  which  In  t\im  is 
governed  by  Bchr6dinger' s  equation. 

C.  qUAinaH  theory  of  third  HARMOlilC  aEREBAnOR  IR  A  TND-LEVEL  SYSTEM 

The  quantum  system  to  be  considered  will  be  taken  to  consist  of  two 
levels  with  energies  and  E^  In  order  of  increasing  energy  where 


the  natural  freqxwncy  Q  is  defined  by 

0  -  (Zg  - 


(3.12) 


as  Illustrated: 


FIG.  3*1— Energy  level  diagram  for  two-level  system. 

The  Hamiltonian  of  the  system  is  assumed  to  consist  of  the  unperturbed 
term  and  an  interaction  term  ,  such  that 

,  (3.13) 

where  the  matrix  elements  of  the  two  are  given  by 

(iPHqIi)  = 

<2iMQl2>  =  Eg  (3.14) 

<lpiQ|2>  =  <2;Hq|i>  =  0 
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and 


<1  |3<'|  2>  -  -  >1  B 

(2  lift'l.  1);  «  -  ^  B 

<lpi  |l>  =  (2\H  |2>  =  0 


where  E  is  the  total  electric  field  present.  The  above  matrix  elements 
eo:e  taken  In  the  energy  representation  In  which  Is  diagonal  and  the 
phases  of  the  unperturbed  wave  functions  cure  taken  so  as  to  make  the  dlx>ole 
matrix  elements  real.  It  Is  assumed  that  the  type  of  interaction  is  induced 
electric  dipole;  hence  we  have  the  zero  diagonal  components  of  51  .In 
order  to  make  this  apply  to  the  magnetic  case  one  would  simply  change  E 
to  H  and  restrict  the  rf  magnetic  field  to  lie  in  a  plane  normal  to  the 
dc  field. 

Equations  of  Motion 

The  equations  of  motion  for  the  two-level  system  in  the  density  matrix 
notation  are 


(p-.  -  p®  )  1  / 


(3.16a) 


.  (P22  ~  P22^  ^  Aj'  It'  \ 

^22'*“  “  .  ^21  ^12  ■  ^12  ^21^ 
Tn  111 


(3.16b) 


^12  ^  "  ^^12  ^^^22  ■  ^11^  ^12 


(3.16c) 


^21  *  ^  ^^21  ■  ^^11  ■  ^22^  ^21 


,  (3.l6d) 


where  T^  Is  the  longitudinal  relaxation  time,  Tg  is  the  transverse 
relaxation  time,  euad  p^  and  p^  are  the  values  of  the  diagonal  elements 


-  4k  - 


of  the  density  netrlx  in  the  absence  of  applied  radiation  fields.'  *  In 
general^  and  p^  will  he  given  by  the  Boltzmann  factor  when  the 

system  Is  In  theival  ocmtact  with  a  reservoir.  If  the  system  is  at  a 
teaqperature  T  ,  thm  ve  have 


4  -  Pjl  “  ,  (3.17) 

where  k  Is  Boltzmann's  ccmstant. 

Let  us  assume  that  there  are  present  two  electric  fields:  the  first, 
cos  ojt  ,  is  the  applied  fundamented  field  at  the  frequency  oj  and  is 
assumed  to  be  strong;  the  second,  E^  cos  ('3ajt  -  0)  ,  is  the  third  harmonic 

field  which,  for  generality,  may  be  arbitrary  in  origin  but  in  the  case 
ultimately  considered,  will  be  generated  from  the  nonlinear  properties  of 
the  molecular  system.  The  total  field  E  is  then 

E  =  E^  cos  oat  +  E^  cos  (3£ut  -  0)  .  (3.l8) 


The  parameters  p  may  be  introduced  in  the  following  manner: 


(3.19a) 

(3.19b) 

(3.19c) 

,  (3.19d) 


brief  explanation  of  the  relaxation  times 
in  Appendix  A. 


and  T-  is  given 


where  now  the  strength  Of  the  interaction,  |jE  ,  is  measured  in  the  units 
of  angular  fi'equency.  TUe  conjugate  is  left  on  (even  though  the  quantity 
is  real)  in  order  to  retain  the  symmetric  form  of  the  resulting  eqviatlons. 

The  interaction  Bamlltonian  is  then 


31,. 


’12 


21 


Pi* 


loib 


*  -lajt 

Pi® 


P3® 


130* 


) 


(3.20) 


Substituting  (3*20)  into  (3*l6)  gives 


Pll* 


^11  ■  ^11 


=  i(p2i  -  ®i2)vPie  +  *  fif  ^  ) 


(3.21a) 


P22 


P22  “  P22 


4/.  .  \/o  ^  -  i3cot  .  -i3a3t\ 

i(Pl2  -  Pai'^Pl  Pi  +  ^36  ) 

(3.21b) 


^12  *  -f  -  -  ‘(^22  - 

"  (3.210) 


»21  *  -f  *  “‘>21  -  ‘(‘>11  -  ‘>22>(V““  b‘"~'  * 


(3.21d) 


We  now  look  for  a  steady  state  solution  to  the  al;  ve  set  of  equations. 
The  usual  procedure  in  solving  these  equations  for  the  interaction  of 
radiation  with  matter  is  to  assume  the  diagonal  components  to  be  constants 
and  the  off-diagonal  components  to  vary  as  e*^  .  Such  sun  assumption 

leads  to  the  usual  linear  or  single  quantm  interactions.  In  order  to 
include  the  effect  of  higher  order  interactions  involving  more  than  a 
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single  quauitum,  a  more  general  solution  Is  required.  Thus,  ve  assume  the 
following: 


12CDt 


+  e-*'**  (3.22a) 


-  xl .  X&  .  x(:"> .- 


laoDt 


22  22  22 


22 


.  X^)  .  X^")  (3.22b) 


=  X^>  .  X^l)  a-“  .  x^)  .  x^3>  (j.g,,) 


-  x(^)  .  x(-^)  e-^  .  x(3)  e‘3“  .  xi;3)  e-^3^,(3.22d) 


21  21 


21 


21 


where  the  diagonal  components  contain  only  even  harmonics  and  the  off- 
diagonal  components  contain  only  odd  harmonics.  The  superscripts  denote 
the  corresponding  frequency  dependence  and  the  X  are  complex  constants. 
Terms  up  to  third  order  are  kept  in  P21  since  these  will  he 

seen  to  give  rise  to  the  third  harmonic  generation  (fifth  order  terms  would 
give  rise  to  fifth  harmonics,  etc.,  but  only  third  harmonic  generation  is 
here  considered).  Although  terms  to  second  order  in  the  diagonal  components 
wo\ild  he  sufficient  to  show  third  harmonic  generation,  there  is  a  frequency 
detuning  effect,  to  he  discussed  later,  whose  correct  magnitude  requires 
the  inclusion  of  the  fourth  order  terms  in  and  .  It  can  he 

shown  that  even  order  terms  in  the  off-diagonal  components  and  odd  order 
terms  in  the  diagonal  comppnents  are  related  and  vanish  identically. 

^^^All  the  coefficients  of  these  terms  could  he  linearly  related  to 
the  average  or  dc  component  of  the  off-diagonal  components  X|2  • 

From  the  random  phase  argxnnent  of  statistical  mechanics  (see  Appendix  A 
for  references)  the  off-diagonal  components  of  p  are  identically  zero 
in  equilibrium;  hence  X^  =  X|^  =  0  . 
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From  the  haelc  reXatlons  governing  the  density  matrix  for  a  two- 
level  system, 

^11  **■  *^22  *  ^  » 

and  (3.23) 

^12  ""  ^21  ' 

the  following  relations hips  between  the  constants  X  may  be  derivedi 


4i* 

4-1 

.  (3.24a) 

(±2) 

11 

(3.24b) 

w . 

11  * 

.4*)  =  o 

(3.24c) 

Ail) 

^2 

II 

(3.24d) 

.  (3.24e) 

After  substituting  the  assumed  solutions  (3*22)  into  the  equations 
(3*21),  using  (3*2Ub)  and  {3.2hc)  and  employing  harmonic  balance,  the 
following  algebraic  equations  are  found  for  the  diagonal  components: 


X°  -  X« 
All  A, 


i  T, 


'••21  ^12  ^  Pi  *  '>^21  *12  '  Pi 


(3.25a) 

+  (4l^  -  ^ll^)  Pc 


■'21 


12 


(a»  -  i)  x(f)  .  (4J;)  -  P3 


(3.25b) 


*  (4l^  "  ^12^^  ^1 
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and 


(20,  ♦  i)  -  (4^)  -  X^^))2* 


(3.25c) 


^^21  ^12 

(W  -  ^)  -  xj^b  ^3  ■»•  -  X^OPi  (3.25d) 


(■“♦I:) 


-  (4^)  -  4^>)pt  -  (43) .  43)),; 


(3.25e) 


and  the  off-diagonal  c(3tQponent8: 

[(n  -  4d)  +  ~]  x[^)  =  A^3^  +  2x[^^  3*  +  2xJJ^  3*  +  2x5i^^  ^3  (3.25f) 

[(ft  +  0))  -  ^]  =  A +  2x[’^^  e,  2X.52)  p*  +  p*  {3>2?g) 

[(n  -  3a))  +  ^]  =  A  +  2^^'  ‘3i  'i  2X^'  (^.?5h) 

[(a  •»•  30))  -  ^]  X^J^  =  A  +  2X.^^^  ,  (3.251) 

where 

AaxJ^-X^  .  (3.26) 


^e  quantity  A  represents  the  average  difference  between  the  diagonal 
components  when  the  system  Is  In  equilibrium  with  the  radiation  fields. 
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The  quantity  may  he  defined  aa  the  value  vhen  the  syatem  Is  In 
thermal  equilibrium  with  its  surroundings; 

•  <3.27) 

When  multiplifcdi  by  the  tpL;  .iiaber  of  particles,  N  ,  the  quantity  IfA 
r'^preeents  the  average  prpi  ;•  ..on  difference  n^  -  where  n^  and 
are  the  nviiter'  cf  ttcle^  l!.<  lower  and  upper  states,  respectively. 

Ey  uair'7  !H.25b;  ,^h  (3.25e),  the  diagonal  components  may 

be  eliminated  from  (3.25f)  ::'OU‘’h  (3*251)  to  yield  a  set  of  equations 
involving  only  the  orf-diagcnax  t'“  niis.  In  order  to  simplify  the  algebra 
somewhat,  it  shall  be  assimed  that 


»  m-  >  ST" 
^2  ^1 


(3.28) 


which  physically  is  equivalent  to  saying  that  the  linewidth  of  the  transition 
is  much  less  than  the  frequency  of  operation.  For  a  gas,  this  implies  that 
the  operating  pressure  is  not  too  high.  Upon  using  (3*28)  in  the  diagonal 
equations  and  assuming  ,  the  off-diagonal  equations  become: 


[“ 


-  O)  + 


+  0)  - 


12 


^21 


-  X, 


(3) 

12 


^<3) 


) 

(3.29a) 

) 


(3.29b) 


30)  +  —  Pj^P*  + 

2u) 


30)  -  -  P,P,  - 

20) 


p  +  —  (x<5;^  -  x<^b  (3.29c) 

A  P3  +  —  (X<p  -  X<^b  .(3.29d) 
^  0) 
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It  viU  be  sbovn  later  that  the  hamonlc  polarization  is  directly 
proportional  to  and  .  Since  we  are  interested  in  those 

conditions  giving  rise  to  a  large  third  harmonic  polarization  we  may 
examine  Eqs.  (3*2|9)  for  conditions  consistent  with  this  goed.  !nie- 
maximization  of  the  X's  will  be  consistent  with  the  minimization  of  the 
determinant  of  the  coefficients  in  (3.29).  Prom  an  examination  of  this 
determinant  we  see  that  there  are  two  frequency  regions  where  such  a 
maximization  may  occur,  namely  o  «  0  and  a>  m  n/3  . 

These  two  situations  correspond  to  the  fundamental  frequency  near 
the  natural  transition  and  to  the  case  where  the  third  harmonic  is  near 
the  natural  frequency.  We  shall  consider  these  cases  separately,  thereby 
simplifying  somewhat  the  mathematics. 


D.  SOLUTIOHS  TO  ISE  QUANTUM  MECHANICAL  EQUATIONS 
1.  Solution  for  Case  I;  to  «»  q/3 

Let  us  assume  that  the  fundamental  or  pump  frequency  is  approxim??heiy 
one-third  of  the  natural  transition  frequency,  a>  «  n/3  ,  ^nd  that  th<= 

fundamental  fields,  measured  by  ,  are  sufficiently  stroqg  to  cause 

appreciable  harmonic  generation,  yet  not  so  strong  that  the  assumption 
<  cc  is  invalid.  The  upper  limit  on  the  fields  imposed  by  the  latter 
restriction  is  approximately 


E  (volts/cm)  <  3,500  v(kMc) 


(3.32) 


-l8 

for  a  system  with  a  dipole  moment  equal  to  one  debye  (lO"  ”  esu)  .  Si-rxctly 
speaking,  the  approximations  involved  in  the  solution  will  be  valid  to  the 
order  (dj^/u)  which  will  be  less  than  0.1  if  we  satisfy  the  coridition 
(dj^/oj)  <  1  .  This  assumption  is  made  in  order  to  make  the  mathematics 


tractable . 

By  combining  (3.29a)  and  (3.29b)  and  solving  for 


and 


> 


> 


ve  have 


in  teras  of 


end  X. 


(3) 

12 


(3.33) 


and 


It  la  to  be  noted  that  althoiigh  X^^  and  X^^^  are  nonvanishing  to 
first  order  in  p^  ,  expr.'iissions  (3.33)  and  (3*3^)  include  terms  up  to 
fifth  order  in  p^  (the  next  nonvanishing  term)-,  and  also  retain  the  term 
l/Tg  compared  to  2  ,  even  thoiigh  by  assumption,  2  »  l/Tg  .  It  is 

necessary  to  Tetain  these  smaller  quantities  in  order  to  evalxiate  expressions 
(3.25a)  (which  will  be  found  to  describe  the  phenomenon  of  saturation)  as 
many  of  the  first  order  terms  cancel.  For  the  evaluation  of  other  quantities 
the  higher  order  terms  may  be  dropped.  Simplifying  (3*33)  and  (3*3^)  ly 
dropping  terms  of  the  order  of  l/Tg  and  keeping  only  first  order  terms 
in  p^  gives 


(3.35) 


a  +  CD 


(3.36) 
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and  the  derived  ^imntlty 


-  2(0  A 

'  2 ' 

a  -  (0 


(3.37) 


Tha  third  hamonlc  components  of  and  may  be  found  by 

substituting  (3.37)  into  (3.29c)  and  (3.29d).  These  are 


,(3) 

^12 


P3- 


2 


ft  -  3(0 


jr _ 

2(0 


(0 


(3.38) 


2  K 


X^3)  =  A. 
^21  ^ 


^3- 


a2  2 

ft  -CD 


ft  3(0  - 


2(0 


(3.39) 


We  are  Interested  in  the  case  where  co  >»  ft/3  •  Examination  of  (3.38)  and 
(3.39)  shows  that  X^^  will  be  larger  than  Xg^^ 
order  (ft  +  3(n)/(ft  -  3m)  »  which  for  these  conditions  will  be  very 


by  a  factor  of  the 


large.  We  thus  drop  X, 


(3) 


.(3) 


'21 


■'12 


in  comparison  to 
(3) 

The  condition  for  maximizatlou  of  X^'  will  correspond  approximately 
to  the  vanishing  of  the  real  part  of  the  resonant  denominator  of  (3.29)  or 


ft  -  3(0  + 


3PiPj 


=  0 


(3.1^0) 


2(0 


By  solving  for  the  frequency  (o  we  obtain 


ft  9 

<0  »  -  (H-  —s 

3  2ft‘^ 


(3.41) 
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or  after  defining  an  effective  natural  frequency  n'  by 


we  have  the  condition  for  optimum  interaction 

=  n73 


(3.42) 


(3.43) 


The  new  frequency  Cl'  la  the  effective  natural  frequency  under  the 
conditions  of  a  strong  applied  field  at  a  frequency  near  ClJ'i  .  In 
terms  of  the  quantity  Cl'  the  components  of  the  density  matrix  at  3u) 
are 


2  p- 


2 

n  -  cj 


and 


«  0 


(3.4M 


(3.44b) 


From  expressions  (3.33)^  (3*34),  (3*44)  and  the  defining  relations  (3.24d) 
and  (3.24e),  the  off-diagonal  components  of  the  density  matrix  are  known. 
Next  we  must  evaluate  expression  (3.25a)  whiclvas  will  be  seen, describes 
the  phenomenon  of  satxiration.  Substituting  (3*33),  (3*34)  and  (3*44)  into 
(3.25a)  gives 

-  4®,  /  4P^P*  _ 

-ii - —  =  i  A  { - 5 - ^  4-  —g  •A  '  ' ■  ■'  ‘  -w  " 

i  \t2(0^  -  <t»^)  (0  -  oj^)*!!  +  0^  (O'  -  3a>)  ] 

2  P3P3  ^2 _ 4T2  4  s*^p^) 

1  +  T|(n'  -  yof  (n^  -  u)®>[i  +  (n' 
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If  we  consider  the  lladting  oMe  where  the  harmonic  fields  generated  enre 
small,  (O®  -  this  expr^aaion  becomes 


*\3 


.  1 4  J-M-  *  .  ■  -- 

i  Tj  |t2(0®  -  (D®)  (a®  -  a?)*tl  +  ^  (O'  -  3l>)‘l 


.  (3.W) 


Now  by  using  the  fact  that  =  1  and  +  X^g  =  1  ,  we  may 

show  that 


X?i.X^=i(A-A«) 


(3A7) 


whsre 


A* 


(3A8) 


6 

’(.'Ith  Xj^^  Xgg  being  the  thermal  equilibrium  values.  Substituting 

(3.^7)  into  (3.^5)  and  solving  for  A  yields 


.1 

A  -•  i?.  <1  ■}. 

I 


,*\3 


T  + 


8  Tg  (P^pp 


Tg(n^  -  (Q^  -  Oi^f  [1  +  0^  (n'  -  3cof] 


,,  ^  Vl _ iViiLl^lV\(3.49) 

1  -t-T^  (a'  -  yaf  {a^  -  o^)  [i  +  (n'  -  3a>)^Jj 


Jpctn  using  the  same  assumptions  leading  to  (3.46)  we  have 

4.4' Lar 

T  HtX  -  »")  (»"  -  »¥  U  .  3|(a' - 


30))^] 


1-1 


'(3.50) 
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with  the  knowledge  of  A  given  by  (3***9)  or  (3*50)  and  the  off- 
diagonal  elements  (3»33)»  (3*3^)  and  we  have  sufficient  knowledge 

of  the  quantum  system  to  ceLLcxaate  its  response  to  the  fundamental  fields 
CLiid  the  txteiiiy  to  which  it  v  .11  generate  the  desired  third  harmonic  fields, 
for  the  case  a>  «  fl/3  • 

2.  Solution  for  Case  II;  o  » 

The  second  case  of  Interest  is  one  in  which  the  frequency  of  the 
fundamental  radiation  is  near  the  natiural  transition  frequency,  O)  w  fl  . 
Here  the  third  harmonic  fields  will  occur  at  approximately  three  times  the 
natural  frequency  in  contrast  to  the  case  previously  considered  where  the 
output  was  near  the  natural  resonant  frequency.  To  find  the  response  of 
the  system  under  these  conditions  we  shall  solve  Eqs.  (3.29a)  through 
(3.29d)  xinder  the  new  assumption  that  tn  •«  along  with  those  previously- 
used,  1/^2  * 

First  of  all  consider  Eqs.  (3-29a)  and  (3.29b)  for  the  fundamental 

cQmtxDnents  of  and  X-,  .  Since  we  have  chosen  ti)  w  £1  it  is  cleeur 

(1^  (1) 

that  X^'  will  be  much  greater  than  and  hence  we  may  neglect  the 

latter  quantity.  Considering  the  solution  to  (3«29a)  to  first  order  only 
(neglecting  the  much  smaller  fifth  order  term)  gives 


T 

n  -  a>  +  — 


0)  Tg 


(3.51) 


By  setting  the  real  part  of  the  denominator  of  (3.51)  equal  to  ze-'-o  and 
defining  the  solution  to  the  resulting  equation  £l"  one  obtains 


where 


PlPf 

Q"  =  ft(l  -±5^) 


(3.52) 


(3.53) 
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The  quantity  n"  Is  the  effective  natural  frequency  of  the  system  under 
the  influence  of  the  radiation  fields  applied  near  the  natural  frequency 
of  the  system. 

l^n  using  expression  (3»52)  for  ,  and  *  0  in  (3*29c) 
and  (3.29d)  we  obtain 


(3.5»^) 


14-  a 


(3.55) 


where  in  all  nonresonant  terms  we  have  set  a>  =  and  dropped  terns  of 

2  2 
the  order  of  p  compared  to  those  of  the  order  of  fi  .  Substituting 

(3.52),  (3-5*^)  and  (3.55)  into  (3.25a)  gives 


^11  ■  ^11  _  2  ^  3  A  Pgpf 


i  T^  1  +  t|(0"  -  a>f  k  n^(n"  -  0)  -  ^) 


(3.56) 


i*q3 


3  AP*Pj^ 


k  a^ia"  -  CO  +  i 

^2 


The  first  term  on  the  right  hand  side  will  predominate,  the  last  two 
being  of  the  order  of  (P^/o)^  smaller  than  the  first.  Using  (3.*<-7) 
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•Qdi  first  turn  t^ym 


if 

^'rrzsss:  ■ 

"  *  x  >  ,^(a  - .  »)* 

where  A*  is  defined  hy  (3.1t0),  With  expressions  (3.^),  (3.51*^),  (3.35) 
end  (3*?7)  ve  have  the  density  matrix  for  the  case  of  puntping  near  the 
natural  frequency. 

Kith  the  knowledge  of  the  density  matrix  in  '^he  two  cases  of  interest 
we  are  now  in  a  position  to  determine  the  reaction  of  the  quantum  syst«n 
on  its  surroundings,  namely  a  resonant  cavity. 

E.  MACROSCOPIC  FOLAIUZASIOI  ABD  CAVPR  REACIIOH 

He  turn  now  to  the  question  of  the  dynamics  of  the  electromagnetic 
(E>M)  field  and  in  particular  to  the  reaction  of  the  quantvmi  system  on 
it.  From  classical  electricity  and  magnetism  we  know  that  the  electro* 
magnetic  fields  may  he  derived  from  •  ,  and  ^  ,  the  scalar  and  vector 
potentials,  or  if  the  charge  density  is  zero  a  gauge  may  be  chosen  where 
9  B  0  and  both  jg  and  H  may  be  determined  from  the  relations 

i  =  -^ltA  •  (3-58) 


The  vector  potential  is  determined  from  the  relation 


(3.59) 


where  £  is  the  total  currant  density  ard  Is  composed  of  several  types 
of  currents  as  follows; 

Ref.  p.  117. 
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(v)  tnw  vl^  ti«ii«pori«tlon  of  tnw  chargM.' 

(b)  Poljgrin^ioo  curvMtat  eurrents  that  arise  trm  the  change 
polaris&tlwt  vith  tlae. 

(c)  Napurtiaatlw  eunnrafeas  stationary  eurr«ats  that  flow  within 
regltms  whl^  ara  laaeeesalhle  to  ohaenratlixi  but  whleh  night 
give  rise  to  net  boundary  or  volume  currents,  due  to  laqwrfeet 
orbit  ca'iCeUatlni  on  an  atoailc.  scale. 

(d)  Ckmvectlve  currents:  currents  due  to  the  notion  of  a  nedlum 
bB  a  whole. 

1.  Macrcr.copic  Polarlaatloo 

Our  Interest  will  be  in  polarisation  currents,  (b),  since  we  are 
considering  electric  dipole  effects  Inbound  states  of  atons  and. molecules, 
not  free  electrons,  and  bulk  motions  clearly  will  not  be  of  Interest.  Ihe 
polarization  current  which  su:ts  as  a  soxirce  for  the  E<lf  fields  is  given  by 

>  (3*60) 

where  P  is  the  polarization  or  the  macroscopic  dipole  moment  per  unit 
voliane. 

The  macroscopic  polarization,  P  ,  will  be  equal  to  the  expectation 
value  of  the  polarization  for  a  single  atom  tines  the  number  of  systems 
per  unit  volume,  N  : 


P=.ll£  =  N<ju>  .  (3.61) 

The  expectation  V8l.vie  of  a  single  dipole  system  (j|i)  may  be  calculated  from 
the  density  matrix  p  by  the  reiaticxi^^^ 

e  Tr  (p^x)  ,  (3.62) 

where  p  and  jj^  are  the  matrices  of  the  density  matrix  and  dipole  mostent 
^^^See  Appendix  A  for  the  derivation  of  this  relation. 
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operators,  respectively.  The  vector  notation  is  to  he  interpreted  as 


<*Sc>  ■ 

) 

with  similar  relations  for  p  and 

y 

the  matrices  are 

/^n  ^12^ 

V21  ^22/ 

and  the  trace  of  their  product  is 

Tr(pp)  = 

By  the  previous  choice  of  the  phases  of  the  wave  functions  we  have 
=  =  p  j  hence  (3.65)  reduces  to 

p  =  p(p^  pg^)  ,  (3.66) 

where  now  we  assume  that  the  direction  of  the  qxisntlty  p  is  that  of  the 
fundamental  field  E  .  Returning  to  the  definitions  of  p^g  and  Pg^ 
in  terms  of  the  X's  and  siibstituting  into  (3*61)  gives 

*  4l^^  e'^^]  .  (3.67) 


(P  P^,) 


,  (3.63) 


i 


p  .  In  the  representation  uspd, 
z 


ro  p 


12 


f  p  “ 


:2i 


,  (3.6M 


*^21  ^21  *^12 


(3.65) 
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By  using  the  follovlng  definitions 


Pi*) 

(3.68a) 

p(-) 

(3.68b) 

pw 

=  Np(x[|^  +  X^^b 

(3.68c) 

p(-) 

3 

=  Np(x[‘3^  +  X^'^b 

(3.68d) 

the  polarization  becomes 

P  =  p(-^)  e^  +  p(-) 

e-^  +  P^+)  e^^CDt  ^  p(-)  ^ 

(3.69) 

displaying  components  at  both  the  fundamental  and  the  third  harmonic. 

The  magnitudes  of  these  components  axe  related  to  the  quantum  mechanical 
properties  by  Eqs.  (3*68). 

2.  Cavity  Reaction 

To  be  specific,  in  the  problem  at  han(^  it  will  be  assumed  that  the 
quantum  system  Is  placed  inside  a  cavity  which  is  resonant  at  both  the 
fundamental  and  third  harmonic  frequencies.  The  resonance  at  the  funda¬ 
mental  Is  provided  in  order  to  achieve  the  required  field  strengths  at 
’’easonable  applied  power  levels.  Provision  of  a  resonance  at  the  harmonic 
matches  the  impedance  of  the  source  to  that  of  the  outside  world,  increasing 
the  power  generated  by  the  nonlinear  action  of  the  sample. 

In  order  to  calculate  the  power  generated  by  the  component  of  the 
polarization  at  we  must  first  determine  the  fields  it  sets  up  in 

the  cavity.  This  may  be  formally  expressed  by  the  relation 


(3.70) 


I 

\ 

i 

\ 


where  X  «  X '  -  IX^  and  the  quantities  B  and  P  vary  as  e  . 

The  quantity  X^  is  called  the  cavity  susceptibility  and  relates  the 
field  in  the  cavity  to  its  source,  in  this  case  a  time-varying  polarization. 

Let  us  assume  that  the  polarization  in  the  cavity  is  a  vector  function 
of  the  coordinates  of  the  cavity,  P(r)  ,  and  let  the  orthonormal  eigen¬ 
functions  of  the  cavity  be  u  (z*)  •  To  simplify  things  somewhat  let  us 

further  assume  that  there  is  Just  one  resonance  near  the  frequency  cn  , 
denoted  cd^  ,  and  the  loaded  Q  of  the  cavity  is  given  by  .  By 

letting  the  field  E(r)  be  given  by 

S(£)=  =  )io^£^  ^  (3-71) 


it  can  be  shown  that 


2  ^  „,l/2 


kma  f^PV 


E  = 


2  2^  ^0 
^0 


(3.72) 


where 


/  £(£>  •  SoU)  ‘ fo •  <3.73) 

cavity 

The  quantity  f^  is  a  filling  factor  and  V  is  the  volume  of  the  cavity.  '*■ 
In  this  case  from  (3.70)  and  (3*72)  we  have 


1 

X 

c 


hrtu?  f^ 


(3.71+) 
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It  should  be  pointed  out  that  In  all  previous  expressions  involving 
the  fields  and  ,  these  are  the  fields  at  a  particular  point  in 

space  and  are  not  to  be  confused  with  the  mrmallzed  quant'' *y 

3.  Output  Power 

The  power  coupled  out  of  a  cavity  where  the  fields  are  generated  from 
within  is  given  by  the  relation 


out 


(3.75) 


where  cUq  is  the  angular  frequency  of  the  output  radiation,  W  is  the 
energy  stored  In  the  cavity,  and  is  the  external  Q  (which  for 

optimum  coupling  equals  the  Q  of  the  source).  The  energy  stored  in 
the  cavity  is 


W  =  ^  /  (E^  +  H^)  dv  ,  (3.76) 

cavity 

where  the  bar  signifies  a  time  average  and  we  assume  e  =  1  ,  p  =  1  , 

the  sams  as  for  free  space.  At  resonance  we  have  E‘  -  =  h2  ,  so  (3*76) 
becomes 


W 


f 

cavity 


(3.77) 


If  we  assume  that  the  third  harmonic  field  is  written  as  the  product  of 
a  function  of  time  (giving  the  amplitude)  and  a  fuiiotion  of  space  (the 
normal  mode  function) 


-  If.  -  «!)]  (,) 


(3.78) 
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vhere  jli  Is  an  arbitrary  phase  factor^  then  ve  have 


''3-S’‘3*3 


(3.T9) 


Combining  (3*75)  and  (3*79)  gives 


Pout^3“)  = 


(3a>)E-  E„ 


(3.80) 


By  using  the  results  derived  so  far  we  shall  now  calculate  the  harmonic 
power  out  in  the  two  cases  considered. 


F.  OUTPUT  POWER  AMD  REACTION  AT  THE  FUIQAMENTAL 

1.  Case  I;  ob  w  fl/3 
a.  Output  Power 

We  shall  now  evaluate  the  haniionlc  power  for  the  case  where  the  applied 
fields  are  at  a  frequency  approximately  one-third  the  natural  transition 
frequency.  In  this  case  combining  (3.68a)  with  (3.33)  and  (3.3^)  and  (3.68c) 
with  (3.^^)  and  (3.^5)  gives,  for  the  dipole  moments, 


2Nn  n  A  ^  (X^J^  -  X^^) 


r>2  2 

n  -  CD 


(3.81) 


P^"^^  =  NjiA 


^  -  . I  ' 

'■3  n2  2 

n  -  CD 


n'  -  3a>  + 


•  (3.82) 
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The  component  of  the  polarization  at  3oJ  ,  >  Is  seen  to  consist 

of  two  terms;  the  first  is  proportional  to  (or  to  E-)  and  the  second 
is  proportional  to  (or  Ep  .  They  have  somewhat  different  character¬ 
istics. 

The  first  term,  proportional  to  ,  is  the  polarization  induced  in 

response  to  the  harmonic  field  E^  •  It  is  the  linear,  first  order 

response  of  the  system  to  fields  at  the  frequency  3sJ  such  as  would  occur 

were  :  a  signal,  at  this  frequency  applied  to  a  system  with  a  natural  frequency 
n'  .As  such,  this  component  of  will  have  the  same  spatial  dis¬ 

tribution  as  the  normal  mode  fields  in  the  cavity. 

The  second  component  of  P^^^  ,  proportional  to  ,  is  the 

polarization  generated  by  the  nonlinear  action  of  the  media;  it  is  this 

canponent  of  the  polarization  which  acts  as  the  source  of  the  third  harmonip 

fields  in  the  cavity.  The  spatial  distribution  of  this  component  will  be 

equal  to  the  cube  of  the  spatial  distribution  of  the  fundamental  fields 

(P-  a  E:;  ,  etc.).  Since  this  spatial  distribution  is  not  the  same  as 

jX  lx  /^\ 

that  of  the  other  component  of  one  must  allow  for  a  different  filling 

factor.  Let  us  denote  the  filling  factor  for  the  first  component  by  f^ 

and  that  of  the  second  by  f^  .  Then  combining  (3*72)  and  (3*82)  gives 


-i0 

e 


i  N^A 


i2q  -ir  i-i 

1  + - ^  (30)  -  0)^)  fi'  -  30)  +  — 

-  o)^  J  L  T2  J 


(3.83) 


where  is  the  cavity  resonant  frequency  near  3o)  and  is  its 

loaded  Q  .  Upon  moving  the  first  term  on  the  right  to  the  left,  using 
the  definition  of  p^  (3* 19c)  and  reeurranging  the  denominator  slightly, 
we  have 


<1  + 


U*  NA  Tg  f  ^ 


a  1  +  (30)  -  0)3)  [1  -  iT  (£1'  -  3o))] 

o)_  J  L  , 


8«n  NA  Tgf^  pj 


i2Q, 


1  + 


0)- 


(30)  -  0)^) 


a  * 

)  1  - 


iT2(a'  -  3o)) 


[o2  2‘ 

n  -  0) 
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)^ing  a  quality  factor  tor  tlw 

aa 


.  (3.85) 

i.ar  abaorptlon  coefficient, 
available.  In  Appendix  D, 


>  (3.86) 


(3.87) 


;  .  ,  C  i'  -  30))] 

-  ,  (3.88) 


1  -  i  Tg  (o'  -  3tD)l 

<ct  ;i  of  aaturatlon  and  is  defined 


(3.89) 


Trm  1%.  (3.^)  WB  tixA 


1 

r  ^  Tg  ] 

(3.90) 

As  explained  In  section  0>1^  where  It  was  derived,  this  expression  Is 
an  approximation  In  the  limit  that  Is  less  than  the  theoretical 
maximum  value,  which  will  he  the  case  In  any  practical  situation.  A 
more  exact  expression  may  be  derived  from  Eq.  (3.^).  In  section  0  a 
more  detailed  discussion  of  saturation  will  be  given. 

By  using  (3.80)  and  (3.88)  ve  find  that  the  power  coupled  out  at 
the  hsrmonlc  Is 


(3a>) 


8®  f|  V 


Qfli  ♦  -  yx>f]  is?  - 


Pout^3“)  - 


2nQ 


2Q 

1  +  1  — ^  (3»  -  ®o) 

_ flOj  _ y 


8  f. 


Q^[l  -  1  Tg(0'  -  3®)] 


(3.91) 


The  above  expressions,  along  with  the  subsidiary  defining  conditions 
(3.19a)/  (3.^)/  (3.85),  (3.&r)/  (3.90)  give  the  value  of  the  harmonic 
power  coiqiled  out  of  the  cavity  In  terms  of  the  operating  conditions, 

(®  ,  B^)  ,  the  properties  of  the  cavity  (Q^  »  »  V  »  /  ^2  * 

and  of  the  quantum  system  (A  ,  ^  ,  Tg  ,  7)  .  It  then  conqpletely 

defines  the  generation  of  harmonic  power  due  to  the  nonlinearities 
in  the  quantum  system  In  tenos  of  the  familiar  linear  properlties  of 
the  system,  let  us  sljqillfy  l^hls  expression  and  put  It  In  a  form  which 
Is  more  mslly  undexmtsndsble  trem  m  physical  point  of  view. 
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It  v«  examine  the  denominator  of  (3.9l)>  using  the  definition  of 
loaded  Q  ,  1/^2  ■  *■  ^  the.  unperturbed  cavity 

Q  ,  we  see  that  it  is  of  the  form 


(3.92) 


which  is  that  of  a  cavity  (Q^)  coupled  to  the  extemsd  world  (Q  ) 

w  C 

with  an  additional  loss  mechanism  Q  _ .  .  The  value  cf  Q  , 

sample  sample 

describes  in  effect  the  "source  impedance"  of  the  sample  at  the  frequency 
3(0  .  The  exact  form  of  these  quantities  in  (3*91)  shows  that  they  are 
reactive  as  well  as  resistive.  !niat  arises  from  the  fact  that  both  the 
cavity  and  the  atom  are  resonant  and  will  have  a  reactive  component  if 
the  frequency  3(o  does  not  coincide  with  their  resonance.  Let  us  assume 
that  the  cavity  Is  tvined  by  varying  its  resonant  frequency  ca^  such  that 
the  reactive  or  imaginary  components  in  the  denominator  cancel.  Then  the 
dencxolnator  becomes 


(3.93) 


Taking  for  the  mcment  3-1  ,  i.e.,  no  saturation,  unity  filling  factor, 

f^  =  1  ,  and  3(0  =  >  the  effective  source  impedance  of  the  sample  is 

described  by  i/q  .  This  quantity  Is  equivalent  to  the  familiar  "magnetic 
Q  ,"  of  the  standard  maser  terminology  and  may  be  defined  In  the  usual 
manner  as 


1 


Q 


sample 


(power  lost  in  sample) _ 

2rt  (frequency)  (energy  stored  in.  cavity)- 


P-.  4. 

lost 


(3.94) 


If  the  system's  population  distribution  corresponds  to  positive  temperatiures 
(more  in  lower  level  than  in  upper),  then  (X^  -  ^^)  ^  0  implies  A>  0 
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laplles  >  0  and  hence  a  net  loss  or  attenuation  of  the  signal  at 

3a>  I  vhereas  if  populations  are  Inverted,  Q  Is  negative  and  the 

0 

systen  acts  as  a  negative  resistance.  The  posslhlUty  then  exists  for 
oscillations  to  occur  If  1/Q  becomes  sufficiently  negative  to  cancel 
tte  losses, 


(3.95) 


This,  of  course.  Is  the  familiar  condition  for  start  oscillations  in  a 
maser  oscillator. 

Let  us  now  cooiblne  the  cavity  Q  and  the  sample  Q  into  an  effective 
source  Q  defined  by 


8  f. 


Q  Q«  Q  [1  +  1^(0"  “  3to)^] 

source  0  s  2'  ' 


(3.96) 


Expre3slon  (3.93)  then  becanes 


(3.97) 


vhlch  Is  maximized  for  Q  =  Q  ;  under  these  conditions  ve  have 

e  source 


(  ^  1 

■  1 . . ,  .1 

2| 

Q  Q 

' 

source  e 

J 

source 


(3.98) 


max 


Referring  to  (3.91)  we  see  that  the  power  out  under  matched  conditions, 

Q  =  Q  ,  is  directly  proportional  to  Q  .If  the  sample 

Q  is  positive,  indicating  a  normal  population  and  hence  loss,  the  source 

Q  is  reduced,  lowering  the  power  out.  If  the  population  is  inverted, 

Q  <  0  ,  then  the  source  Q  will  Increase  and  the  output  power  will 
6 

increase  proportionally.  Ve  thus  see  that  the  question  of  population 

-  69  - 


inversion  le  iaiportant  in  determining  the  magnitude  of  the  emitted  power 

through  its  effect  on  the  lource  Q  but  is  not  required  in  order  to 

achieve  harmonic  geuerationi  From  a  qualitative  point  of  view  we  see 

that  the  system  bebavv^s  first. as  a  nonlinear  element,  generating  third 

harmonic  fields,  and  second  in  a  linear  manner,  either  amplifying  or 

attenuating  the  signal  depending  on  whether  it  is  in  an  Inverted  or  normal 

population  state.  Clearly,  the  latter  llneeur  process,  to  be  Important, 

reqxilres  that  an  allowed  transition  be  near  the  output  frequency  3io  as 

demonstrated  by  the  resonance  dencmdnator  in  the  sample  term.  In  the 

case  herein  considered  where  3a)  <»  C  this  is  satisfied. 

We  shall  not  pursue  the  question  of  population  inversion  any  further. 

It  will  be  assumed  that  the  system  is  in  a  normal  poj  datlon  state, 

Q  >  0  ,  and  that  the  external  Q  is  matched  to  the  source.  Then  using 
8 

the  definition  of  ,  Eq.  (3.91)  for  the  power  out  becomes 


'  2  source  1 

128rt[l  +  1^(0'  -  30))^}  (£1^  -  (D^)^ 


(3.99) 


Finally,  if  we  assume  that  conditions  sure  satisfied  such  that  ya  =  Si'  , 
there  is  no  saturation  so  8=1  ;  and  that  the  unloaded  Q  of  the 
cavity,  ,  is  much  larger  than  the  simple  Q  ,  ,  the  expression 

for  the  power  out  simplifies  to 


e) 


4  2 

^s^l 


128n 


(3.100) 


Before  discussing  in  detail  the  expressions  for  output  power  we  shall 
calculate  the  reaction  of  the  sample  to  the  fundamental  fields  and  then 
perform  the  calculations  for  the  case  where  o)  n.  0  ,  which  will  give 
similar  results  to  the  above. 
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h*  Rawttion  TeaaMml^ 

fiM  pemr  suffaicd  to  tbo  carlty  at  the  fundoaental  mqt  be  dieil]|^ted 
via  three  aain  «eebealau.  Flret,  there  la  the  cavity  loss  vhlch  la  due 
to  ^  finite  conductivity  hfvthe  valla  and  la  repreaented  by  the  unloaded 
quality  factor  of  the  cavity,  ^  .  Second,  there  nay  be  a  loaa  due  to 
abaorption  of  the  fundaBental  by  the  atonic  ayatem  vhen  there  la  an  atomic 
reamance  nearby.  For  the  caae  at  hand,  cd  «  o/3  ,  thla  loaa  vlU  be 

very  anall.  k>th  of  theae  loaa  mechanlans  are  first-order  In  nature,  the 
power  absorbed  being  proportional  to  the  energy  denjplty  .  They  may 
be  lumped  together  in  an  effective  firat  order  Q  ,  Q'  ^  •  Itae  third 

source  of  loss  la  third-order  in  nattire  and  is  that  power  going  via  the 
nonlinearity  into  the  third  order  process.  This  component  of  the  power 
supplied  by  the  fundamental  generator  goes  into  cavity  losses,  molecular 
absorptlcm,  and  output  power,  all  at  a  frequency  3(a  •  This  form  of 
loss,  reflected  by  the  nonlinearity,  is  familiar  from  the  theory  of 
parametric  amplifiers. 

To  evaluate  the  effective  Q  of  this  coupling,  we  may  proceed  in  a 
manner  similar  to  that  yielding  the  outpxxt  power  at  3(n  •  By  using 

(3»72),  (3.8l)  and  allowing  for  an  external  source,  we  have 


1  2HpQ[A  ^  -  X^^)] 


JIL 

2  2  * 
fl*"  -  0)  ♦ 


■t-  (external  source  term) 

(3.101) 


(•a)  (3) 

Upon  substituting  for  and  X^'  ,  rearranging  the  expression,  and 

dropping  higher  order  terms,  we  obtsdn 


f  1  ^fos  r  ^ 

®i  [(fl^  -  < 


il?l 


2(a^  -  o)"") 


(external  soiirce  term)  . 

(3.102) 


(3) 


The  effective  Q  of  the  third  order  losses,  Q 


Is  then  seen  to  be 


where  f  is  the  filling  factor  at  the  fuadaaentali  is  the  eaegpl# 

Q  ,  8  l8  the  laturatlon  paraneter.  ^  expected,  the 

magnitude  of  the  reflected  loss  la  d^^ndent  on  the  fundamental  field 
atrwagth  (8^^)  malting  the  loan,  l^/Q  ,  go  aa  .  The  maximum 

value  of  i/Q'^'  will  occur  for  8*1  ,  f  ■  1  ,  and  ^  j 

it  is 


(3.104) 


Since  represents  the  strength  of  the  perturbation  and  0  is  of 

the  order  of  ,  the  vinperturbed  Hamiltonian,  we  see  that  l/Q''^'  is 

of  the  order  tines  i/Q  .  For  practical  fields  [see  Eq.  (3*32)] 

0  8 

this  factor  will  be  small,  demonstrating  the  necessity  of  either  an  extremely 
high  cavity  Q  or  a  large  value  of  1/Q  in  order  to  obtain  any  sort  of 
efficiency;  otherwise  most  of  the  input  power  will  be  wasted  in  cavity 
losses.  Such  a  consideration  is  seen  to  limit  the  usefulness  of  this  as 
a  traveling  wave  device  where  the  field  strengths  will  be  quite  low  for 
reasonable  power  levels. 

We  shall  now  proceed  to  csdculate  the  expected  power  in  second  case, 
where  the  fundamental  fields  are  approximately  resonant  with  the  trsmsltion 
frequency  oc  aj  n 

2.  Case  II;  o)  *<  Q 
a.  Output  Power 

Having  gone  through  a  similar  calculation  for  the  other  case,  we  shall 
follow  the  same  procedure,  adding  cconents  where  there  are  differences. 

Prom  (3.68a}  and  (3*52)  we  have 


(3.105) 


•nA  froK  (3.68q)^  (3*9^)  «d  (3>99}>  we  have 


» »g  Pi 


oil  -  1  T^'  -  (0)]J 


(p.l06) 


hy  substituting  this  in  the  expression  for  the  generation  of  cavity  fields 
(3«72)>  asBianing  that  0  »  I/T^  ,  and  tuning  the  cavity  so  as  to  tune 
out  the  reactive  component  of  ,  we  have 

.j  fl  3  f,  8  1  8  »  f _  p? 

—  ♦ - - - 2  - - 2 - - *  (3.107) 

L^3  J  ^ 

where  f^  and  f.  are  the  filling  factors^  Q  is  defined  in  (3*89)/ 

X  s 

is  the  combination  of  unloaded  cavity  and  external  Q's  ,  and  8 
for  this  case  may  be  derived  from  (3*57)  and  is 


8  = 


1  ♦  Tglo"  -  mf 


(3.108) 


By  using  expression  (3.8o)  for  the  power  out,  we  obtain 


'’outO”) 


(3.109) 


2 

In  this  case,  because  of  the  factor  (1/^2)  >  contribution 

to  the  source  Q  from  the  sample  will  be  very  small  and  it  may  be 
neglected  in  caiiq;>arlson  to  the  cavity  Q  ,  for  reasonable  values  of 
the  latter.  Ihysiealily  this  is  due  to  the  fact  tl^t  the  output  ftsijuency 


-  73  - 


oorrMpondi  to  thzvf  tlm«  the  oaturftl  tronsltlon  fr«fu«ney  and  thKf  vlU 
be  negUgibXe  Belf-ebeorptlon,  In  contrast  to  the  prevlons  case  vhere  the 
output  oeeurz^  near  the  natural  transition  and  there  was  absorption. 
Optlsdsing  the  output  by  setting  *  Qq  >  bave 

- S' - —  (3.110) 


«  •  2^^  [1  +  TgCo"  -  a>f] 


Upon  assuming  0"  «  m  ,  6-1  ,  we  have  the  apprcslnate  expression 


rf  4%4^ 


2^3  ..3 


(3.111) 


b.  Reaction  at  the  Fundamental 

As  before^  there  will  be  three  main  contributions  to  the  loss  at  the 
fundamental:  the  cavity  losses,  the  first  order  absorption,  and  the  third 
order  losses  which  are  reflected  by  the  nonlinearity.  The  cavity  losses 
will  be  the  same.  For  this  case,  however,  the  first  order  losses  in  the 
quantum  system  will  be  much  larger  than  the  third  order  losses,  due  to  the 
proximity  of  the  natural  resonant  frequency  to  the  "p\mip"  freqiwncy.  The 
value  of  the  Q  due  to  the  first  order  losses  in  the  quantum  system  is 


1  fS  1 

0^  1  +  T2(ft''  -  00)^ 


(3.112) 


As  vlrtxially  all  the  third  harmonic  power  is  dissipated  in  the  cavity  and 
in  the  useful,  load,  the  third  order  Q  will  be  such  that  it  will  give  a 
loss  equal  to  twice  the  output  power; its  specific  value  will  net  interest 
lu  here. 

^^^When  the  output  is  matched,  Q,  ■  ^  power  dissipated  in  the 

cavity  equals  the  output  power  and  hence  the  total  power  dissipated  at  3a) 
equals  twice  the  output  power. 
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Our  iatuTMt  in  tb«  valuts  of  the  wal*eu3iar  Q's  ,  iKjth  at  tba 
fund  tat  al  aaft  at  tbo  hazamlc,  can  ta  aaaa  «h«  ooa  eonaidart  tlM 


Of  gattlag  the  aoat  haraonie  poaar  out  tor  a  slaaa  li^ait  powar  laval.  B* 
9rcl>34i  tt  (UffWant  In  ttm  two  easea. 

"For  tlia  flrat  caaa  oonaidarad,  o  it  q/3  ,  tba  quantua  agrataa  inmanta 
a  larga,  first  order  loss  at  tha  bamonle,  vMla  at  the  fisidaaeatal  there  la 
(»Ijr  a  Buch  saaller  third  order  loss.  An  Incraaaa  in  the  cavity  Q  at  3a 
is  seen  to  be  less  valuable  than  a  correspemding  Increase  in  the  Q  at  to 
since  the  effective  naxlam  Q  at  3o>  is  detensined  by  the  ssaqple  Q  , 

At  the  fundaaental,  however,  the  field  in  the  cavity,  ,  will  be 
determined  primarily  by  the  cavity  Q  ,  ,  as  the  third  order  losses 

of  the  quant lu  system  will  in  general  be  less.  In  this  liait,  4‘i  - 
and  the  output  power  for  a  given  input  power  will  vary  as  the  cube  of  the 
fundamental  cavity  Q  . 

On  the  other  hand,  when  o  «  0  ,  the  quantm  system  "loads"  the 

fundamental,  limiting  the  effective  input  Q  to  that  presented  by  the 
sample}  very  little  is  gained  by  Increasing  the  cavity  Q  beyond  this 
value.  Referring  to  Eq.  (3.110)  we  see,  however,  that  the  output  is 
proportional  to  the  cavity  Q  at  the  third  harmonic,  and  an  Increase 
in  its  value  will  bring  a  corresponding  increase  in  the  power  out. 

In  other  words,  if  we  consider  our  system  to  be  either  "transparent" 
or  "opaque,"  depending  on  whether  it  is  nonabsorbing  or  absoxblng,  we  find 
that  if  it  is  transparent  at  the  fundamental,  bamonlc  generation  is  enhanced 
by  providing  a  high  Q  resonator  at  the  fundamental,  while  if  it  is  opaque 
at  the  fundamental,  the  best  performance  will  occur  for  a  high  Q  at  the 
output.  If  the  system  were  transparent  at  both  frequencies,  then  improving 
the  Q  values  at  both  fundamental  and  harmonic  would  result  in  improved 
operation. 


It  should  be  noted  that  in  the  case  where  own  even  though  the  systmn 
absorbs  more  power  by  virtue  of  its  rescmance  near  o  this  power  is  first 
order  in  nature  and  does  not  help  at  all  in  contributing  to  the  generation 
of  harmonics.  It  is  to  be  resuasbered  that  only  the  third  order  losses 
(at  o)  contribute  to  the  harmonic  generatioa  process.  In  fact,  the  first 
order  losses  may  reduce  the  level  of  output  power  throtigh  the  phenomencan  of 
saturation  to  be  discussed  next. 
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0.  BkTmsim 


IQ  the  preceding  lectimte  iB«3tt(»  la  aade  of  the  phenoqeacm  of  aatuzwtloa* 

1 

In  the  cmtext  of  thia  problem  It  deaerlhes  the  tendency  of  the  poithUetlon 
density  of  the  xqpper  and  lover  levels  to  equaUse  tinder  the  Influehea  at 
strong  rf  fields.  A  parameter  8  has  been  defined  which  describes  this 
effect  by  relating  the  difference  in  population  density  under  dynastic 
conditions,  A  ,  to  the  difference  in  populaticm  density  under  conditions 

of  thermal  equilibrium,  df  .  When  the  population  densities  1?^  and 

0  ^ 

1^  have  their  equilibrium  values,  S  takes  Its  maxlmvim  value  of  unity. 

As  the  levels  tend  to  equalize,  A  tends  to  zero  and  hmtce,  S  tends  to 

zero. 

In  referring  to  the  equations  for  the  output  power  (3.99)  and  (3.110), 
we  see  that  in  both  cases  the  power  depends  on  the  square  of  S  .  Thus, 
as  the  levels  tend  to  equalize,  and  S  0  ,  the  power  output  falls  below 

the  nonsatiurated  value. 

When  a  system  tends  to  saturate,  this  implies  that  its  Internal  energy 
has  Increased  since  the  average  energy  of  the  system  (due  to  this  degree  of 
freedom)  Is 


E.  +  Ep  «l 

.  -t - £  -  _  A 

2  2 

where  E^  and  Eg  are  the  energies  of  the  two  levels.  This  Increase  In 

the  energy  will  occur  as  a  result  of  the  net  absorption  of  radiation  by 

the  syston.^^^  In  addition  to  Its  coupling  to  the  radiation  field  the 

system  is  presumed  to  be  In  contact  with  a  thermal  bath,  this  mechanism 

tending  to  restore  thermal  equilibrium,  with  the  systmn  giving  up  energy 

(2) 

to  the  bath  at  a  rate  given  by  the  inverse  of  the  quantity  When 

^^^It  is  implicitly  assumed  that  we  sure  considering  a  situation  where 
the  populations  of  the  two  levels  are  not  inverted,  corresponding  to  a 
positive  temperature. 

(21 

'  'The  quantity  T,  Is  defined  as  the  Inverse  of  the  rate  at  which 
the  system  thermallzes.'^' 
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the  net  ro\,e  of  energy  Input  exceeds  the  rate  at  which  the  system  can 
dispose  of  this  energy  (to  lattice  vibrations  In  a  solid  and  kinetic 
energy  or  other  d^prees  of  fi'eedom  In  a  gas)  then  the  system  begins  to 
saturate.  He  shsdl  now  investigate  this  effect  In  both  cases  considered 
and  see  how  It  limits  the  output  power. 

1.  Satvgatlon  When  aa  •»  n/3 

For  this  case  the  system  may  absorb  power  under  two  conditions.  There 

will  be  a  first  order  effect  corresponding  to  the  absorption  of  power  at 

the  fundamental  in  the  tail  of  the  absorption  curve.  For  aarrov  lines  one 

would  expect  this  to  be  small  since  the  applied  frequency  Is  approximately 

one-third  of  the  natiural  or  resonant  frequency.  The  power  absorbed  in  this 

2 

case  would  be  proportional  to  the  incident  fundamental  power  or  . 

The  second  term  In  the  saturation  Is  due  to  the  absorption  of  the  harmonic 

power  by  the  system  by  virtue  of  the  proximity  of  Its  resonance  Q '  to 

2 

the  output  frequency  3oi  .  This  power  will  be  proportional  to  or 
to  E^  .  We  thus  expect  the  saturation  parameter  to  contain  two  terms, 

^22  S 

one  proportional  to  E^  (or  and  the  other  proportional  to  E^ 

(or  p^.)  .  Referring  to  Eq.  (3*90),  we  see  this  Is  so. 

Let  us  now  calculate  the  field  strengths,  p^  ,  at  which  the  saturation 

takes  the  value  S  =  l/2  .  Assuming  3to  =  fl'  ,  the  first  order  term  wo\lLd 

cause  saturation  when 


Tg  n 


1 


or 


Tg 


while  the  third  order  term  would  cause  sat  mat! on  for 


(3.113) 


P,  - 


16  Tj^Tg 


(3.lii^) 
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For  the  assimptlon  w  Tg  >  valid  for  a  gae  at  microwave  fre9.uenclee, 
and  the  assimptlon  H  »  l/Tg  *  the  third-order  term  le  seen  to  be  the 
primary  source  of  saturation. 

At  the  cmset  of  satiuratlon  the  source  Q  ,  (3.96)/  tends  to 
and  the  expression  (3*99)  for  power  out  varies  as  .  For  the  predominant 
third  order  saturation,  Eq.  (3.99)  tends  asymptotically  toweu^s 


4 

(l  +  const.  E^) 


which  shows  that  the  power  out  actually  decreases  in  the  limit  of  large 


2.  Saturation  When  to  x» 

In  the  second  case,  o  «  fl  ,  the  strong  fundamental  fields,  near 
resonance  and  causing  first  order  absorption,  will  dominate  the  saturation. 
From  Eq.  (3.108)  the  condition  for  saturation  Is 


^1  “ 


(3.115) 


In  this  case,  under  the  conditions  of  saturation  Eq.  (3.110)  becomes 


e: 


2  2 

(l  +  const* E^) 


which  shows  that  even  under  these  conditions  the  power  out  Is  proportional 
to  the  power  In,  and  does  not  decrease  as  In  the  first  case.  Tn  this  .^.Bit 
the  expression  for  the  power  out  becomes 


out 


(3a))  =  — 


2  2 

3  Qq  f 2  E^  V 


„3  „2 
Q  n-'  T.  T- 
s  12 


(3.116) 
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3.  Comnenta  on  Sat\>ratlon 

Sy  coB^aring  the  saturation  levels  for  the  two  cases  of  Interest, 

(3<ll4)  and  (3«113}>  ve  have 

-  ki*  1?  t|  .  (3.U7) 

(ei)a  '  ^ 

As  JIT^  and  JlTg  are  usually  much  greater  than  unity,  we  see  that  the 
occurrence  of  saturation  comes  earlier  In  the  case  of  pumping  near  resonance. 
In  the  above  It  has  been  Implicitly  assumed  that  to  =  n '  or  oj  =  fl"  . 

A  further  look  at  the  exact  expressions  for  satiiration,  (3-90)  and  (3.108), 
shows  that  by  operation  off  resonance.  Si'  ,  Si"  -f  (a  ,  urie  effeetp  of 
saturation  can  be  reduced.  However,  (3.99)  and  (3.110)  show  that  corre¬ 
sponding  to  this  reduction  in  saturation  there  is  a  reduction  in  the  output 
power.  In  the  operation  of  a  device  based  on  these  principles  the  choice 
of  CD  relative  to  Si'  or  Si'  must  be  based  on  the  various  parameters 
of  the  problem,  such  as  power  available,  the  harmonic  power  desired  and 
the  properties  of  the  circuit. 

One  conment  which  can  be  made,  however,  is  that  in  all  the  expressions 
for  satiuration,  (3.113),  (3.11*^)  and  (3.115)»  a  speeding-up  of  the  relaxa¬ 
tion  process,  i.e.,  a  reduction  in  ,  will  result  in  an  increase  in 
the  value  of  at  which  saturation  occiursj  this  will  result  in  higher 

available  output  powers.  Hence,  in  any  material  which  is  to  be  used  for 
harmonic  generation  in  this  manner,  we  desire  a  short  or  spin-lattice 

relaxation  time. 

H.  DISCUSSION 

In  the  preceding  sections  of  this  chapter  the  mathematical  theory  of 
harmonic  generation  in  a  two-level  system  has  been  discussed.  Two  partlcu- 
leu:  cases  have  been  studied;  the  first,  in  which  the  frequency  of  the 
third  harmonic  output  is  near  the  natureO.  resonant  frequency;  the  second, 
where  the  applied  fundamental  frequency  is  near  the  natvu:al  transition. 

In  both  cases  the  quantum  system  behaves  as  a  nonlinear  element  generating 
a  component  of  the  macroscopic  polarization  at  a  frequency  three  times  that 
of  the  input.  This  time-varying  polarlzatlcxi  acts  as  a  source  for  the 
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electromagnetic  fields  which  nay  be  calculated  from  Maxwell's  equations. 

The  magnitude  of  the  polarisation  and  the  fields  generated  may  be  expressed 
entirely  In  terms  of  the  kr.cvn  parameters  of  the  quantum  system  and  Its 
surroundings;  no  phenomenological  nonlinearity  must  be  assumed.  Thus  It 
Is  seen  that  the  nonlinearity  Is  Intrinsic  to  the  quantum  system  and  may 
be  evaluated  In  terms  of  known  quantities.  Some  of  the  aspects  of  the 
results  obtained  will  now  be  discussed. 

1.  Resonant  vs  Honresonant  Behavior 

Although  only  two  particular  cases  have  been  considered,  namely  when 
the  applied  frequency  Is  such  that  m  «  n/3  or  to  n  ,  It  is  more  gener¬ 
ally  true  that  the  existence  of  the  nonlinearity  is  Independent  of  the 
frequency  oi  .  This  may  be  seen  from  the  fact  that  in  the  over-all  inter¬ 
action  the  total  energy  In  the  radiation  field  Is  constant  (three  fundamental 
photons  being  destroyed  and  one  at  three  times  the  frequency  being  created) 
and  hence  the  energy  of  the  quantum  system  Is  conserved. The  fact  that 
the  system  contains  discrete  natviral  frequencies  (in  this  case  only  one) 
manifests  Itself  in  the  theory  by  means  of  the  resonant  denominators,  which 
in  the  two  cases  considered  are  1  +  '  -  3a))  and  1  +  -  o)) 

A  general  theory  would  contain  two  such  denominators,  one  for  the  fundamental 
and  one  for  the  harmonic,  and  would  appear  as  the  product  of  the  two  above. 
The  two  cases  considered  eu*e  thus  only  special  cases  of  the  more  general 
result  in  which  we  have  chosen  the  frequency  to  minimize  one  or  the  other 
of  the  resonant  terms. 

The  important  thing  to  note  here  is  that  the  nonlinearity  is  Intrinsic 
to  the  system  and  independent  of  the  frequencies  applied.  The  relation 
between  the  operating  frequency  and  the  natural  transition  frequency  merely 
alters  the  magnitude  of  the  effect  through  the  resonance  denominators.  When 
the  operating  frequencies  are  all  far  from  quantum  transitlTona.  then  the 
system  looks  similar  to  a  nonlinear  reactive  element.  Wien-operatlng  near 
^  a  resonance,  the  system  looks  like  a  ccmiblnatlon  of  nonlinear  resistance 

^^^If  the  quantum  system's  energy  were  increased  or  decreased,  the 
change  would  have  to  correspond  to  the  difference  between  its  energy  eigen¬ 
values.  Such  a  change  would  place  restrictions  on  the  radiation  frequencies 
causing  the  transition,  namely  E  ±  a)^  *  (E^  -  E.  )/h  ,  where  +  indicates 
absorption,  -  indicates  emission,  and  E^  and  E^  are  the  final  and 
initial  energies  of  the  systaa. 
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and  reaetanee  but  ot  a  larger  nagnltuAe  due  to  the  reduction  of  the 
resonance  term.  In  the  miorowaye  region  /tite  aaylimnn  strenth  of  the 
nonlinearity  will  usually  be  saiall  enough  that  resonant  operation  will  be 
required.  In  general,  the  question  of  resonant  vs  nonresonant  operation 
will  depend  on  the  particular  aspects  of  the  overall  problem. 

Another  aspect  of  this  problem  la  the  qviestlon  of  saturation.  Under 
resonant  operation  the  system  will  absorb  energy,  thus  causing  saturation; 
for  nonresonamt  operation  the  saturation  Is  reduced  along  with  the  magnitude 
of  the  nonlinearity.  The  Importance  of  saturation  Is  again  a  function  of 
the  particular  problem,  being  most  severe  for  high  power  applications. 


2.  Dependence  on  the  Parameters  of  the  System 

The  dependence  of  the  ovitput  power,  Eqs.  (3.99)  and  (3.109),  on  the 
various  parameters  of  the  system  are  summarized  below: 

(a)  Dipole  moment.  l%e  output  power  varies  as  p.  ,  a  factor  of 

4  ,  2 

M  coming  from  1/Q  .  From  this  It  Is  clear  that  magnetic  dipole 

8 

transitions  will  give  negligible  effects  compared  to  electric  dipole 
transitions. 

(b)  Llnewldth  and  concentration.  The  ccxicentratlon  along  with  the 

llnewldth  enter  through  the  factor  Q  or,  s^.ternatlvely,  through  the 

8 

absorption  coefficient  y  .  They  occur  as  the  product  NTg  ,  which  for 
a  given  transition  Is  approximately  constant.  Thus  increasing  the  con¬ 
centration  N  ,  will  leave  Q  constant  while  Increasing  the  llnewldth, 

8 

and  conversely. 

The  quantity  Tg  also  occurs  In  the  resonance  denominators  as  well 
as  In  the  expressions  for  satvuratlon  (along  with  the  longitudinal  relaxa¬ 
tion  time  Tj^).  Here  a  broader  line  (smaller  Tg)  will  result  in  a  broader 
frequency  response  as  well  as  in  a  reduction  In  saturation. 

(c)  Population  difference  and  the  Boltzmann  factor.  In  both  cases 

the  output  power  is  proportional  to  the  square  of  the  popvilatlon  difference 

-2 

throxwh  its  dependence  on  Q  .  This  squeured  dependence  points  out  the 

s 

Important  fact  that  the  sign  of  A  or,  adtematlvely,  normal  vs  Inverted 
popiilatlon,  does  not  play  a  fundamental  role  In  the  nonlinear  process.  It 


'  'In  the  first  case  considered,  when  f 
cavity  Q  Is  large,  the  output  varies  as 


■  1  ,  f  T  “  1 

.  This  Is  a  special  case. 
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does,  however,  have  a  eiJbeidiary  affect  through  the  first  order  loading 
of  the  cavity  which  is  proportional  to  the  first  power  of  A  .  A  nrapaal 
population  distribution,  A  >  0  ,  causes  the  systeni  to  present  a  positive 
loss  to  the  cavity  and  lowers  the  Q  ;  an  Inverted  population  presents 
a  negative  loss  and  raises  the  Q  . 

There  exists  possibilities  of  using  inverted  populations  to  enhance 
the  overall  nonlinear  effect  but  they  will  not  be  given  here. 

The  nagnltude  of  A  may  be  found  from  Boltzmann  statistics,  which 
for  no  saturation  is 


1  v(kMc) 

A^A  - - ■  .  (3.118) 

2kT  Uo  T(  K) 

For  room  temperature  and  wavelengths  arovmd  1  cm,  this  is  a  feietor  of 
the  order  of  l/400  compared  to  a  maximum  value  of  unity.  Operation  at 
higher  frequencies  or  lower  temperatures  will  result  in  an  Increase  in 
A  ,  and  hence  an  Increase  in  the  output  power. 

(d)  Filling  factor.  The  output  varies  as  the  square  of  the  filling 
factor  fg  (3.73),  which  is  to  be  compared  to  a  lineeur  dependence  in  first 
order  processes.  This  places  a  strong  Importance  on  good  circuit  design 
by  which  fg  may  be  Increased.  The  other  filling  factor  f^  helps  de¬ 
termine  the  source  Q  (for  the  first  ceise),  and  should  be  minimized  con¬ 
sistent  with  a  maximization  of  fg  .  For  a  gas  completely  filling  a 
structure,  f^  =  1  . 

3.  Pynamlc  Shift  of  the  Hatural  Resonant  Frequency 

Cne  aspect  of  this  problem  which  does  not  appear  in  usual  formulations 
of  the  problem  is  the  effective  dettmlng  of  the  natural  frequency  of  the 
quantum  system.  In  the  two  cases  considered. the  new  frequencies  were 
denoted  n'  and  tl"  ,  and  are  given  by  (3.^2)  and  (3.53)*  This  detuning 
is  of  second  order  in  the  perturbation,  ^  ,  and  under  the  usual  assump¬ 

tion  of  small  pert\irbatlons  is  negligible.  In  the  field  of  magnetic  reso¬ 
nance  this  shift  is  known  as  the  Bloch-Slegert  shift.  In  terms  of  the  dc 


-  82  - 


field  Hq  and  the  rf  field 


It  Is  given  by 


(1) 


(3.119) 


where  ft'  Is  the  effective  frequenqy  and  ft  Is  the  frequency  corresponding 
to  the  unperturbed  eigenstates.  By  vising  values  of  dc  and  rf  fields  typical 
to  spectroscopy  or  maser  applications  this  shift  Is  of  the  order  of  one  part 

7 

In  10  or  less.  In  our  application  with  electric  dipole  transitions  and 
strong  fields  the  detvinlng  may  be  an  appreciable  fraction  of  the  natural 
frequency . 

At  first  such  a  concept  may  seem  strange,  as  the  usual  interpretation 
of  the  radiation  problem  is  that  of  quantum  Jumps  between  the  discrete 
energy  levels  of  the  system.  Such  a  picture,  howsver,  rests  on  the  assump¬ 
tion  that  the  perturbation  is  small  and  is  effective  only  in  causing  transi¬ 
tions  between  these  unperturbed  levexs.  Wlien  the  perturbation  is  strong 
the  only  rigorous  method  of  solution  is  to  go  back  and  solve  the  problem 
of  a  quantum  system  strongly  coupled  to  the  electromagnetic  field.  Under 
such  a  procedvure  it  should  not  seem  surprising  that  conditions  of  optimum 
interaction  would  differ  from  small  signal  theory.  For  the  electric  dipole 
case  such  a  shift  might  be  looked  upon  as  a  sort  of  dynamic  Stark  effect. 

From  the  point  of  view  of  the  problem  at  hand  the  variation  of  the 
effective  nat viral  frequency  with  power  level  Implies  that  given  a  fixed 
frequency  of  operation,  there  is  only  one  power  level  at  which  optimum 
interaction  will  occur.  Vfhether  or  not  this  variation  of  the  natvural 
freqviency  is  significant  depends  on  the  linewidth  of  the  transition.  The 
effect  will  be  less  for  broader  lines. 

^^^Thls  is  precisely  the  same  relation  as  Eq.  (3.53).  It  may  be 
derived  as  in  the  text  or  by  Bloch's  original  arguments.  Another  simple 
method  is  to  take  the  linearly  polarized  field  and  separate  it  into  two 
circularly  polarized  components.  By  diagonalizing  the  Hamlltonlsui  contain¬ 
ing  the  dc  field  and  the  counter-rotating  term  and  averaging  the  latter, 

Eq.  (3.119)  is  derived  as  the  difference  between  the  new  eigenvalues. 
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4.  An  Example;  Harmonic  Generation  in  QaBeoua  HHj 

Me  shall  now  consider  an  example  using  the  Inversion  transition  in 

. 

a.  !nie  Pressure  Dei>endence  of  the  Spectrum 

The  inversion  spectrum  of  annonla  is  made  up  of  many  lines,  each 

corresponding  to  a  different  rotational  mode  of  the  molecule.  At  low 

pressures  these  lines  are  resolved,  each  having  its  own  value  of  absorption 

coefficient.  As  the  pressure  is  Increased  each  of  these  lines  broadens 

with  half -width  of  30  Mc/mm  Hg,  while  preserving  the  magnitude  of  its 

absorption  coefficient.  When  these  lines  begin  to  overlap,  the  overall 

absorption  coefficient  at  a  given  frequency  will  be  the  sum  of  the  indl- 

vidiMil  contributions .  At  high  pressures  the  many  lines  will  be  unresolved 

and  will  appear  as  a  single  broad  line  with  an  absorption  coefficient 

which  is  nearly  the  sum  of  the  Individual  absorption  coefficients.  Thus 

Q  ,  and  hence  the  level  of  power  out,  should  Increase  with  pressure  up 
s 

to  a  point  where  the  llnewidth  of  an  individual  component  becomes  of  the 
order  of  the  si>aclng  between  the  lines.  For  ammonia,  this  occurs  at  a 
pressure  of  about  200  mm  Hg,  corresponding  to  a  half  width  of  6  kMc. 

b .  Saturation 

For  most  gases  with  transitions  in  the  microwave  range  we  have 
T^  =  Tg  Since  NTg  is  a  constant  and  N  «  pressure,  increasing 

the  pressure  will  reduce  T^^  and  Tg  and  hence  will  reduce  saturation. 

The  effects  of  saturation  should  then  be  most  predominant  at  low  pressures. 

By  taking  ah  operating  pressure  of  300  mm  Hg  where  the  line  has  nearly 
its  maximum  absorption  coefficient, ^  7=Txl0^cm^  and  by  using 
the  following  data  for  NH^  : 

^  =  1.47  X  10”^^  esu 

fi  =  2jt  X  2if  a4c 

£a/  =  l/2«  Tg  =  30  Mc/mm  Hg  , 


Try 


See  Townes  and  Schawlow, 


this. 


Ref.  (59),  PP.  352, 


for  a  discussion  of 
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ve  find 


^  -  l.li  X  10"^ 


and 


(n  Tg)  *  2.67 


The  relations  for  output  power  (3.100)  emd  (3*111)  are 


f2 

Poy^|. (watts)  =  4.7  X  10"^°  —  E^(esu) 

^1 


where  the  third  hannonic  frequency,  3co  ,  is'  near  2k  IMc,  and 


Po^t(watts)  =  3.1  X  10’^°  f|  E^(e8u) 


where  3a) »  72  Idle  and  *  5,000  .  Plots  of  the  theoretical  output 

power  per  c.c.  of  cavity  volume  for  the  two  cases  are  plotted  in  Fig. 

3.2,  where  the  effects  of  saturation  have  been  included.  In  Fig.  3.3 

the  saturation  parameters  S  ,  (3.90)  and  (3*108)  are  plotted  for  the 

two  cases  of  Interest.  For  both  plots  the  quantities  are  plotted  against 
2 

Et  which  is  proportional  to  the  power  supplied  at  the  fundamental. 

A  similar  calculation  using  HCN  ,  (7=9xl0^cm  ,^=3x10  , 

n  =  2jt  X  88.6  ld4c)  gives 


•’out 


0.55  X  10*^^  B^(e8u) 


% 


-  85  - 


-  87  - 


I.  CONCWSIONS 


In  this  chapter  a  detailed  calculation  of  the  process  of  third  harmonic 
generation  In  a  two-level  syet«s  has  been  performed.  The  value  of  the  non¬ 
linearity  has  been  related  to  the  constants  of  the  quantvmt  system  and  Its 
surroundings,  which  here  Is  a  cavity.  The  effects  of  saturation  and  the 
losses  In  the  form  of  Q's  have  been  derived  from  the  same  formalism. 

The  dependence  of  the  nonlinearity  on  the  parameters  of  the  system  has  been 
discussed  and  an  example  has  been  given  using  the  well  known  inversion 
transition  In  ammonia. 

The  particular  systoa  studied^  (the  two-level  system)  Is  the  simplest 
of  all  quantum  mechanlceO.  systems  and  Is  clearly  only  an  approximation  to 
any  reel  situation.  It  does  demonstrate,  however,  that  even  such  a  simple 
system  can  show  nonlinear  properties.  Other  quantum  systems  containing 
three  or  more  levels  will  display  similar  nonlinear  properties  whose  magnitude 
and  frequency  range  will  depend  on  the  particular  system  Involved.  For 
example,  in  Chapter  V,  a  three-level  system  Is  considered  as  a  second 
harmonic  generator.  Other  possibilities  exist  for  a  two-level  system; 
second  harmonic  generation  may  be  performed  by  the  use  of  a  dc  blas,^^ 
and  a  form  of  parametric  amplification  Is  possible. We  shall  next 
consider  the  parametric  process. 


See  Chapter  IV 
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chafur  IV 


PARAMETRIC  EFlfiCTS  IN  A  TWO-LEVEL 
ELECTRIC  DIPOLE  SYSimt 


A.  IHTROIXJCTION 

In  the  preceding  chapter  it  was  shown  that  a  two-level  electric  dipole 
system  behaved  In  a  nonlinear  manner  and  Its  use  as  a  third  harmonic  gener¬ 
ator  was  calculated.  Here  another  application  of  nonlinear  phenomena, 
namely  pareunetrlc  amplification  or  oscillation,  Is  considered  using  the 
same  two-level  system.  This  process  has  one  very  Important  difference 
from  harmonic  generation  processes.  For  the  latter  there  exists  an  out¬ 
put  (however  small)  for  any  level  of  Input.  To  obtain  parametric  oscilla¬ 
tion  (sometimes  called  subfrequency  oscillation)  the  level  of  the  applied 
"pump"  signal  must  exceed  some  critical  value  determined  by  the  nonllne- 
eorlty  and  the  circuit  parameters.  This  value  Is  known  as  the  threshold. 

We  shall  calculate  the  threshold  for  such  a  process,  relating  It  to  the 
circuit  Q  and  to  the  parameters  of  the  two-level  system. 

The  quantum  mechanical  system  to  be  considered  consists  of  two  levels 
with  energies  E^  and  E^  ,  where  the  natural  frequency  Is  defined  as 
n  *  (Eg  -'<S^)/ll  .  The  Interaction  of  the  system  with  the  radiation  fields 
will  be  taken  to  be  via  an  Induced  electric  dipole.  In  the  energy  repre¬ 
sentation  the  perturbing  Hamiltonian  will  then  contain  no  diagonal  compo¬ 
nents,  as  a  result  of  which  higher  order  processes  connecting  the  two 
states  must  contain  an  odd  number  of  photons. 

There  are  two  Interesting  types  of  parametric  processes  possible 
using  such  a  system.  In  the  first,  four  frequencies,  oc^  , 
are  present.  Energy  conservation  requires  (Oj^i-ec^  +  ec^-i'a^-O  , 
where  a>  >  0  for  absorption  and  co  <  0  for  emission.  By  setting 

and  providing  cavity  resonances 
at  0^  and  ,  a  pazametrlc  amplifier  Is  possible  where 


20^ 


®s  *®1 


(4.1) 
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In  lueh  a  ayatem  there  will  be  eleetrooagnetle  flelda  at  ,  o)^  and 

CD  .  Thla  la  ahown  achenatlcally  In  Fig.  it'.l. 

8 

In  the  aeeond  proeeaa  the  atoaile  reaonanee  acta  aa  the  Idler, 

CD,  B  0  ,  and  no  eleetroiaagnetlc  field  need  be  preaent  at  cd,  .  Again 


by  setting 


SL 


the  frequency  condition  la 


2®_ 


CD  +  0 

8 


(4.2) 


This  is  essentially  a  Raman  type  process  in  which  the  excitation  Is  by  two 
photons,  and  Is  ahown  In  Fig.  ^.2.  Here  a  cavity  resonance  need  only  be 
provided  at  cd  In  order  to  obtain  the  pareunetrlc  process.  It  Is  this 
latter  process  which  will  now  be  described. 


FIG.  ^.1--Energy  level  diagram  for 
the  first  type  of  para¬ 
metric  process; 

2od  ■  CD  +  0), 
p  s  1 


FIG.  U. 2- -Energy  level  diagram  for 
the  second  type  of  para¬ 
metric  process; 

2co^  »  CD.  +  0  . 

p  s 


B.  THE  EQUATIONS  OF  MUTION 

The  equations  of  motion  for  em  Induced  electric  dipole  may  be  derived 
in  the  following  manner: 

^^^Theae  equatlopa  are  derived  In  a  report  by  B.  T.  Jaynes. Their 
equivalence  to  the  density  matrix  formulation  la  demonstrated  In  Appendix  C. 
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1.  A«suine  only  two  energy  levels  are  involved^  so  that  the  state  of 
the  syste®  can  he  represented  sis 

-  a^(t)lu^>  +a2(t)|u^  ,  (^.3) 

are  the  eigenfunctions  of  the  unperturbed  system. 

2.  The  Hamiltonian  le 


♦ 

where  |  u^)  and  |  u^ 


where  31^  la  the  unperturbed  Hamiltonian,  Is  the  dipole  operator 

and  E  Is  the  radiation  electric  field. 

3.  In  matrix  notation,  these  sure  defined  as 


qp 


^^121 


;  M 


Eg 


where  u^g  Is  the  matrix  element  of  between  the  two  states,  smd 

and  Eg  su'e  the  eigenvalues  of  the  unperturbed  Hamiltonian. 

4.  The  dipole  moment  Is  given  by  the  effective  value  of  the  dipole 
operator,  P  » 

5.  Time  derivatives  are  calculated  from  the  ccumiutoitor  with  the 
Hamiltonian, 


^  -  IS  *  ('^•5) 


6.  From  Sqs.  (4.4)  and  (4.^)  the  following  equations  of  motion  are 
obtained: 


and  W  la  proportional  to  th«  energy  atored  In  the  dipole  ayatem,  referred 
to  the  equal  population  condition  aa  the  zero  point  of  energy. 

Equation  (4.6)  deacrlbea  the  dipole  moaent,  P  ,  aa  an  harmonic 
oscillator  driven  by  the  electric  field  through  a  coupling  constant  pro¬ 
portional  to  W  .  If  the  populations  are  equal;  the  coupling  to  the 
electric  field  la  zerO;  which  means  that  the  dipole  does  not  "see"  the 
radiation  field.  The  maximum  coupling  to  the  field  occurs  when  the 
population  difference  Is  a  maximum. 

Equation  (4.7)  equates  power  absorbed  to  the  rate  of  change  of  stored 
*  • 
energy;  P  la  equivalent  to  a  current  density  so  the  BP  Is  the  power 

delivered  to  the  dlpoleu:  system  by  the  radiation  field;  W  ,  In  turn.  Is 

the  rate  of  change  of  energy  stored  In  the  dipolar  system.  To  accovint 

for  loss  mechsuilsms  that  might  exist,  Eqs.  (4.6)  and  (4.7)  are  modified 

to 


(4.9) 
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idMr*  Tq  It  ttot  rtltnUtn  tiat  titetltitl  with  tte  daeay  of  tba  dlpala 
aonant,  aad  la  tfaa  ralaxation  tlM  for  tba  daeajr  of  atorad  aoargi'* 
n»  tf  la  tba  ataadx-atata  raltia  for  the  unperturbed  anei^,  ablch,  for 
a  afvtw  la  ttaeraal  afilLllbrluB,  la  glvaa  hr  the  BoltzBaan  dlatrlbutlon. 

lq.\]atloBt  (t.S)  an4  (V.9)  lavolve  the  three  varlablea  V  ,  P  and 
1  ,  aad  ao  an  additional  equation  la  neeeaaary  for  a  unique  aolutioaf. 
Proa  Slater 'a  aorauLl  node  expanalon  for  the  flelda  In  a  reamiaat  eavltjr, 
Bq.  (4.10)  can  be  obtained: 


B 


iwBa^fPV^ 


2 

“a 


CD  CD 


(D  1 


(4.10) 


where  cd 
a 


reaonant  frequency  for  the  unperturbed  cavity, 
loaded  quall'^  factor  for  the  cavlty> 


f  •  the  filling  factor, 

V  -  cavity  volume,  and 

B  *  number  of  Interacting  dlpolea  per  unit  volume. 

It  la  MBumed,  In  Eq.  (4.10),  that  one  normal  mode  field  predomlnatea, 
and  cga  unite  cure  uaed. 


C.  PARAMBIRIC  ANPLIFICAHOR  OR  OSCILLATIOB 

A  ateady-atate  aolutlon  to  Bqa.  (4.8),  (4.9)  and  (4.10)  can  be  obtained 
by  letting 
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where  P  ^  ■  P2  and  I  ^  •  K  ,  aa  required  to  give  real  valuea  for  1 
-a  Of  -a  cx  ' 

and  P  .  .  Each  value  of  a  correapondi  to  one  frequency  present. 

ly  cosiblnlng  Iqs.  (4.8)  and  (4.9)  with  the  ejqpansions  above  and  using 
harmonic  balance ^  we  get 


where  ■  moment  at  frequency  -  (o)^  +  ‘o.)  >  and 


JT 

-  - 

0  e 


-Iv-A 


(4.12) 


The  term  is  proportional  to  the  average  value  of  the  population 
difference j  and  Is  different  from  W  because  of  the  saturation  caused 
by  the  presence  of  the  fields. 

Considering  now  the  case  of  a  cavity  containing  active  dipoles  and 
resonant  at  two  frequencies  oo^  emd  ,  let  us  define  another  frequency, 

^  "tbe  relationship 


2©  ■  ©  +  ©.  .  (4.13) 

p  SI 

The  quantities  ©  ,  ©  and  ©.  ,  are,  respectively,  the  puop,  signal 

P  8  X 

and  Idler  freqviencles  of  the  paraunetrio  system.  Since  the  cavity  does 
not  resonate  at  ©^  ,  then  E^  is  zero,  but  P  has  a  eoDq>onent  at 

that  frequency. 

The  anqplltude  of  the  punq)  field,  E^  ,  required  for  parametric 
oscillation  Is  reduced  as  the  pump  frequency  approaches  il  .  For  this 
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reason,  let 


-  q(1  +  c) 

•  n(l  +  X) 

ffi  ■  fl(l  +  26  -  X)  ,  (^.1^) 

S 

where  6  ,  X  «  1  .  Sie  term  X  le  Introduced  to  represent  the  shift 
of  the  natural  frequency  a  due  to  the  presence  of  the  strong  punq? 
fields.  From  Eqs.  (^,11)  and  (U.12),  we  get 


6  n 


(4.15) 


—k 

2£(€  -  X)/  ‘ 
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In  the  derivation  of  Bqa.  (4.16),  (4.17)  and  (4.18)  the  follovlng 
aesunptlons  have  been  made: 

(a)  c  -  X  » 

(b)  -  Tg 

(c)  The  pump  anplltudee  are  imich  more  Intense  than  the  signal  emd 
Idler  angtlltudes  so  that  only  first  order  terms  In  signal  and  Idler  have 
been  Inclxided. 

Equation  (4.l8)  expresses  the  satviratlon  effect.  For  low  values  of 

pump  field,  ve  have  »  V  ,  and  as  the  amplitude  of  the  punqp  field 

u  o 

increases,  approaches  zero,  which  means  that  the  popipatlons  tend  to 

equalize.  A  normal  population  distribution  corresponds  to  a  negative 
value  for  ,  and  similarly  a  populatiwu  inversion  means  that  Is 

positive.  For  a  normal  population,  the  range  of  Wq  Is  0  >  Wq  >  .' 

There  Is  a  relationship  between  6^  ,  e  and  X  so  as  to  optimize 

the  parametric  effect.  This  relationship  c£in  be  obtained  by  considering 
the  determinant,  A  ,  formed  from  the  coefficients  of  the  dipole  moments 
and  P  .  ; 

8 


A  -  4fi^  [b^  -  X(X  -  2c)]  + 


J4fi  j(€  -  X)*^  +  e 


€  -  X 


(4.18) 


The  magnitude  of  A  becomes  a  minimum  near  the  condition  for  which 

6^  »  X(X  -  2c)  ,  (4.19) 


which  means  that  Eq.  (4.19)  expresses  the  condition  for  optimizing  the 

2 

dipole  moments.  Using  the  value  for  B  given  by  Eq.  (4.19),  P_  Is 

6 

given  by 
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Iquating  the  laeglnazy  teme  that  reeult  froa  auhatltutlng  the  abofe 
expreeelon  for  into  la.  the  eirc\ilt  equation,  we  obtain  the 
oeclllation  condition  in  terse  of  the  loaded  Q  of  the  reeonator: 


1 


(e  -  X)^ 

1 


X(X  -  2e)‘ 
1  + - 2 - 

€ 

(4.21) 


where 


1 

H 


fi  /  2tf  k(l 
i  /  a 


(4.22) 


Nt  *  population  of  the  lower  state 
J 

»  energy  stored  In  the  cavity  "  • 


The  derivation  of  Eq.  (4.21)  aseuaes  that  P  and  E  are  uniform 
over  the  interaction  region.  A  negative  sign  is  used  in.Eq.  (4.22)  to 
maka  H  positive,  SIS  <  0  for  the  normal  population  distribution. 

The  real  terms  that  appear  from  the  substitution  of  Eq.  (4.20)  into 
Eq.  (4.10)  are  cancelled  by  choosing  the  appropriate  value  for  ,  the 
resonant  frequency  of  the  cavity. 

By  lotting  X  ■  V«  »  *4*  (4-21)  becomes 


H  [(1  -  xf  +  l]  [l  -  x] 


(4.2)3) 


Prom  Eq.  (4.19)  it  is  seen  that  6^/e^  •  x(x  -  2)  ,  so  that  the 

permissible  range  for  x  is  x  >  2  or  x  <  0  .  Since  the  right  hand 
side  of  Eq.  (4.23)  mist  be  positive,  the  region  of  interest  for  x  is 
X  <  0  ,  which  msauM  that  c^  and  a>^  are  both  eitl»r  larger  or  smaller 
than  fl  . 
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Iftifttlon  (k,23i)  war  Is*  vrlttw  la  ttxw  of  tho  on-rotMMaeo  aaMtur> 
ated  absorption  eosffleiwt  glwn  in  Appendix  D: 


vhex*  e  is  the  velocity  of  light.  Assuming  a  filling  factor  of  unl^ 
for  the  resonant  cavity,  we  then  find  that 


(^^.25) 


where  X  is  the  wavelength.  Combining  Bias.  (4.23)  and  (4.23)  results  in 


rXQj^  [(1  -  xf  +  l]  [l  -  x] 


f(x)  . 


(4.26) 


The  function  f(x)  has  a  maximum  value  approximately  equal  to  0.13 
which  means  that  we  must  have  4x/7XQ^  <  0.13  for  oscillations  to  occur. 
Tor  a  rotational  mode  of  oscillation  In  a  gas  at  room  temperature  at 
X  »  1  mm  ,  a  strong  absorption  line  might  have  7  •  lO"'*’  cm''*’  ,  which 

means  that  >  (4i(/l.3)  x  10^  .  If  Is  specified,  then  x  is 
given  by  Eq.  (4.26).  As  em  example,  with  ■  1.2  x  10^  ,  we  find  that 

X  ~  -  1  In  terms  of  e  ,  this  means 

-  0(1  +  e) 

-  n(l  *  «) 

cD_  -  n(l  +  3€) 


Vx(x  -  2)'  Id  - 


l« 
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For  *  10"^®  ,  we  have  -  lO^e  kv/cm  .  With  c  -  0.01  ,  then  we 
find  Ep  "  10  kv/cm  and  f g  -  fp  -  6  kMc/s  Is  the  difference  between 
the  signal  and  pump  frequencies. 

1).  CONCLUSIONS 

The  preceding  discussion  has  indicated  a  manner  in  which  a  two-level 
system  can  function  as  a  parametric  amplifier.  It  differs  from  usxial 
parametric  amplifiers  In  that  the  quantum  system  provides  the  idler^^^  and 
th"  pump  frequency  equals  half  of  the  sum  of  the  signal  and  idler  frequenr 

clue.'  ' 

By  providing,  a  dc  bias,  operation  can  occur  where  to  »  oj  +  cu 

P  8  i 

other  ochemes  using  more  than  two  levels  offer  the  possibility  of  parametric 
procoBseo  which  may  prove  more  efficient  than  the  one  here  presented. 


^^^In  this  respect  It  is  similar  to  Suhl'a  ferrlmagnetu*  amplifier, 
which  can  use  a  mode  of  the  sample  for  the  idler. 

(2) 

'  'Subharmonic  pumping  of  a  parametric  amplifier  lias  been  uv*i'f  Mured 
b,y  Mortenson.  63 
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CHAprat  V 


!rHE  THEORY  OF  SECOHD  HABMOHIC  GEHERATIOE  IN  A  THBEE*LEFEL  SISTEH 

In  this  chapter  ve  shall  be  coneemed  vith  the  Interaction  of 
radiation  vlth  a  quantum  mechanical  system  consisting  of  three  levels. 

The  system  may  have  only  three  levels  In  the  frequency  region  of  Interest, 
as  In  the  case  of  an  S  «  1  spin  system,  or  If  there  are  more  than  three 
levels  ve  may  consider  those  most  strongly  coupled  to  the  fields.  Such  a 
three-level  system  can  act  as  a  three-frequency  mixer,  a  special  case 
being  second  harmonic  generation. 

Using  the  density  matrix  approach,  the  source  of  the  nonlinear  effect 
Is  pointed  out  and  the  problem  of  harmonic  generation  Is  considered  In 
detail.  The  dependence  of  the  nonlinearity  on  the  parameters  of  the 
system  Is  considered. 

A.  EQUATIOWS  OP  MOTION 

Let  the  system  under  consideration  consist  of  three  energy  levels 
E^  ,  Eg  ,  and  E^  ,  In  order  of  Increasing  energy  and  define  the 
natiiral  frequencies  as  In  Fig.  ^.1. 


31 


E. 


21 

J_ 


PIO.  5.1— Energy  level  diagram  for  a  three-level  system. 
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fbe  HavUtaoiaa  MQr  te  tttiMo  u 

+  ,  (5.2) 

/ 

vhere  3i  is^the  perturbatloi  resxiltlqg  from  the  explication  of  external 

fields  aM  Is  the  matrix  element  of  the  perturbation  taken  between 

levels  1  and  i  .  The  energy  representation  Is  used  aM  It  Is  assumed 

that  there  exist  n<»izero  matrix  elements  connecting  all  of  the  levels. 

Ve  may  write  the  equatl(xis  of  motlm  of  the  system  using  the  density 

matrix  formulation  In  the  standard  manner.  We  shall  etsstime  that  the 

longitudinal  relaxation  times  are  equal  while  the  transverse  times  may 

be  arbitrary.  The  eussunptlon  of  equal  longitudinal  times  simplifies  the 

eq'jations  for  the  diagonal  components  and  the  resulting  expressions  for 

the  saturation  effects.  Such  an  assumption  is  physically  Justified,  as 

64  . 

it  is  usually  the  ease  that  such  times  are  equal.  Further,  we  are  not 
here  Interested  in  effects  which  depend  in  any  way  on  such  discrepancies 
in  relaxation  times;  their  effect  would  be  merely  to  alter  the  details 
of  the  operation  and  would  In  no  way  affect  the  principles  Involved. 

The  equations  for  the  three-level  system  are: 


) 


^21  ^3  ^31  ‘  *^12  ^1  ■  '^13  ^31] 

“  [^21  ^12  *  ^3  ^32  "  ^21  ^  ■  *^23  ^32]  ^5. 2b) 
“  [^31  ^13  ^32  ^23  ■  ^"3 "  ^32  ^23] 


Try 


See  Appendix 


A. 
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aad 


^12 

12 

^23  ^  ^23 

^13  ^  ^13 


^  “21^12 "  B 1 

[(P22  ’  *  ^32  ^3  •  P13  ^32] 

1  ^32p23  "  H  1 

[{P33  -  p22^^"3  *  ^13  ^21  "  ^21  ^3^ 

^  ^31^13  “  'S  ■ 

[(P33  -  Pii)^3  ■*■^23^"  ^12  ^23! 

(5^2d) 

(5.2e) 

(5.2f) 


Let  us  examine  these  equations.  The  first  three  describe  the  behavior 
of  the  diagonal  conq)onents  which  are  essentially  the  probability  of 
occupancy  of  the  three  levels.  In  the  absence  of  an  external  perturbatl<»i^ 
31  «  0  ,  these  quantities  take  on  their  equilibrium  values.  If  31  0  , 

then  they  are  modified,  giving  rise  to  saturation  effects.  The  off- 
diagonal  equations  describe  the  coupling  of  the  states,  and  It  Is  these 
which  describe  the  absorption  and  emission  of  radiation.  Let  us  specif¬ 
ically  consider  the  equation  for  the  coupling  of  the  2-3  transition.  The 
first  term  on  the  right  hand  side  of  the  ei^presslon  for  Pg^  gives  the 
direct  coupling  due  to  the  pertxirbatlon  3ig^  which,  to  be  effective,  must 
be  almost  rescmant  with  the  frequency  .  In  addition  to  the  direct 

ccvpllng  there  are  two  additional  terms  which  give  rise  to  indirect 
coupling.  The  first,  3ig^  ,  shows  that  If  there  Is  coupling  of  the 
1-3  transition,  ^  ,  then  the  perturbation  3ig^  will  mix  with 

this  to  give  a  resultant  coupling  of  the  2-3  transition.  The  last  terrx 
Is  of  the  same  nature.  Since  p^^  will  be  seen  to  be  proportional  to 
3^^  ,  this  coupling  term  Is  proportional  to  3ig^  3^^  and  Is  seen  to  be 

of  second  order.  Such  a  process  involves  two  photons  and  Is  nonlinear 
In  character.  The  other  two  equations  are  of  the  same  nature,  displaying 
both  linear  and  second  order  coupling. 

We  may  now  qualitatively  understand  the  behavior  of  this  system  as 
a  frequency  mixer.  Let  us  suppose  that  we  apply  two  radiation  fields  at 
freq\iencles  m  and  ov,  .  Let  us  further  postulate  that  they  are  ap- 
proximately  resonant  with  the  1-3  and  1-2  transitions,  respectively.' 


'  'In  general,  any  two  freqTiencles  can  be  mixed,  but  unless  they  are 
approximately  rescmant  with  the  natural  frequencies  the  effect  will  be 
small. 
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These  pertuihatlcns  will  then«  In  a  linear  fashion,  drive  the  system  saA 
give  nonzero  values  to  the  ooifjpoaents  and  .  These  will  then 

nix  with  the  perturhatlocs  and  3^^  and  give  a  driving  term  for 

at  a  frwffmney  a>^  ■  oi^  -  .  As  shown  In  Chapter  II,  we  nay  cal¬ 

culate  the  component  of  the  dipole  vZKMtent  at  the  fre^Luency  (o^  ,  and 
since  It  Is  proportional  to  Pg^  we  see  that  the  fields  at  frequencies 

(Cl  Slid  (Up  will  mix  to  generate  a  cooponent  of  the  dipole  laoment  at  the 

^  fl ) 

difference  frequency.'  ^  If  we  now  provide  a  cavity  resonance  at  this 

(2) 

frequency,  we  can  extract  power  axid  the  system  acts  as  a  mixer.'  ' 

A  similar  type  of  parametric  pr(x:e8s  has  been  considered  by 

33 

Anderson.  There  he  applies  a  pump  field  resonant  with  and  allows 

for  cavity  resonances  at  (c^  m  Og^  and  (Ug  m  0^2  +  ®2  “  ®3l  ' 

Above  a  given  threshold  punq>  power  level  determined  by  the  Q's  of  the 
cavities  one  can  achieve  a  form  of  parametric  anpllflcatlon  or  oscillation 
at  the  frequencies  oi,  and  (u_  .  It  Is  similar  In  operation  to  Sxihl's 

ferromagnetic  amplifier.  * 

We  shall  now  consider  the  specific  case  of  second  harmonic  generation 
using  the  same  three-level  system. 

B.  BES(»fANT  HAFMOinC  GENERATION 

We  shall  start  by  assuming  that  we  have  a  three-level  system  in 
which  the  intermediate  level  with  energy  E„  lies  approximately  midway 
between  the  Outer  levels  which  have  energies  E^  and  E^  .  We  shall  assume 
that  there  exist  matrix  elements  connecting  all  the  levels  and  that  the  longi 
tudinal  relaxation  times  are  equal.  We  shall  assume  magnetic  dipole  inter¬ 
actions  and  let  there  be  a  strong  field  applied  at  the  fundamental  frequency 

H((d)  «  cos  ((ut  -  ^)  , 


^  ^It  should  be  noted  that  there  will  eilso  be  generated  a  comgponent 
of  dipole  moment  at  the  sum  frequency  but  because  of  the  resonant  nature 
of  the  quantum  system  this  will  be  very  small,  as  it  does  not  corre8p(xid 
to  any  of  the  natural  transition  frequencies  of  the  system. 

(2) 

'The  magnltxide  of  the  power  out  will  be  pr(9ortional  to  the 
product  of  the  power  at  frequencies  (Ul  and  (Ug  and  »?ill  be  nonzero  for 
any  finite  value  of  this  product. 
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and  a  weaker  aeeood  baziaoDle  field  whleh  Is  gmerated  lay  tbe  nooUitear 
process,  given  by 

H(2a>)  ■  Hg  cos  (2tDt)  . 


Ve  need  only  consider  perturbations  which  are  almost  resonant  with  the 
various  transitions,  as  the  effect  of  nonresonant  perturbatic»is  Is  very 
small. We  may  then  write  our  perturbing  Hamiltonian  as 

jt;, .  -  B,  oo.  («*  -  K)  .13^  - ») .  -  /)] 

4  .  -  oos  (<nt  -  ,1)  .  -  >>)]  I  .  (5.3) 


'  U-iq  H_  • 

’fis  -  -  1*13  COB  art  -  -  "T* 


12tDt  -i2u3t 
e  +  e 


Let  us  now  define  the  quantities  which  measure  the  strength  of  the 

perturbations  In  dimensions  of  frequency 


^23^1 


^13  «2 


(5.^) 


and  the  conjugate  quantities 


^21  -  ^*12 

^32  ■  ^23 

6  ■  B* 

31  ^13 


^  'The  ^magnitude  of  nonresonant  perturbations  is  fovind  to  be  of 
order  of  ^  /oi  ,  while  resonant  perturbations  are  of  the  order  of  iH  /idn  , 
where  An  is  the  linewidth.  Since  in  aiiy  reasonable  situation  cu  »  An 
we  may  neglect  the  nonresonant  terms. 
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9 


Vtm  a  earaful  axHdnatloa  of  (9;2),  v«  tee  that  the  "natuzal*' 
freqittBneies:  of  ^e  Tarloua  coqpoaents  of  the  density  oatrlx  an 


t  P22  * 

Bjj  ,  «U  "dc” 

^12  ~  “21 

“ai”  -“ei 

^23  ~  *^32 

^32  ”^32 

P13  ~  «3i 

C! 

1 

a 

Since  only  those  perturbatlws  vhlch  are  almost  resonant  vlU  have  any 
effect,  «e  assume  a  solution  of  the  form 


\  \  ttot 

^11  “  ^11  *  ^12  “  ^12  * 


^22*^22  ^  *‘23*^23*  *  ■  (5*5) 

-  1 

•  A__  ,  P, «  ■  e 

33  33  13  13 


vhere  the  X  are  constants,  which  are  in  general  complex.  By  using 
harmonic  balance,  our  eqxiatlons  become 

^^11  ‘  ^11^  “  ^  ’'l  ^^12  ^21  ^3  ^31  ‘  ^21  ^2  "  ^31  ^3^ 

^^22  ■  -  1  \  Ogi  ^12  *  ^23  ^32  "  ^12  ^21  "  ^32  ^23^ 

^^33  ’  ^33^  -  i  (P31  ^3  +  P32  ^23  ‘  ^3  ^31  ”  ^23  ^32^ 

|(fl2i  "  ®)  +  ^12  •  ‘  ^22^  ^12  *  ^32  ^13  "  ^13  ^32] 

|(«32  -  ®)  +  :;^J  ^23  ■  [^^22  ■  ^33^  ^23  *  ^13  ^21  "  ^21  ^3] 

j(0^  -  2»)  +  ^  j  X^^  -  |(X^^  -  X^j)  +  Pg^  X^  -  p^  Xg^j  . 


(6a) 

(6b) 

(6c) 

(6d) 

(6e) 

(6f) 
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Ve  oBy  now  proceed  to  eolye  these  e4.usti<»s  using  the  foUawlog  assunptlons. 
First,  asBuias  that  the  frequency  <a  is  such  that  2cd  »  0^  ,  and  taJte 
the  fundsaental  fields  to  be  such  larger  thsn  the  hanaonie  fields.  Then 
ve  have  for  the  two  eonvooents  eutid  X^^  their  values  as  detemined 

by  first  order  interactions: 


,  '  ^22^  ^12 
12 


(5.7) 


(Xgg  -  X^^)  ftg. 


(5.8) 


By  substituting  these  in  Eq.  (^.6f)  and  using  the  fact  that  Sen  >  , 

ve  have  the  relation  for  X^^  : 


’‘13  ■  ‘  ’i3  '‘is  <*‘33  ■  *  L  :V\ 

V  ■'12/V  ^23> 


1  1 

1  6(X^^  -  -  ^^22  "  ^3^ - ■  ^22^  »  (5.9) 

‘^12  “^23 


vhere  6  ■  -  ca  ■  -  (n^  -  a*)  .  We  nov  substitute  these  in  the 

eqiuitions  for  the  diagonal  cooiponents,  giving 


e  ^12  ^21  (^11  “  ^22^ 

^11  ^11  ^  ’^l  5  _ 


X  -  X  a  1  T 

33  33  ^  ^1 


^•ao  (^00  "  ^aa) 


-  C.C. 


,  (5.10) 


(^22  "  ^22^  *  "  (^U  "  ^11^  "  (^3  '  ^33^ 
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where  «•  hare  Inelu&ed  cnly  the  larfeat  tame,  and  where  e.c.  is  the 
coevlex  eooj^ugate  ^piantltir.  Bgr  performing  aone  algebra,  we  find  the 
following: 


where 


(A-i-  2B)  X*^)A 

(1  +  2A)(1  +  2S)  -  AB 


(5.11) 


(1  +  A)(l  +  B)  -  (1  +  A)X®  -  (1  +  B)XJ^ 


(1  +  2A)(1  +  2B)  -  AB 


(5.12) 


Xrt«  ■ 

33 


(1  4-  2A)(B  t  X*^)  -  B(A  +  ^) 
(1  +  2A)(1  +  2B)  -  AB 


(5.13) 


hi  “  ^33 


X^  (14-  3B)  -  X^^  (n-  3A)  (A  -  B) 
(1  +  2A)(1  +  2B)  -  AB 


,  (5.14) 


2  ^12  ^21 


6^  + 

12 

L  12J 

(5.15) 


2 


^12  p  -7- 
L  23J 


(5.16) 


C.  MACROBCOFIC  MAGIlEriZ/UriOir  AMD  CAVITT  RBACIKMT 
1.  Magnetization 

Given  the  density  matrix,  we  can  proceed  to  calculate  the  values  of 
all  the  observables  of  the  system.  In  particular,  we  are  Interested  in 
the  eoapcnents  of  the  magnetic  dipole  moment.  The  prescription  for  finding 
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(5.17) 


the  value  of  ib  flhaervahlei  Q  ,  l» 

y'. 

(q)  ■  Tr  (pQ)  . 

We  find  for  the  dipole  eoemt  per  unit  voloae  or  for  the  oagaetlzatlan  at 
the  harmonic 


«  K  (p +  P13  »*3i)  t 


(5.18) 


and  at  the  fundamental, 

®  (Ppi  Pi  9  Pi  9  Ppi  +  P23  P32  ^  ^32  ^23^  ^  (5*19) 

uhere  N  is  the  number  of  systems  per  unit  volume.  The  second  harmonic 
conqponent  of  M  will  determine  the  amount  of  output  power  generated, 
while  the  conponent  at  the  frequency  co  determines  the  amount  of  . 
resistive  loading  at  the  fundamental. 

2.  Cavity  Reaction 

We  may  determine  the  fields  generated  in  the  cavity  by  the  magnetiza¬ 
tion  from  the  relation 


H  «  ^  M  , 
c 


where  the  cavity  susceptibility,  ,  is  given  by 


4]f» 


f 


X  -  0)®  +  1  — 


(5.20) 


(5.21) 


Here,  H  and  M  vary  as  and  we  have  the  definition 

M  f  -  J  JU(t)  .  a^(x)  dV  ,  (5.22) 

oavity 
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vhere  J^Oc)  noma},  node  function  for  the  cavity.  It  Is  eissumed 

that  there  Is  only  one  mode  near  the  frequency  m ,  and  Is  the  loaded 
Q  of  the  cavity  mode.  As  described  In  Chapter  III,  ve  must  allow  for 
different  filling  feetora  for  the  two  cea^ponents  of  ^  substituting 

(^.9)  Into  (?.l8)  and  using  this  In  conjunction  with  (!$.20)  and  (^.21), 
we  have 


where  now  we  have  assumed  that  the  cavity  is  on  resonance,  Is 

defined  as  the  filling  factor  for  the  component  of  proportional  to 


,  fg  Is  the  filling  factor  for  the  component  of  Mg  proportional 
,  and  Qg  Is  the  loaded  Q  of  the  cavity.  By  defining  the 


magnetic  Q  of  the  ssjnqple  at  2u>  ,  ,  as 
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Is  a  saturation  parameter  defined  as 


where 


11 


(5.26) 


D.  CHOICE  OF  6 

Let  us  now  suppose  that  we  are  In  a  position  to  vary  the  value  of 
the  energy  Eg  relative  to  a  fixed  value  of  E^  -  E^  .  In  a  paramag¬ 
netic  solid  this  can  be  done  by  varying  the  angle  and  magnitude  of  the 
magnetic  field  relative  to  the  crystal  axes.  We  shall  look  for  the  value 
or  values  of  6  which  will  maximize  the  nonlinearity  of  the  system.  This 
maximization  will  occur,  when  for  a  fixed  we  maximize  Hg  ,  [Eq.  (5.25)], 

Let  us  first  assume  that  the  levels  are  equally  spaced,  B  =  0  .  In 

this  case  the  term  in  brackets  reduces  to 


(\i  -  >■22)  -  ^ 

23  12 


^^22  ”  ^33^ 


(5.27) 


If  we  now  assume  further  that  the  linewldths  of  the  two  transitions  are 
eq\ial,  Tg^  ss  >  substitute  for  ,  Xgg  ,  and  X^^  from  (5.1l)> 

(5.12),  and  (5.13)f  and  expand  the  Boltzmann  term  for  the  equilibrium 
values,  keeping  only  the  linear  term, we  get  for  (5.27) 


1 

- 2:  (B  -  A) 

Tig2kT 


(5.28) 


We  see  that,  under  these  assumptions,  the  nonlinear  term  is  nonzero  only 
if  B  ^  A  ,  which  Implies  ^  ^*^23^  ’  B  f  A  ,  we  see 

that  the  term  for  B  =  0  depends  on  the  difference  between  quantities 
and  hence  will  not  contribute  significantly  to  the  effect.  We  thus  set 
(5.27)  eqvial  to  zero  and  the  expression  for  Hg  becomes 


H 


2 


^1  ^13 

Se  J 


fg  V* 


ft  ^12  ^23 


11 


^13  Sn2 


(5.29) 
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vhlch  is  maxlndzed  for  &  ■  V‘^12  *  maxlaum 

value  for  Hg  occurs  vhen  Eg  lies  one  llnewldth  from  the  value 
(E^  +  E2)/2  ,  or  alternatively,  when  the  fundamental  alghal  at  o  sees 
both  the  transitions  at  and  0^2  at  their  hiuLf -power  points.  Upon 

setting  6  =  iA]^2  '  ®2 


where  we  have  used  l/Qg  «  l/Q^  +  1/q^  .  From  (5.32)  we  see  that  the 
apparent  source  Q  Is  a  parallel  combination  of  the  cavity  Q  and  a 
modified  magnetic  Q  ,  the  modification  involving  the  filling  factor 
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and  Including  the  effects  of  saturation.  The  source  Q  Is  then 


1 


Q 


source 


(5.33) 


The  question  of  populatlcxi  inversion  is  the  same  as  in  the  case  of  the 

t\<o-level  system:  if  the  population  of  levels  1  and  3  is  in  a  normal 

state,  then  the  quantum  system  presents  los%  reducing  the  source  Q  ; 

if  levels  1  and  3  are  Inverted  then  becomes  negative,  the  source 

Q  (5.33)  is  Increased,  and  the  possibility  for  oscillations  exists. 

Assuming  that  the  external  coiq)ling  is  chosen  to  be  optimum, 

Q  =  Q  ,  our  expression  for  the  power  becomes 

e  source  *  ^ 


P(2a5) 


2  2 
0)  f  ^  V  S 


(^x 


13  ^^12  ^23  12'  "1  ^source 

256,  4  4 


(5.34) 


We  may  consider  the  above  under  two  special  situations, 
low  power  levels  where  there  is  no  saturation,  =  1 
tive  magnetic  Q  ,  f  determines  the  source  Q 
the  case  for  efficient  operation,  we  have  for  the  power 


If  we  are  at 
,  and  the  effec- 
,  as  would  be 


P(2u)) 


0) 


(i* 


12 


256* 


‘13 


Sa2 


f 


while  at  high  power  levels  the  saturation  reduces  the  loading  of  the 
sample  and  we  have  Eq.  (5*3^)  ^ith  *  Qq 


F.  SATURATION 

The  phenomenon  of  saturation  occurs  when  the  populations  of  levels 
1  and  3  tend  to  equalize  and  is  described  by  the  saturation  parameter 
®13  defined  in  Eq.  (5.26).  The  exact  value  of  must  be  found  from 

(5*1^)  and  will,  in  general,  be  complicated.  If,  however,  we  make  the 
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Bliig>llfylng  assumption  that  A  -  B  ,  l.e., 


(5.36) 


When  the  fundamental  field  strength  Is  such  that 


<  > — ^ — r 

“^12  "^l  ^^12 


(5.37) 


then  satxiratlon  will  start  to  occur.  If  we  examine  the  expression  for 
the  output  power  (5.35)  and  ,  (5.24),  we  see  that  the  power  Is  pro- 

portional  to  that  we  desire  large  and  oi"  narrow 

resonance  lines.  Prom  (5.37)  see  that  to  reduce  saturation  we  desire 
as  small  a  value  of  as  possible.  These  two  requirements  point 

out  the  important  fact  that  we  desire  as  small  a  value  of  as  Is 
possible,  consistent  with  a  narrow  linewldth.  This  requirement  is  in 
sharp  contrast  to  maser  theory,  where  saturation  and  a  large  value  of 
are  desired. 

Assuming  that  we  reach  saturatioiv  by  dropping  the  1  in  the 
denominator  of  (5.36)^  find  for  the  saturated  output  power 


01  V  r 

^3  4 


(5.38) 


or  by  substituting  for  ,  Bq.  (5.24), 


(4  -  >‘33)' 


(5.39) 
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G.  DEPENDENCE  ON  PARAMETERS  OF  THE  SYSTEM 


From  this  relation  for  the  saturated  power  and  the  equivalent 
expression  for  the  low-level,  nonsaturated  case  (5"35)>  ''e  can  examine 
the  dependence  on  the  parameters  of  the  atom  and  the  cavity. 


1.  Filling  Factor 

First  of  all  we  see  that  the  output  power  is  proportional  to  the 
volume  of  the  cavity,  V  ,  and  to  the  square  of  the  filling  factor,  fg  , 
defined  in  Eq.  (5.22).  The  dependence  on  the  square  of  the  filling  factor 
emphasizes  the  Importance  of  good  circuit  design. 


2.  Dipole  Moments 

The  low-level  expression  (5«35)  shows  that  the  output  is  proportional 

to  the  square  of  the  dipole  moments  of  the  1-2  and  2-3  transitions  and 

2 

is  Independent  of  the  moment  of  the  1-3  transition  —  the  denom¬ 
inator  Is  cancelled  by  the  same  term  In  .  The  value  of  the  moment 

must  be  large  enough  to  make  *  however.  We  see  that 

we  desire  large  dipole  moments,  especially  at  the  1-2  and  2-3  transitions. 

2 

In  the  saturated  case  the  output  is  proportional  to  and  is  Inde- 


and  p, 


23 


We  thus  see  that  for  high  power  operation, 
For  this  case  u...  and 


1*23  do 


pendent  of 

a  large  is  exxremexy  aeairaDj.e.  For  tnis  case  |i^g 

determine  the  applied  jKJwer  level  necessary  to  cause  saturation  from 
(5.3T)-  The  dependence  on  dipole  moments  is  seen  to  be  critical  in  both 
cases.  If  paramagnetic  substances  were  to  be  used,  this  strong  dependence 
on  p  wovild  point  to  the  use  of  materials  with  large  effective  spins. 


3.  Relaxation  Times  and  Concentration 

In  (5.35)  for  the  unsaturated  case  T^  ,  the  longitudinal  relaxa¬ 
tion  time  does  not  appear.  The  transverse  relaxation  time  appears  in  the 
2 

form  —  the  latter  factor  is  obtained  from  .  The  concen¬ 

tration,  N  ,  comes  in  linearly  through  Q^g  ,  and  so  at  low  power  levels 
a  maximization  of  Nt^  will  result  in  the  maximum  nonlinearity.  At 
saturation  power  levels, the  ratio  t  /t^  points  out  the  desirability 
of  making  T^  and  t  as  nearly  equal  aa  possible  (we  have  the  restriction, 
T^  >  T  ,  from  the  definition  of  these  quantities).  This  less  strict 
condition  on  T^  should  make  possible  the  use  of  a  wider  class  of  mater¬ 
ials  than  maaer  applications. 


Under  saturation  conditions  ve  find  that  the  power  18  proportional 
to  the  square  of  the  quantity  Nt  .  For  high  power  operation  Nt  should 

3 

be  maximized;  for  low  power  levels,  Nt  ,  Since  t  generally 
decreeises  with  concentration,  a  large  change  In  N  usually  results  In  a 
less  significant  change  In  the  product  Nt  ,  although  perhaps  It  results 
in  a  more  important  change  in  Nt-'  .  We  here  neglect  entirely  the 
effects  of  cross -relaxation. 

Temperature  Dependence 

The  last  ina,^or  dependence  Is  that  of  temperatiire.  In  the  formulation 
It  arises  both  in  the  temperature  dependence  of  the  relaxation  times  and  in 
the  Boltzmann  factor.  In  the  latter  case  we  have 


3  ]£T 


> 


assuming  three  levels  only.  The  Boltzmann  factor  comes  in  through  the 
magnetic  Q  ,  yielding  a  lower  value  of  and  hence  a  stronger  non¬ 

linear  effect  for  lower  temperatures.  By  again  assuming  low  power  levels, 
where  (5«35)  applies,  we  see  that  P(2tu)  varies  inversely  an  the  first 

power  of  T  .  For  satxirated  conditions  (5-39)>  the  Boltzmann  factor 

_2 

introduces  a  factor  of  T  but  the  temperature  dependence  of  T.^  will 
in  most  Instances  more  than  cancel  the  T  .  Hence,  low  teiiq>erature 
operation  han  the  effect  of  increasing  the  magnitude  of  the  nonlinearity 
for  low  power  levels  while  increasing  the  effects  of, saturation  for  high 
power  operation.  For  high  power  operation  no  advantage  is  gained  by 
operating  below  room  tenqjerature  unless  a  material  is  used  in  which  either 
T^  is  independent  of  temperature  or  where  t  also,  increases  with 
decreasing  temperature. 

Low  teraperatvure  operation  may,  however,  be  useful  in  another 
application,  namely  mixing.  If  it  were  desired,  for  example,  to  down- 
convert  from  a  frequency  o)^  to  a  lower  frequency  at  low  signal 

levels,  then  a  three-level  system  with  suitably  placed  energy  levels 
could  be  used.  Low  tenqierature  operation  would  Increase  the  magnitude  of 
the  nonlinearity  as  well  as  reduce  circuit  noise.  The  details  of  this 
process  will  not  be  presented  here. 
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5.  Nonresonant  Behavior 

Throughout  this  chapter  we  have  assumed  that  2uj  ■  .  Such  a 

condition  is  not  required  for  hannonic  generation  hut  does  serve  to 
Increase  the  magnitude  of  the  nonlinear  effect  by  reducing  the  resonance 
denominator.  By  operating  on  resonance^  the  nonlinearity  possesses  a 
resistive  character,  while  off  resonance  it  would  be  mainly  reactive. 

Another  aspect  of  off  resonance  behavior  is  concerned  with  the 
choice  of  6  .  Previously  it  was  shown  that  &  »  maximum 

nonlinearity  when  saturation  was  ignored.  If  saturation  is  included,  we 
must  compare  the  variation  in  the  nonlinearity  given  by  (5.29), 


H 


2 


(5.^) 


with  the  more  exact  saturation  expression  for 


1 


(5.^1) 


In  the  limits  of  saturation  we  have 


Hg  «  6 


12 


(5.^) 


which  shows  that  the  saturation  power  should  Increase  with  increasing  & 
because  the  fundamental  absorption  responsible  for  saturation  is  reduced. 

In  the  preceding  discussion  it  has  been  assumed  that  the  fundamental 
field,  ,  is  given  and  is  continuously  variable.  In  actual  practice, 

from  an  applications  point  of  view,  it  will  usually  be  the  case  that  the 
level  of  the  fundamental  power  will  be  fixed.  It  is  then  desirable  to 
maximize  the  output  power  for  a  given  input  power  by  varying  &  . 
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To  do  such  a  maxliBlzatlon  natheaatlcslly  vould  be  an  extremely  Involved 
task  due  to  the  conQillcated  dependence  of  the  efficiency  on  6  . 

Specifically  dependent  on  6  are  the  Q  at  the  fundamental  the 

saturation  (5-36)  which  affects  the  source  Q  (5*33)>,'fcbe  fiudamental 
Q  (5* ^5)  and  the  output  power,  (5*39)  and  finally  the  dependence  not 
specifically  Included  In  the  analysis,  Such  a  maximization  would 

be  most  easily  accomplished  In  practice  by  experimentally  varying  5  , 
for  example,  by  altering  the  direction  and  magnitude  of  an  applied 
magnetic  field. 


H.  FUNDAMENTAL  ABSORPTION  AND  CAVITY  LOADING 

From  Eqs.  (5.19),  (5-7),  and  (5.8),  we  have  the  magnetization  at  the 
fundamental : 


SXgg)  (^33  "  > 

(5.^3) 


where  Is  the  component  of  varying  as  e^^  ,  and  we  have 

again  assigned  “  ^23  ^12  “  ^23  ’  first  term  will,  as 

before,  be  small,  leaving 


2 

^12 


N  H.  T 


1  12  ^^33 
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r2  2  X 

&  T^) 


(5.44) 


By  using  the  cavity  susceptibility,  we  find  that  the  fundamental  loading 
of  the  cavity  by  the  sanple  is  given  by  a  Q 


1 


f  S, 
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V  ^  4) 


(5.45) 
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where 


1 


Snl 


2n  N 

» 


(5. '^6) 


and  f  is  the  filling  factor.  From  this  we  see  that  the  loading  is 
reduced  if  the  filling  factor  is  small  or  if  saturation  has  occurred,  or 
if  the  system  is  detuned,  i.e.,  large  & 


I.  RATIO  OF  FIRST  AND  SECOND  ORDER  MAGNETIZATIONS 

Upon  taking  the  ratio  of  the  magnitudes  of  the  first  and  second 
order  magnetizations  under  the  same  assungjtions  leading  to  the  derivation 
of  the  power  expressions,  we  find 


«1 


H, 


(5.J+7) 


For  a  magnetic  dipole  we  have  pi  ~  10 


~  10*20 


so 


Ml 


If  a  material  had  a  linewidth  corresponding  to  =  0.5  x  lo'^  ,  or 

around  one  megacycle,  then  for  =  1  oe  the  second  order  magnetization 

is  comparable  to  the  first.  For  an  electric  dipole  system  with  its  much 

-l8 

larger  dipole  moments  u  ~  10  »  we  would  have  for  the  ratio  of  the 

polarizations 


^2  8  4 

-*10  T  - 
^1  =1 
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Energy  Levels  (kMc) 
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If  were  10  '  sec,  then  for  ■  0.1  esu  or  30  volts/cm  ,  we 

have  ^2  **  ^1  ■  ^  electric  dipole  system  with  the  desired  selection 

rules  and  energy  levels  would  thus  display  an  enormous  nonlinear  effect. 

J.  EXAMPLE 

As  a  brief  example  let  us  consider  the  use  of  ruby  as  the  active 
substance.  A  plot  of  the  ground  state  S  »  3/2  energy  levels  for 
at  90°  to  the  c-axis  is  shown  In  Fig.  5*2,  where  the  curve  Is  taken  from 

65 

the  analytical  work  of  Chang  and  Slegman.  The  matrix  elements  are 
calculated  to  be 

^12  “  ^23  “ 

^^13  =  , 

20  3 

where  is  the  Bohr  magneton.  We  shall  take  N  =  0.25  10  spins/cm"^; 

-8 

then,  from  Strandberg,  the  linewidth  corresponds  to  t  =  0.4  X  10  sec. 
The  operating  frequency  is  taken  to  be  9*5  kMc  for  the  fundamental  and  19.O 
kMc  for  the  harmonic.  With  the  assumed  spin-lattice  relaxation  times  of 
lO"®  sec  at  300°K,  lO"^  sec  at  TT°K  and  4  X  lo"^  at  4.2°K  and  an 
unloaded  cavity  Q  of  5OOO,  the  output  power  as  a  fiinction  of  fundamental 
field  strength  is  plotted  for  the  three  temperature  ranges  In  Fig.  5- 3* 

The  increase  in  low-level  power  at  low  temgperatures  and  the  effects  of 
satiiratlon  are  clearly  seen  in  this  plot. 

K.  CONCLUSIONS 

We  have  seen  that  the  three-level  system  can  act  as  a  frequency 
converter  and  in  particiilar  as  a  harmonic  generator.  The  dependence  on 
the  parameters  of  the  system  has  been  considered  and  in  particular  the 
desire  for  a  short  spinr lattice  relaxation  time  has  been  pointed  out.  . 

For  microwave  frequencies  the  system  does  not  seem  as  well  suited  for 
high  power  applications  as  does,  for  example,  a  ferrite,  but  at  higher 
frequencies,  using  materials  with  large  spin  and  large  crystal  field 
splittings  or  electric  dipole  transitions,  this  system  may  find  use. 

Perhaps  the  most  useful  application  of  this  system  would  be  els  a  low-level 
mixer  which  was  mentioned  previously. 
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CHAPTER  VI 


EXPERZHEMTAL  EVIDENCE  OF  THIRD  HARMONIC  GENERATION  IN  A 
TWO-LEVEL  SYSTEM 

The  experiment  described  in  this  chapter  was  undertsiken  to  complement 
the  theoretical  work  of  Chapter  III.  The  primary  objective  of  the  experi¬ 
ment  was  to  observe  the  harmonic  generation  effect,  predicted  by  the  theory 
and  previously  unobserved,  and,  if  observable,  to  obtain  as  much  quantitative 
and  qualitative  information  as  possible.  The  experiment  was  successful 
inasmuch  as  the  effect  was  Indeed  observed,  giving  harmonic  powers  in  excess 
of  10  milliwatts.  The  qualitative  results  generally  agreed  with  predictions, 
but  the  quantitative  results,  in  particular  the  output  power,  were  only  in 
fair  agreement  with  theory. 

A.  THE  OPERATING  SUBSTANCE,  NH^ 

The  material  used  in  the  verification  of  the  possibility  of  harmonic 
generation  In  a  two-level  system  was  ammonia  gas.  Specifically,  use  was 
made  of  the  Inversion  transition  occurring  near  kMc.  The  use  of  ammonia 
had  several  advantages;  first,  it  is  readily  available  and  fairly  easily 
handled;  second,  there  is  extensive  knowledge  of  its  spectrum  due  to 
previous  spectroscopic  studies;  third,  its  spectrum  displays  a  very  strong 
absorption  line  in  the  neighborhood  of  2k  kMc  which  is  a  convenient  frequency 
range  with  the  equipment  available. 

It  was  further  decided  to  perform  the  experiment  for  case  I^^^  where 
the  pump  was  applied  near  one-third  of  the  natural  frequency  of  2k  kMc, 
putting  the  fundamental  in  the  range  of  X-band  and  the  third  harmonic  at 
K-band. 

B.  THE  CAVITY 

In  spite  of  the  very  strong  absorption  line  in  ammonia,  the  magnitude 
of  the  predicted  nonlinearity  is  such  as  to  require  a  resonant  system  at 

^^^See  Chapter III. 
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both  the  fundamental  and  the  harmonic;  the  former  Is  required  in  order  to 
get  the  required  fundamental  field  strengths  for  reasonable  power  levels, 
and  the  latter  is  required  since  the  harmonic  power  output  is  proportional 
to  the  cavity  Q  at  the  harmonic.  The  cavity  must  then  be  resonant  at 
o)  and  3n5  ,  preferably  with  high  Q's  at  both  frequencies,  and  must 

have  a  region  in  which  the  electric  fields  at  the  two  frequencies  are 
strong  and  essentially  parallel. 

The  design  decided  upon  was  a  re-entrant  cavity  similar  to  a  fore¬ 
shortened  quarter-wave  line  with  the  interaction  region  near  the  gap.  A 
test  cavity  was  built  in  which  both  the  height  of  the  cavity  and  the  length 
of  the  post  were  variable.  Cold  test  measurements  were  made  at  various 
frequencies  and  plots  of  cavity  height  vs  post  height  were  made  for  a 
fixed  resonant  frequency.  This  was  done  both  for  frequencies  near  the 
fundamental  and  near  the  third  harmonic. 

Finally,  cavity  dimensions  giving  simultaneous  resonances  at  oj  and 
3(u  were  found  by  comparing  the  two  sets  of  iso-frequency  curves.  An 
intersection  satisfying  the  desired  conditions  was  found  as  the  inter¬ 
section  of  the  foreshortened  3X/*+  and  foreshortened  9'^/h  at  the 
fvuidamental  and  harmonic,  respectively.  A  cross  section  of  the  cavity 
is  shown  in  Pig.  6.1  along  with  the  pertinent  dimensions. 

The  coupling  to  the  cavity  at  both  the  fundamental  and  harmonic  was 
achieved  by  the  use  of  magnetic  loops.  These  in  turn  were  fed  from  a 
small  coaxial  line  with  an  OD  of  0.09^  in.  It  was  necessary  to  construct 
a  coax ial-to -waveguide  transition  for  the  harmonic  in  order  to  be  able 
to  work  at  K-band  in  the  waveguide.  The  coupling  loss  at  X-band  was  1.5  db 
and  at  K-band  was  5  db,  giving  a  net  loss  of  greater  than  6  db. 

In  order  to  permit  the  evacuation  of  the  cavity,  "0"  rings  were 
inserted  in  the  joint  between  the  base  and  the  top,  as  shown  in  Pig.  6.1. 

The  major  source  of  vacuum  problems  in  the  cavlLy  resulted  from  leakage 
along  the  center  conductor  of  the  coaxial  coupling  loops .  The  vacuum 
conditions  were  satisfactory  for  short  term  operation,  but  over  long  periods 
of  time  enough  air  would  leak  in  to  cause  breakdown.  This  breakdown  will 
be  discussed  later. 
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FIG.  6.1- -Cross  section  of  the  cavity. 
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C.  EXPERIMEMTAL  SETUP  AND  OPERATION 

A  Bchematic  of  the  experimental  setup  is  shown  in  Fig.  6.2.  The 

power  soxirce  was  a  2J51  magnetron  operated  with  a  2  iisec  pulse  at  50  pps, 

-4 

giving  a  duty  cycle  of  10  .  After  passing  through  a  directional  coupler 

to  lower  the  power  level  the  signal  was  put  through  a  low-pass  filter  in 
order  to  remove  magnetron  harmonics.  It  was  then  passed  through  an 
attenuator  which  was  used  to  vary  the  level  of  the  applied  power.  The 
maximum  available  power  was  about  ■  -eak.  Both  the  indicent  and  reflected 
power  were  monitored  by  crystal  leoectors  and  an  E-H  tuner  was  used  to 
match  into  the  cavity. 

On  the  output  side  (K-band)  the  coaxial  to  waveguide  transition  was 
followed  by  a  slide-screw  t\aner,  a  frequency  meter,  a  precision  attenuator 
and  a  crystal  dete^'tor. 

In  all  the  experiments  performed,  the  gas  in  the  cavity  was  at  a 
constant  pressure  and  was  not  circulated.  In  order  to  determine  if  the 
source  of  the  effect  was  some  general  property  of  the  gas,  dry  nitrogen, 
argon  and  air  were  also  used  in  the  experiment. 

D.  EXPERIMENTAL  RESULTS 

The  experiment  was  performed  using  four  gases;  ammonia,  air,  dry 
nitrogen,  and  argon.  Only  the  ammonia  gave  any  harmonic  generation  under 
the  desired  operating  conditions.  The  cavity  was  matched  to  the  wave¬ 
guide  at  the  fundamental  so  that  the  reflected  power  was  down  10  to  20 
db  from  the  Incident .  Under  these  circumstances  output  powers  on  the 
order  of  10  mw  were  observed  at  the  third  harmonic. 

1.  Breakdown 

Operation  at  high  power  levels  was  frequently  limited  by  breakdown 
of  the  gas  in  the  cavity.  The  input  power  level  at  which  the  cavity 
broke  down  was  a  function  of  the  gas,  the  operating  pressure  and  the 
condition  of  the  cavity,  especially  the  input  loop.  Once  breakdown  had 
occijrred,  the  cavity  had  to  be  disassembled  and  cleaned  before  it  could 
be  used  again.  In  some  runs  maximum  power  (l  kw)  could  be  applied  to  the 
ammonia  without  breakdown,  while  in  others  slight  air  leaks  caused  break¬ 
down  at  low  power  levels.  That  the  breakdown  was  due  to  the  air  could  be 
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C.  EXPERIMEMTAL  SETUP  AND  OPERATION 

A  schematic  of  the  experimental  setup  is  shown  in  Fig.  6.2.  The 

power  source  was  a  2J51  magnetron  operated  with  a  2  ^lsec  pulse  at  50  pps, 

-4 

giving  a  duty  cycle  of  10  .  After  passing  through  a  directional  coupler 

to  lower  the  power  level  the  signal  was  put  through  a  low-pass  filter  in 
order  to  remove  magnetron  harmonics.  It  was  then  passed  through  an 
atteniiator  which  was  used  to  vary  the  level  of  the  applied  power.  The 
maximum  available  power  was  about  1  kw  peak.  Both  the  indicent  and  reflected 
power  were  monitored  by  crystal  detectors  and  an  E-H  tuner  was  used  to 
match  into  the  cavity. 

On  the  output  side  (K-band)  the  coaxial  to  waveguide  transition  was 
followed  by  a  slide-screw  tuner,  a  frequency  meter,  a  precision  attenuator 
and  a  crystal  detector. 

In  all  the  experiments  performed,  the  gas  in  the  cavity  was  at  a 
constant  pressure  and  was  not  circulated.  In  order  to  determine  if  the 
source  of  the  effect  was  some  general  property  of  the  gas,  dry  nitrogen, 
argon  and  air  were  also  used  in  the  experiment. 

D.  EXPERIMENTAL  RESULTS 

The  experiment  was  performed  using  four  gases;  ammonia,  air,  di'y 
nitrogen,  and  argon.  Only  the  ammonia  gave  any  harmonic  generation  under 
the  desired  operating  conditions-  The  cavity  was  matched  to  the  wave¬ 
guide  at  the  fundamental  so  that  the  reflected  power  was  down  10  to  20 
db  from  the  incident.  Under  these  circumstances  output  powers  on  the 
order  of  10  mw  were  observed  at  the  third  harmonic. 

1.  Breakdown 

Operation  at  high  power  levels  was  frequently  limited  by  breakdown 
of  the  gas  in  the  cavity-  The  input  power  level  at  which  the  cavity 
broke  down  was  a  function  of  the  gas,  the  operating  pressure  and  the 
condition  of  the  cavity,  especially  the  input  loop.  Once  breakdown  had 
occurred,  the  cavity  had  to  be  disassembled  end  cleaned  before  it  could 
be  used  again.  In  seme  runs  maximum  power  (l  kw)  could  be  applied  to  the 
ammonia  without  breakdown,  while  in  others  slight  air  leaks  caused  break¬ 
down  at  low  power  levels.  That  the  breakdown  was  due  to  the  air  could  be 
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FIG.  6.h — Experimental  resiiltB  for  high  pressure  I’cmge. 
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E.  DISCUSSION 

1.  Qualitative 

In  general,  most  all  of  the  experimental  results  qualitatively  agree 
with  the  theoretical  predictions.  First  of  all,  for  a  fixed  low  input 
power  level  the  output  power  increases  with  pressure  up  to  around  100  mm  Hg, 
after  which  it  Is  relatively  independent  of  pressure.  This  is  in  agree¬ 
ment  with  the  fact  that  the  overall  intensity  of  the  ammonia  line  increases 
with  pressure  due  to  the  overlap  of  the  broadened  individual  ccanponents 
as  explained  in  Chapter  III,  section  H.4.  This  increase  should  continue 
until  the  width  of  an  individual  component  approximately  equals  the  overall 
frequency  spread  of  these  components.  At  the  pressure  of  100  mm  Hg,  the 
full  linewldth  is  6  kMc,  which  incliides  most  of  the  strong  ccmponents  of 
the  line,  thus  giving  a  good  agreement. 

Since  the  operation  corresponds  to  3aj  «  0  ,  saturation  should  be 

caused  by  the  qiaantum  system  absorbing  the  generated  harmonic  power.  In 
this  case  the  saturation  denominator  is  of  the  form  (l  +  const  T^T2E^) 
where  =  T^  and  is  inversely  proportional  to  pressure  and  is  the 

fundamental  field  strength.  The  satiiration  should  then  be  minimized 
(saturation  denominator  minimum)  for  higher  pressiues.  A  comparison  of 
the  experimental  curves  shows  this  to  be  true,  the  effects  of  saturation 
being  minimum  at  the  highest  pressure  examined,  300  mm  Hg. 

2.  Quantitative 

A  qxiantative  comparison  of  the  experimental  results  with  theory,  in 

particular  the  magnitude  of  the  harmonic  power,  is  quite  difficult  because 

of  the  large  number  of  iinknown  factors  in  the  experimental  setup.  The 

most  critical  unknown  is  the  filling  factor,  fg  ,  which  depends  upon 

the  magnitude  and  orientation  of  the  fundamental  and  harmonic  fields  within 

the  cavity.  Because  of  the  use  of  higher  order  modes  and  the  irregular 

geometry  of  the  cavity,  the  filling  factor  is  not  known  even  approximately. 

This  is  further  canpllcated  by  the  fact  that  the  assumption  that  the 

interaction  occurs  only  in  the  gap  region  appears  unjustified .  Since  the 

2 

output  power  depends  on  fg  ,  this  uncertainty  is  magnified.  We  shall 
then  only  be  able  to  speak  in  order-of-magnit\ide  terms. 
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CeJ-culatlng  the  field  strength  In  the  gap  region  from  a  knowledge 
of  R/Q  ;  Q  ,  and  the  dimensions  of  the  gap  gives  a  value  equal  to 
200  esu  at  an  Input  power  level  of  1000  watts.  Assuming  the  Interaction 
to  occur  In  the  gap,  as  originally  anticipated,  the  saturation  should  be 
much  larger  than  observed .  Further,  the  power  expected  f rcan  f lelds  of  200 
esu  is  much  larger  than  observed.  Evaluating  the  fields  in  the  other  parts 
of  the  cavity  from  a  uniform  field  approximation  yields  numbers  much  more 
in  line  with  the  observed  output  power  levels  and  saturation.  This  then 
points  to  the  distinct  possibility  that  the  power  Is  generated  throughout 
the  cavity  and  not  Just  in  the  interaction  region.  From  the  results 
obtained  by  use  of  this  cavity  it  is  impossible  to  decide  exactly  where 
the  interaction  takes  place. 

By  assimlng  some  order  of  magnitude  values  for  the  quantities  involved, 

we  may  make  a  comparison  with  theory.  Assuming  E,  =  50  esu  for 

3  ^ 

=  1000  watts,  V  =  6  cm'^  ,  =  0.1  ,  f =  1  ,  the  power  expected  from 

the  cavity  is  (from  Fig.  3.2,  Chapter  III)  P(3aj)  =  250  mw  ,  compared 
to  an  observed  value  of  10  mw.  Including  the  coupling  losses  of  over 
6  db,  the  disagreement  is  near  6  db,  which  is  well  within  the  limits  of 
the  above  assumptions. 

F.  CONCUJSIONS 

From  this  work  we  can  conclude  that  all  evidence  points  to  the 
existence  of  the  nonlinear  effect  predicted.  The  results  do  not,  however, 
allow  a  detailed  comparison  with  theory  because  of  the  many  experimental 
uncertainties  involved.  In  order  to  more  fully  understand  the  effect  and 
to  compare  experiment  with  theory,  a  much  more  detailed  experiment  should 
be  performed.  For  such  an  experiment  the  cavity  should  be  of  a  regular 
geometrical  shape  so  that  the  field  configurations  are  known  analytically. 
In  this  way  both  the  filling  factors  and  the  magnitude  of  the  fields  may 
be  found.  One  possible  cavity  which  might  be  used  is  one  in  which  both 
modes  are  TEM,  such  as  the  shorted  coaxial  ].ine.  For  this  case  the  filling 
factor  may  approach  0.25.  In  a  practical  application,  the  Fabry-Perot 
resonator  would  appear  as  the  best  cavity,  yielding  a  high  Q  and  large 
volume . 
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CHAPTER  VII 


EXPERIMENTAL  EVIDENCE  OF  SECOND  HARMONIC  GENERATION 
IN  A  THREE-LEVEL  PARAMAGNETIC  SYSTEM 

In  Chapter  V,  a  theory  vas  derived  which  suggested  that  it  should 
be  possible  to  obtain  a  second-order  nonlinearity  in  a  three-level  system. 
Such  a  nonlinearity  would  allow  the  mixing  of  two  frequencies  and 

cug  to  obtain  their  sum  and  difference  frequencies,  a  special  case  of  which 
is  second  harmonic  generation.  The  detailed  aspects  of  this  process  were 
considered  in  Chapter  V.  Briefly,  when  the  three  levels  aore  approximately 
equally  spaced  and  the  selection  rules  allow  transitions  between  all 
three  levels,  then  second  harmonic  generation  may  occur  when  the  output 
frequency,  2cu  ,  is  equal  to  the  frequency  spacing  of  the  extreme  levels, 
(E^  -  •  By  using  a  paramagnetic  ion  as  the  active  substance  and 

by  varying  the  magnitude  of  the  applied  dc  magnetic  field  as  well  as  its 
direction  relative  to  the  crystal  axes,  the  above  reqviirements  can  be 
satisfied.  Since  this  effect  had  not  previously  been  experimentally 
observed,  a  preliminary  experiment  designed  as  an  attempt  to  observe  the 
effect  was  performed.  This  experiment  and  the  resiilts  will  now  be 
described. 

A.  MATERIAL;  RUBY 

The  material  used  in  the  experiment  was  ruby.  This  was  chosen 
because  of  its  availability  and  the  fact  that  its  spectrum  is  well  known. 

65 

Chang  and  Slegman  have  calculated,  from  the  known  spin  HamJ.ltonian, 

the  energy  levels  and  matrix  elements,  and  the  latter  have  been  experi- 

66 

mentally  checked  by  Ammann.  These  calculations  of  Chang  and  Slegman 
were  used  to  choose  the  approximate  operating  point  and  no  preliminary 
spectroscopy  was  done  to  check  their  calculations.  The  ruby  sample  used 
was  borrowed  from  R.  Morris,  and  was  a  cylinder  of  pink  ruby  10  mm  in 
diameter  and  2.5  mm  hl^.  The  c-axis  of  the  crystal  was  subsequently 
found  by  X-ray  diffraction  techniques  to  lie  2*5°  from  the  normal  to  the 
crystal  face. 
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FIG.  7-1 — Ruby  sample  used. 


B.  OPERATING  POINT 

Because  of  the  limited  availability  of  rf  pover  sources  the 
operating  frequency  was  chosen  to  be  9-5  fundamental,  and  I9.O  kMc 
harmonic.  From  Chang  and  Siegman,  an  operating  point  in  ruby  satisfying 
the  conditions  mentioned  above  and  corresponding  to  the  frequency  9*5  kMc 
should  occur  with  the  direction  of  the  dc  magnetic  field  at  90°  to  the 
optical  or  c-axls  with  magnetic  fields  of  the  order  of  2000  oersteds. 
Denoting  the  ruby  c-ajcls  as  z  and  the  plane  of  Hq  as  the  x-z  plane, 
then  Hq  Is  along  the  x-axls.  The  strongest  matrix  elements  for  the 
fundamental  Interaction,  corresponding  to  the  1-2  and  2-3  transitions, 
are  found  to  be  along  the  y-axls,  while  those  for  the  harmonic  or  1-3 
transition  are  along  the  dc  field  or  x-axls,  the  same  relative  orientation 
as  in  the  case  of  second  harmonic  generation  in  a  ferrite. 

C.  THE  CAVITY 

Because  of  the  anticipated  low  level  of  the  observed  heurmonic,  it 
was  decided  to  perform  chc  experiment  in  a  microwave  cavity  resonant  at 
both  fundamental  and  harmonic  in  order  to  enhance  the  magnitude  of  the 
effect.  It  was  further  decided  to  place  the  sample  against  the  wall  of 
the  cavity  where  the  rf  magnetic  fields  parallel  to  the  wall  woxild  be 
large.  Because  of  the  spatial  orthogonality  of  the  fundamental  field 
and  the  resultant  harmonic  magnetic  dipole  moment  Ii^(2aj)  ,  it  was 
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necessary  to  choose  cavity  modes  vhlch  vould  have  orthogonal  H  fields, 
both  of  which  were  parallel  to  the  wall  where  the  sangple  was  to  be  placed. 

A  piece  of  rectangular,  X-band  wavegiilde  shorted  at  one  end  and  aper¬ 
ture  coupled  at  both  the  fundamental  and  second  harmonic  was  used  as  the 
cavity  (see  Fig.  1.2).  At  the  fvmdamental,  the  cavity  was  resonant  in 
the  mode,  while  at  the  harmonic  several  closely-spaced  modes  were 

available  for  use.  Looking  from  the  X-band  aperture,  the  and 

TM^2^  were  among  the  suitable  modes.  In  cold  test  without  the  ruby 
sample,  the  frequencies  of  these  modes  were  found  to  lie  where  they  were 
expected.  With  the  ruby  sample  in  the  cavity,  however,  the  modes  near 
the  second  harmonic  were  perturbed  ccxislderably  and  the  precise  deter¬ 
mination  of  the  mode  conflguratioiie  was  not  possible.  Since  the  sample 
did  not  affect  the  X-band  resonant  frequency  significantly  (lying  in  the 
weak  electric  field  region)  it  was  easy  to  arrange  the  dimensions  for 
fundamental  resonance.  In  order  to  obtain  the  K-band  resonance  a  movable 
short  was  attached  to  the  cavity  and  the  sample  was  placed  on  the  wall 
near  the  short.  Due  to  the  large  dielectric  constant  of  the  ruby  (e  lO) 
the  K-band  fields  were  drawn  into  this  region  and  the  tuner  was  made  more 
effective.  (Since  this  waveguide  was  cut  off  to  the  fundamental,  it  did 
not  affect  the  X-band  resonance.)  By  this  means,  a  2:1  ratio  of  the 
resonant  frequencies  was  achieved.  The  loaded  Q's  at  fundamental  and 
harmonic  were  both  300- 

D.  EXPERIMENTAL  SETUP 

A  schematic  drawing  of  the  experimental  setup  is  shown  below  in 
Fig*  7*3*  Tbe  power  soiirce  was  a  2J51  magnetron  operated  in  a  pulsed 
mode  with  a  pulse  repetition  rate  of  60  pps  and  a  pulse  length  of  O.38  usee. 
A  shorting  switch  was  used  to  measure  both  the  incident  power  and  the 
power  reflected  from  the  cavity  xmder  operating  conditions.  A  low  pass 
filter  with  a  rated  40  db  rejection  was  used  to  reduce  the  harmonics 
generated  from  within  the  magnetron.  A  Varian  magnet  was  used  to  provide 
the  necessary  magnetic  fields  and  a  sinQJle  biased  crystal  detector  was 
used  for  detection  of  the  harmonic  output  power.  Urider  operating  con¬ 
ditions  the  cavity  was  matched  in  at  the  fundamental  so  that  the  reflected 
power  was  12  db  below  the  incident. 
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FIG.  7.2— Cutaway  drawing  of  the  cavity. 
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FIG.  7-3 — Schematic  of  experimental  setup. 


E.  RESULTS 


The  following  are  the  qvialitatlve  results  of  the  experiment: 

(a)  Under  the  application  of  several  kilowatts  of  X-band  power  at 
9.503  kMQ  heirmonic  powers  on  the  order  of  tens  of  microwatts  were  observed 
at  19.006  kMc. 

(b)  For  a  fixed  level  of  Input  power  the  magnitude  of  the  output 
power  was  critically  dependent  on  the  magnitude  of  the  applied  magnetic 
field  and  its  angle  relative  to  the  crystal  axis.  No  power  was  observed, 
to  within  the  level  of  the  detection  scheme,  when  the  magnet  was  off 
resonance.  The  width  of  these  resonances  was  typically  20-30  oe. 

(c)  As  the  angle  Q  of  the  magnet  relative  to  the  crystal  axis 
was  varied,  the  magnetic  field  at  which  optimum  interaction  occurred 
changed  in  qualitative  agreement  with  a  curve  of 

E^  (0  ,  Hq)  -  E^  (e  ,  Hq)  =  const. 

(d)  For  a  fixed  angle  0  there  appeared  three  values  of  magnetic 
field  at  which  harmonic  power  occurred,  two  near  I9OO  oe  and  one  near 
1700  oe.  The  latter  response  was  small  comtpared  to  the  other  two. 

F.  INTERPRETATION 

First  of  all,  from  the  fact  that  the  harmonic  power  observed  was 
critically  dependent  on  the  magnitude  of  the  magnetic  field  and  its 
orientation,  it  is  concluded  that  a  nonlinear  process  within  the  material 
itself  is  responsible  for  the  harmonic  generation.  From  the  relatively 
good  agreement  between  the  experimental  operating  conditions  (particularly 
the  value  of  H^)  and  those  predicted, it  is  felt  that  the  nonlinear 
process  described  in  Chapter  V  is  responsible. 

The  reasons  for  the  relatively  low  efficiency  are  felt  to  be  the 
following:  First  and  foremost  is  the  fact  that  the  filling  factor,  f  , 
of  the  ruby  is  most  probably  very  small.  This  arises  from  the  fact  that 
the  ruby  sample  was  placed  in  the  cavity  so  as  to  allow  for  enough  tuning 
at  K-band  to  satisfy  the  2:1  resonant  frequency  requirements,  and  not 
placed  BO  as  to  atteD5)t  to  maximize  the  filling  factor.  Since  the  output 
power  depends  on  the  sqxiare  of  the  filling  factor,  a  small  value  of  that 
quantity  woxild  severely  reduce  the  power  out.  From  the  saii®le  dimensions 
and  the  cavity  size  the  filling  factor  is  at  most  O.O3  euid  probably  near 
0.01. 
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The  second  reeuson  for  lower  efficiency  has  to  do  with  saturation. 
Althoxigh  the  sample  may  not  be  oriented  correctly  relative  to  the  funda¬ 
mental  and  harmonic  fields  for  optimum  harmonic  generation  (small  f  ), 
it  will  still  absorb  the  fundamental  pcuer  and  will  tend  to  saturate.' 

It  is  felt  that  saturation  had  started  to  occur  prior  to  achieving  detec¬ 
table  harmonic  power  levels  at  input  power  levels  «  200  watts.  In  a 
more  detailed  experiment  it  would  be  necessary  to  know  the  cavity  fields, 
in  the  presence  of  the  ruby  sample,  quite  precisely  in  order  to  optimize 
the  filling  factor. 

it  shoiild  be  noted  that  the  requirements  on  the  orientation  of  the 
cavity  fields  are  more  strict  here  than  in  the  case  of  the  maser.  This 
is  a  result  of  the  fact  that  in  a  maser  the  pump  signal  saturates  its 
transition  followed  by  the  amplifying  effect  at  the  signal  transition. 

As  such,  both  frequencies  need  only  independently  see  strong  transitions 
in  the  crystal,  typically  in  the  x-y  plane.  In  the  case  of  harmonic 
generation,  the  re-radiation  of  2to  photons  occurs  d|lmultaneously  with 
the  absorption  of  two  photons  at  to  ,  and  hence  the  relative  orientations 
in  space  of  the  fundamental  and  harmonic  cavity  fields  is  critical. 

Since  the  nonlinearity  and  hence  the  magnitude  of  the  harmonic 
generated  depends  strongly  on  the  location  of  the  middle  level  relative 
to  the  midpoint  between  the  outer  levels,  and  this  in  turn,  for  a  fixed 
00  ,  is  dependent  on  the  operating  point  (6  ,  H^).,  the  power  out  will  be 
a  critical  function  of  the  operating  point.  Unfortunately  the  location 
of  the  c-axls  was  not  known  to  within  3°  so  no  precise  knowledge  of  the 
operating  point  was  available. 

Taking  into  account  the  above  factors,  it  is  felt  that  a  careful 
optimization  would  result  in  several  orders  of  magnitude  Improvement  in 
the  efficiency. 

We  may  explain  the  appearance  of  three  operating  points  where 
harmonic  generation  occurred  from  the  following  argument:  If  the  expres¬ 
sion  for  the  harmonic  power  were  written  in  general  terms  where  we  do 
not  assume  that  2u>  =  (E^  -  E^)/b  ,  then  the  expressions  for  output 

power  would  contain  a  general  denominator  of  the  form 

{2oi  -  (cn  -  ^2i)  ■  ^32^  '  presents  the  possibility  of  three 

points  of  resonant  behavior,  namely  when  2cx>  =  ,  cn  ■  , 

0)  =  ,  for  fixed  cu  .  The  dependence  on  angle  and  magnetic  field 
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of  these  three  operating  points  Is  only  approximately  known  and  would 
have  to  be  Investigated  In  detail  experimentally.  Qualitatively,  however, 
for  cn/2x  ■  and  0  **  90°  ,  the  resonance  conditions  2m  >  0^ 

and  (u  -  should  be  satisfied  for  magnetic  fields  differing  only 
subtly  In  value,  while  oj  ■  will  be  satisfied  for  a  much  lower 
value  of  Hq  •  These  predictions  are  a  consequence  of  an  analysis  of 
the  known  data  on  ruby.  The  experimental  data  are  again  qualitatively  In 
agreement  with  these  predictions.  Here,  as  In  the  determination  of  the 
magnitude  of  the  power  output,  a  precise  knowledge  of  E^(d  ,  H)  Is  needed 
for  any  measure  of  the  quantitative  agreement. 

G.  CONCLUSIONS 

From  the  results  of  this  experiment  we  can  conclude  that  the  predicted 
nonlinear  effect  does,  in  fact,  exist,  but  that  Its  magnitude  appears 
to  be  smaller  than  anticipated.  An  exact  quantitative  ccsooparison  was  not 
possible  due  to  the  large  uncertainty  in  some  of  the  parameters  involved. 
Such  a  quantitative  experiment  should  be  performed  both  to  better  check 
the  theory  and  to  determine  the  feasibility  of  possible  applications.  In 
such  an  experiment  such  things  as  temperature  and  concentration  dependence 
as  well  as  the  effect  of  the  location  of  the  center  level  should  be 
examined.  In  particular,  the  prediction  that  optimum  harmonic  generation 
occurs  when  the  levels  are  not  eqtwily  spaced  shovild  prove  an  interesting 
check  of  the  theory. 

Other  experimental  work  along  this  line  could  involve  the  use  of 
different  materials  with  larger  zero-field  splittings  for  use  at  higher 
frequencies  as  well  as  the  examination  of  materials  with  faster  spin- 
lattice  relaxation  times  to  reduce  the  effects  of  saturation.  Considera¬ 
tion  of  processes  other  than  harmtxaic  generation,  such  aa  mixing,  should 
also  prove  Interesting. 
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CHAPTER  Vni 


CONCLUSIONS  AND  RECCIQ4ENDATI0NS 

The  purpoae  of  this  reseai'ch  has  been  to  study  the  phenomena  of 
nonlinear  or  multiple  quantm  effects  and  to  discuss  possible  applications. 

We  have  found  that  a  q\aantum  mechanical  system  with  discrete  energy  levels 
can  behave  in  a  nonlinear  as  well  as  in  a  linear  manner.  The  nonlinear 
effects  occur  as  a  result  of  iateractlons  which  Involve  more  than  a  single 
quantum  of  electromagnetic  radiation,  whereas  the  linear  effect,  used  in 
masers,  is  a  single  qiiantum  effect.  These  multiple  qiiemtum  phenomena 
then  make  it  possible  to  Use  a  suitable  quantum  system  for  nonlinear 
applications . 

A  general  discussion  of  multiple  quantum  phenomena  and  the  mechanism 
of  indirect  coupling  responsible  for  the  effect  was  presented  in  Chapter  II. 
Following  this  discussion,  three  specific  applications  were  considered: 
harmonic  generation  in  a  two-level  systm;  parametric  effects  in  a  two- 
level  system;  harmonic  generation  in  a  three-level  system.  Although  these 
are  admittedly  simple  quantum  systems,  being  at  best  an  approximation  to 
the  real  case,  their  solutions  do  yield  insight  into  more  complicated 
systems .  Further,  in  many  cases  only  a  few  levels  will  really  be  of 
Importance  in  an  interaction  and  the  problem  can  be  simplified  to  a  two-, 
three-  or  perhaps  a  four-level  problem  soluble  by  the  methods  presented 
here. 

The  general  character  of  these  phenomena  differ  from  single  quantum, 
masqr-like  phenomena.  In  particular  it  is  fo\ind  thiat  poptLLation  inversion 
is  not  required  for  many  of  the  phenomena  and  that  the  saturation  effects 
are  generally  to  be  avoided.  These  requirements,  along  with  the  strong 
dependence  on  the  magnitude  of  the  dipole  moments  involved,  make  the 
criteria  for  the  choice  of  materled-s  different  than  for  maser  applications. 
This  may  make  possible  the  lise  of  new  materials  and  will  require  considerable 
materials  reseeurch. 
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The  ultimate  application  of  these  principles  vlU  probably  be  In  the 
millimeter  wave  to  optical  region  of  the  spectrum^  where  at  present  suitable 
nonlinear  materials  are  not  known.  Using  the  nonlinear  effects  discussed 
here,  mixing,  modulation,  demodulation  and  other  applications  should  be 
X>osslble  at  these  frequencies. 

Because  of  the  fact  that  the  study  of  this  field  is  relatively  new, 
little  has  been  done.  There  then  exists  much  work  to  do  in  the  future. 

In  particular,  more  experimentation  Is  needed  in  order  to  verify  the 
qxiantltatlve  aspects  of  the  theory.  Should  these  experi^mental  results 
appear  as  promising  as  the  theory  indicates,  then  memy  and  various 
applications  can  be  considered.  When  this  state  is  reached  an  intensive 
study  of  materials  will  have  to  be  made  to  find  suitable  ones  for  the 
various  desired  applications. 

In  summary,  this  new  field  appears  quite  interesting  scientifically, 
both  from  the  point  of  view  of  providing  a  better  understanding  of  radiation 
processes  and  frcxn  the  application  standpoint.  Much  research  is  needed, 
however,  before  the  ultimate  feasibility  of  application  is  known. 
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APPENDIX  A 


THE  DENSITY  MATRIX 

It  is  the  purpose  of  this  appendix  to  familiarize  the  reader  with  some 
aspects  of  the  density  matrix,  in  particular  its  description  of  the  inter¬ 
action  of  radiation  with  matter.  No  attempt  will  be  made  to  prove  the 
relations  stated;  for  those  interested  in  proofs  and  further  details, 
references  67-71  nay  be  consulted. 

Let  us  suppose  we  have  a  series  of  possible  states  of  the  atom 
which  we  may  expand  in  terms  of  a  complete  set  of  kets  In) 


^i  =  X  ^n 


.  (A.l) 


The  expectation  value  of  an  operator  Q  in  the  state  is  defined  as 


<Q>i  = 


.  (A.2) 


In  terms  of  the  kets  |n)  this  may  be  written 


n,m 


=  I 


where 


®in  ®^im  ^nm 


n,m 


,  (A.3) 


nm 


.  (a.M 


=  <n|Q|m> 


If  the  atom  may  be  in  one  of  a  group  of  states  with  a  probability 
CD^  ,  where  we  have  the  obvious  restriction 


I 


(A. 5) 


then  the  expectation  value  of  the  operator  Q  in  the  arbitrary  state  i|r 
is 

<\|f|Q|i>  =  ^  0)^  <Q>^  ,  (A. 6) 

i 


where  we  have  taken  into  account  the  statistical  behavior  of  matter  as 
well  as  the  intrinsic  statistical  nature  of  quantum  mechanics,  Eq.  (A.3)- 
By  substituting  Eq.  (A. 3)  in  Eq.  (A. 6),  we  obtain 

<«  =  y.'L  “i  “L  “i»  V  • 

i  n,m 

Define  the  density  matrix  p  as 

p  =  )  CD.  a*  a.  ,  (A. 8) 

mn  i  in  Im  »  \  / 


and  then 


n,m 


By  performing  the  sum  over  n  ,  we  have 

<«>  =  X 

m 

=  Tr  (pQ) 


(A. 9) 


(A. 10) 


(A. 11) 
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vblch  gives  the  prescription  for  finding  the  expectation  value  of  an 
observable. 

Ve  shall  state,  without  proof,  two  addltlcnal  properties  of  the 
density  matrix  which  ve  shall  use: 

(a)  Tr  p  =  1  , 

or 


I 


(A. 12) 


(b) 


(A.13) 


The  first  relation  states  that  the  sun  of  the  diagonal  elements  is  unity. 

If  we  should  use  the  energy  representation  where  our  kets  are  the  elgenkets 
of  ,  the  uni>ert\irbed  Hamiltonian,  then  the  n\inbers  may  be 

Interpreted  as  occui>atlon  nuniberB  of  the  various  states,  n  ,  normalized 
to  a  single  molecvile.  In  this  representation  it  Is  easily  seen  that  each 
diagonal  element  Is  non>negatlve;  this  Is  true  In  general,  however,  for 
any  representation. 

The  second  relation  merely  states  that  the  density  matrix  is  Hermitian 
This  simplifies  solving  for  the  density  matrix  as  it  reduces  by  half  the 
number  of  off-diagonal  elements  to  be  found. 

We  should  like  to  now  find  the  time  development  of  the  density  matrix. 
Let  us  divide  the  Hamiltonian  into  three  parts: 

J{  =  +  V  .  (A.ll^) 

The  first  term  Is  the  Hamiltonian  of  the  atom,  including  the  possibility^ 
of  static  magnetic  fields,  and  Is  time-independent.  The  second  term,  ^ 
shall  represent  the  perturbation  due  to  electromagnetic  radiation.  The 
third  term,  V  ,  represents  the  interaction  between  the  atoms  and  the 
lattice  and  between  the  atoms  themselves  and  will  give  rise  to  spin-lattice 
atkl  spin-spin  relaxation. 
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By  explicitly  writing  down  the  components  of  the  density  matrix, 
taking  their  time  derivatives,  and  using  SchrSdlnger's  equation,  the 
equation  of  motion  for  p  can  be  shown  to  be 

Ihp  =  [3i,p]  ,  (A. 15) 

where  [Jt,p]  Is  the  commutator  of  these  quantities.  We  shall  now  assume 
we  are  in  the  energy  representation, 

(Jl*)  =  E  B  ,  (A.i6) 

and  shall  for  the  time  being  neglect  the  interaction  V  in  the  Hamiltonian. 
The  equations  then  become 

n  =  m  : 

^*^nn  ^  ^nk  ^kn  ’  '^nk^kn^  '  (A.l?) 

k 

and 

n  ^  m  : 

Ihp  =  (E  -  E  )p  +  )  (Ji  V  Pw  -  P  ,  •  (A.l8) 

nm  n  m'  nm  '  nk  km  nk  km'  ' 

k 

These  equations  suffer  from  the  fact  that  they  do  not  give  the  correct 
equilibrium  value  for  p  in  the  limit  Jt  -♦0  .  It  is  shown  in  texts 
on  the  siibject  that  in  thermal  equilibrium  the  density  operator  is  given 

ty 


exp  (- 

p  =  - .  (a.19) 

Tr  expC-Jt^/kT) 


From  this  expression  we  see  that  the  off-diagonal  elements  are  zero 
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while  the  diagonal  elements  are  given  by 


P 


nn 


exp  (-  E^/kt) 

^  exp  (-  E^/kT) 
m 


(A. 20) 


From  Eqs.  (A. 17)  and  (A.l8)  it  is  seen  that  if  the  system  is  perturbed 

from  equilibrivm  and  the  perturbation  is  removed,  then  will  remain 

constant  while  p  will  oscillate  sinusoidally  with  a  frequency 
nm 

01  =  (E  -  E^)/fi  which  is  not  in  agreement  with  (A. 19).  The  reason  for 
this  is  that  the  relaxation  mechanisms  which  give  rise  to  thermal 
equilibrium  have  been  neglected.  As  the  exact  nature  of  the  perturbation 
V  is  not  known,  its  effects  are  usually  phenomenologically  added  via 
longitudinal  and  transverse  relaxation  times,  T^^j  and  ,  respectively. 

With  the  transverse  relaxation  times  added,  the  off-diagonal  equations 
become 


m  : 


ih 


i«p_  +  (E„  -  E^)p^  +~P 
nm  m  n  nm  t 


nm 


nm 


I 


p  -  P  ^  ^w,)  f  (A. 21) 
nk  '^km  nk  km''  *  '  * 


where  to  satisfy  hermiticity  we  require  .  By  the  insertion  of 

this  relaxation  tenn,  p  is  seen  to  decay  to  zero  when  3i  =0  ,  as 

nm 

required  by  (A. 19). 

For  the  diagonal  equations  we  introduce  the  trsnsition  probability, 

,  defined  as  the  probability  per  unit  time  that  the  system  will  make 
a  transition  fran  state  i  to  s'tate  J  due  to  relaxation  processes .  The 
diagonal  equations  then  becone 


^  »  I  “to  ■  »■«.  “nk) 

k 


X  ^^nk  Pkn  "  Pnk^kn^ 


(A. 22) 
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In  words,  these  two  forms  of  relaxation  times  have  the  following  Interpre¬ 
tation: 

(a)  The  transverse  relaxation  time  x  gives  the  decay  time  of  the 

^  J 

term  which  Is  a  measure  of  the  correlation  between  states  1  and  j 

and  is  the  usual  Tg  of  paramagnetic  resonance  terminology.  It  Is  propor¬ 
tional  to  the  inverse  of  the  linewldth  of  the  transition  i  -»  j 

(b)  The  longitudinal  relaxation  time  T . ,  is  a  measure  of  the  time 
required  for  the  transition  1  J  to  thermalize  with  the  lattice  and  is 
the  T^  of  paramagnetic  resonance. 

We  may  make  a  simplification  of  Eq.  (A. 22)  if  we  assume  all  longitudinal 
relaxation  times  eq\ial  to  T^  .  In  this  case  we  have 

Ifi^  +  ^  (p  -  P®  )  =  y  V  P,,  -  P  u  K  )  •  (A. 26) 

T,  '  nn  nn'  ^ '  nk  ^kn  nk  kn'  '  ' 

^  k 

Equations  (A. 21)  and  (A. 22)  constitute  the  density  matrix  formulation 
of  the  interaction  of  radiation  with  matter.  In  principle,  when  solved, 
they  will  give  us  all  the  information  we  need  to  analyze  the  problems 
herein  considered. 
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APPESDIX  S 


BELAII(»  EBTWEBl  THE  DEHSITT  MATRIX  EQUATIONS 
OF  I«3TION  AND  THE  RATE  EQUATIONS 

Consider  an  M-level  quantum  mechanical  system.  The  equations  of 
motion  for  this  system  In  terms  of  the  density  matrix  are 

n  ■  m: 

M 

liV  -  1«  >  (p. .  W.  -  p  W  .  ) 

^nn  l_^  kk  kn  '^nn  nk' 

k«l 

M 

*  I  ‘■kn  -  ‘■d. 

k-1 

n  y!  m: 

H 

liip  +  (E  -  E  )  p  +  ~  p  «  y  (3i  p  -  p  M  )  ,  (B.2) 

nm  '  m  n  nm  t  nm  J,  '  nk  km  nk  km  ^  / 

nm 


where  the  matrices  are  taken  with  respect  to  the  elgenkets  of  the  iinper- 

turhed  Hamiltonian,  ,  (energy  representation).  The  quantity 

is  the  probability  per  unit  time  that  the  system  will  go  from  state  k 

to  n  due  to  relaxation  mechanisms.  The  quantities  t  m  j  are  the 

nm  mn 

transverse  relaxation  times  which  are  related  to  the  llnewldth  of  a 
Lorentzlan  line  by  t  >  l/jt  v  ,  where  ^  v  is  the  full  llnewldth 
at  half  Intensity. 

Let  us  assume  that  the  applied  perturbations  vary  sinusoidally  in 
time  and  that  no  more  than  one  of  these  freqiiencles  Is  near  any  allowed 
transition.  We  may  then  write  the  perturbation  In  the  form 


+  e 


-la>  ,  t 
nk 


) 
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(B.3) 


(1) 

tihere  and' 

(B.4) 

Under  these  assunptioas  the  solutions  to  (B.2)  to  first  order  In  the 
perturbation  are 


,  iott  t 

W  (P  -  p  )  e  ® 


nm  urn  nn' 


nm 


ft  (0  -  a>  )  +  -i- 

'  mn  mn'  t 


(B.5) 


mn 


Substituting  (B.5)  into  (B.l)  gives 


M 

^nn  ^  ^  ^^nn  ^nk  ”  ^kk 
k«l 


Wkn)  -  I- 


^‘^kk  "  ^nn^ 


2t 


nk 


k«l 


1  +  V  *  “ 

nk'  nk  nk' 


(B.6) 


By  defining  the  quantities 


in 


2t 


nk 


nk 


1  +  (n  -  CD  )‘ 

nk  '  nk  nk' 


(B.T) 


the  equations  for  the  diagonal  coiis>onents  become 

M 

p  +  )  (p  w  .  -  p. ,  w.  )  +  >  r  .  (p  -  p. . )  »  0  ,  (b.8) 

nn  nn  nk  kk  kn'  /,  nk  '  nn  kk'  ’  '  ' 

k-1  k=l 

vhere  the  quantities  represent  the  transition  probabilities  due  to 

the  applied  radiation  fields. 

^^^By  convention  ve  have  and  cd^^  ■  -cUj^  . 
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If  ve  ccoQSlder  a  sjrstaB  eoatainlng  R  Identical  Independent 
syateas,  then  the  populatlDC  density  or  the  number  of  systems  In  the 
level  1  Is  given  Toy 

“l  ■  ***11  • 

Multiplying  (B.8)  by  N  and  using  (B.9)  gives 
M  M 

^  ^  I  <«1  ^'ik  -  \  '^kl)  ^  Z  "ik  (-1  -  \)  -  0  »  (B.IO) 

k-1  k-1 

2 

vhlch  Is  recognized  as  the  familiar  rate  equation. 

Thus  the  rate  equations  are  equivalent  to  the  diagonal  equations  of 
the  density  matrix. 
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AFEENDIX  C 


UERIVATION  OF  JAXMES'  NEOCLASSICAL  EQUATIONS  PR(»4  THE  DENSITY  MATRIX 

In  showing  the  equivalence  between  Jaynes'  semlclassical  equatlons^^ 
and  formulation  in  terms  of  the  density  matrix  we  shall  begin  with  the 
equations  of  motion  for  a  two-level  quantum  system  neglecting  the  effects 
of  relaxation  terms.  Let  the  energy  levels  be  eis  shown  in  Fig.  C.l, 


0  .  I  (s^  -  E^) 


FIG.  C.l — Energy  level  diagram  for  a  two-level  system. 

and  take  the  form  of  the  Hamiltonian  to  be 

^12  =  ® 


In  the  energy  representation  the  equations  of  motion  become 

Pii  =  -  S  -  P12) 

P22  ^  "  ii  ^Pl2  "  ^21^  (,c.2b) 
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and 


p^2  -  i  n  Pl2  =  -  ^  (Pgg  -  p^)  (C.2c) 

Pgl  +  1  fl  Pgl  «  -  ^  (Pll  -  pgg)  .  (C.2d) 

The  expectation  valixe  of  the  dipole  moment  M  is  given  by 
**  ^  ^^op^  ^  ^  ^*^12 

where  p  =  const  =  p^g  =  Pg^  .  The  energy  of  the  unperturbed  Hamiltonian 
relative  to  a  zero  of  energy  lying  midway  between  the  levels  is  denoted  by 
W  and  is  given  by 


W  =  =  I  h  n  (P22  -  Pll)  ”  (c-^) 

The  first  derivatives  of  these  quantities  are  found  to  be 

M  =  H  (Pi2  ^21^ 

=  i  n  p  (p^2  -  Pg^)  ,  (C.5) 

and 

W  =  I  ft  n  (pgg  - 

=  i  E  0  p  (p^2  -  p23^)  ,  (C.6) 

where  we  have  used  the  eq'uatlons  of  motion  (C.2)  to  eliminate  the  time 
derivatives.  Prom  (C.5)  and  (C.6)  we  may  derive  the  first  of  Jaynes' 
equations : 


W  =  E  M 


(C.7) 
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Taking  the  time  derivative  of  (C.5),  we  have 

M  =  i  n  (p^2  - 

=  -  p  (Pl2  P21)  ■"  (%)  ^  (Pll  ■  P22)  ' 

and  by  substituting  from  (C.3)  and  (C.4)  we  obtain 


These  two  equations,  (C.7)  and  (C.8),  constitute  two  of  Jaynes '  neoclassical 
equations.  The  first,  (C.7)/  states  that  the  rate  of  increase  in  the 
energy  of  the  quantum  system,  W  ,  is  equal  to  the  product  of  a  field, 

E  ,  and  a  current  M  .  The  second,  (C.8),  describes  the  time  develop¬ 
ment  of  the  dipole  moment,  M  ,  which  is  seen  to  be  the  equation  of  an 
undamped  harmonic  oscillator  driven  by  the  field  E  through  a  coupling 
constant  proportional  to  W  .  These  form  a  set  of  coupled  nonlinear 
differential  equations  and  have  been  proposed  by  Jaynes  as  a  method  of 
studying  the  interaction  of  radiation  with  matter. 

The  third  equation  is  simply  a  description  of  how  the  resultant 
dipole  moment  M  generates  the  fields  E  and  may  be  written  in  many 
forms.  The  one  used  in  this  paper  is  simply  to  define  a  complex  cavity 
susceptibility  (X^  =  -  IX^”)  and  then  to  relate  the  dipole  moment 

M  to  the  field  E  by 


E  =  ^M  .  (C.9) 

c 

We  may  phenomenologically  add  relaxation  terms  to  (C.7)  and  (C.8)  in 
the  following  way.  First  consider  (C.7)  which  describes  the  energy  of 
the  system,  and  let  the  characteristic  time  in  which  it  can  exchange 
energy  with  the  "lattice,"  or  thermalize,  be  T^^  .  Then  in  the  absence 

of  a  perturbation,  E  =  0  ,  the  system  will  return  to  thermal  equilibrium 
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if  ve  add  a  term  (W  -  W*)/tj^  ,  giving 


W  + 


W  -  W* 


-EM 


(C.IO) 


The  equation  for  M  ,  (P*8)>  nay  be  altered  to  account  for  relaxation 
mechanisms  if  ve  add  a  term  SW/Tg  in  the  form  of  a  damping  term.  This 
predicts  a  decay  time  of  Tg  for  the  dipole  moment.  The  modified  equation 
is 


(C.ll) 


If  ve  add  the  same  terms  to  the  density  matrix  ve  get  the  modified 
equations  including  relaxation 


(“u  -  “’ll' 


“u* 


ifi 


(p2i  ’  Pi2^ 


(C.12a) 


Pl2  -  i  n  p^2  Tg  ^12  "  ■  «  ^^22  "  ^11^  '  (C.l2b) 


where  we  have  included  only  one  diagonal  and  one  off-diagonal  equation. 
From  these  equations  we  find  the  following  derived  equations 


W  +  =  E(M  +  i  M) 

ll  ig 


(C.13) 

(C.l4) 


Upon  ccanparing  these  with  (C.IO)  and  (C.ll),  we  see  that  they  differ 

slightly  in  that  they  have  additional  terms  proportional  to  l/Tg  and 

(l/Tg)  .  This  is  not  unexpected  since  the  terms  were  added  phenomenologically 
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and  cannot  be  expected  to  be  valid  when  their  effect  is  large.  If  Tg 

is  not  too  small,  we  see  that  the  additional  terms  are  negligible  since 

11  -Q 

the  natural  frequencies  ^  10  arid  Tg  is  typically  10  .  By 
neglecting  small  quantities  we  see  that  (C.13)  and  (C.li^)  are  the  same 
as  (C.IO)  and  (C.ll). 

In  studying  the  interaction  of  radiation  with  a  two-level  system 
we  can  either  use  the  density  matrix  equations  (C.12)  or  Jaynes'  equations 
with  relaxation  added  (C.IO)  and  (C.ll). 

72 

It  should  also  be  mentioned  in  passing  that  Feynman,  et  al.,  have 
shown  that  for  the  two-level  system  one  can  write  the  equations  of  motion 
in  general  as; 


where 


and 


and 


r  =  0)  X  r 

^  ^  ^ 

} 

(C.15) 

£  =  (j*!  ,  rg  ,  r^) 

9 

(C.i6) 

ci  =  (cD^  ,  a>2  ,  05^) 

9 

(C.17) 

ri  =  p^g  +  Pg^ 

rg  H  i(p2^  .  p^2) 

(C.l8) 

^3  ~  ^22  '  ^11 

> 

“l  =  ^^12  ^21'/® 

0^2  =  ^(Vg^  -  Vi2)/« 

(C.19) 

=  9. 


No  use  will  be  specifically  made  of  these  equations  in  the  form  (C.15); 
they  are  merely  mentioned  for  completeness. 
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APPEHDIX  D 

DERIVATION  OF  THE  ABSORPTION  COEFFICIENT 

If  we  have  an  electromagnetic  field  propagating  with  the  velocity 
c  ,  the  energy  flow  is  given  by  the  Poynting  vector 


S  =  A  E  X  H 
A.  4n  ^  ~ 


which  for  E  perpendicular  to  H  is 


Isl 


(D.l) 


(D.2) 


Performing  a  time  average  res\ilts  in 


where  we  have  assumed  E  =  E^  cos  ost 

If  we  have  a  sample  which  has  a  resonance  at  this  frequency,  the 
density  of  power  absorbed  in  the  absence  of  saturation  is 


T, 


P  = 


2fi 


2  t,2 

ji  E  o 

—  (ergs/cm"^) 


(D.4) 


Considering  a  volume  1  cm  square  and  dz  long,  the  power  absorbed  is 


dP 


e  2  2 

oilA®  Tg  Eq 

2fi 


dz  (ergs) 


(D.5) 
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The  absorption  coefficient  is  defi^'i'ni  as 


y  = 


where -dP/dK  is  the  power  absorbed  per  unit  length.  Using  (D.3) 
in  (D.6)  gives^^^ 

(X)  lut  NA*  Tg 

y  =  -  -  , 

c  ft 


or 


7 


05  _1_ 


where  Q  is  the  sample  Q  defined  in  Chapter  III. 
s 


f 


^^^See  Gordy,  et  al.,^^  Eqs. 


(I^-.IO),  (^.l4),  and  (4.15a). 


(d.6) 

and  (D.5) 

(D.7) 

(D.8) 
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